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Uber die Anwendung der Plattentheorie zur Bestimmung 
von Warmespannungsfeldern | 


E. Tremmel, Salzburg 


Zusammenfassung. Es werden die in prismatischen Kérpern mit einfach zusammenhangenden 
Querschnitten und freien Mantelflachen von ebenen Temperaturfeldern hervorgerufenen Spannungs- 
zustande untersucht. Die den Randbedingungen geniigende Airysche Spannungsfunktion ent- 
spricht der Biegeflache einer dem Querschnitt kongruenten, unverdrehbar gelagerten Platte, 
auf die eine der Laplace-Ableitung des Temperaturfeldes proportionale Belastung wirkt. 


Fir einen prismatischen oder zylinderformigen, in der Richtung z erstreckten 
Kérper mit einfach zusammenhingendem Querschnitt sei ein zweidimensionales 
instationares Temperaturfeld #(x, y, t) durch die Oberflachen- und Anfangsbedingungen 
sowie gegebene zeitabhingige Warmequellenverteilungen W(z, y, t) festgelegt. Aus 
der Lésung der Fourierschen Warmeleitungsgleichung 


(1) 


in der mit a und A Temperaturleitzahl bzw. Warmeleitfaihigkeit bezeichnet sind, folgt 
bei Beriicksichtigung der mit den beliebigen Konstanten c,, ¢,, c; allgemein ange- 
schriebenen Oberflichenbedingungen 


“O(a, y,t 
6, (6, t) +n | | = 05 Ae, 8) (2) 
wobei s die von einem beliebigen Punkt aus geziihlte Bogenlange und n die Normale 
der Randkurve bedeuten, sowie der Anfangsbedingungen 


(=OF PSO a, Y) (3) 


fiir das gesuchte Temperaturfeld 
9 = O(a, y, 1). (4) 


Zur Bestimmung des von dieser Temperaturverteilung hervorgerufenen Spannungs- 
zustandes ermitteln wir nach dem von E. Melan und H. Parkus! angegebenen Ver- 
fahren, das nachstehend kurz beschrieben werden soll, zunichst das thermisch- 
elastische Verschiebungspotential ®; seine Ableitungen nach den Koordinaten- 
richtungen sind Lésungen der auf Temperaturwirkungen erweiterten elastischen 


Grundgleichungen 
este 2(1+) a0 


A neon ae Be (5) 
es gilt somit 
Lane (6) 
Ou, be 


1 E. Melan-H. Parkus: Warmespannungen, 8. 6ff. Wien. 1953. 
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wobei die Verschiebungskomponenten parallel zu den Achsen «x; des kartesischen 
Koordinatensystems in tiblicher Weise mit wu, und die kubische Dilatation mit 


3) 
Ou; 
$ =e Ox; 
bezeichnet seien, ferner » bzw. « Poissonsche Zahl bzw. Temperaturausdehnungs- 
koeffizienten bedeuten. Die Einfiihrung von (6) in (5) liefert die Differentialgleichung 


am lt+p 
AD 5 &O, (7) 
aus der das thermisch-elastische Verschiebungspotential bestimmt werden kann. 
Fir die Spannungen erhalt man aus (6) unter Beriicksichtigung der Verzerrungs- 


Verschiebungsgleichungen 


1 / eu, are, 
Coe 9 ee = oD (8) 
sowie des Hookeschen Gesetzes 
; I 
oun = 20 lee + qeyy a 8) ba, (9) 
worin das Kronecker-Symbol mit 6,;, bezeichnet ist (5:0— ; ote 7 schlieBlich 
CD 
oi, = + 26(q Ge — AP xn). (10) 


Die gemaB (10) bestimmten Spannungen werden aber im allgemeinen den Randbedin- 
gungen nicht geniigen. Fordern wir das Verschwinden der Normal- und Schub- 
spannungen o,,, und o,, an der freien Oberflache, in unserem Falle demnach 


Onn (8) = 9, (11) 

o.n (8) =0 (12) 
lings der Randkurve s, dann haben wir eine Airysche Spannungsfunktion 

F226 0 (13) 
mit 

AAF = 0 (14) 


zu suchen, welche die zur Tilgung der aus (10) abgeleiteten Randwerte o,,,, und 6, » 
erforderlichen Spannungen 


= — Onn: (15) 
on capa et Crm (16) 
liefert. Werden die aus (13) folgenden Spannungen 
= PQ 
Gip 26(5 = 426,,) (17) 


den aus dem thermisch-elastischen Verschiebungspotential mit (10) abgeleiteten 
Werten iiberlagert, so ist der endgiiltige, die Randbedingungen erfiillende Temperatur- 
spannungszustand durch 


Orn = Cin +65, = 2G| 
gegeben, wofiir wir mit 


 (Q—- ®) 
Ox, OX), 


—A(Q—9)6,,| (18) 


VO (19) 
ktirzer 
OV 
on = — 26 ( in guy 4 dix) (20) 
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schreiben wollen. y ist daher eine Spannungsfunktion, die — wie aus (19) unter 
Beachtung von (14) sowie der aus (7) abgeleiteten Differentialgleichung 
l+y 
AAG = =—— « Ad (21) 
hervorgeht — der inhomogenen Bipotentialgleichung 
AAY = —= +" x A8 (22) 
—y 


mit den lings des freien Randes geltenden (11) und (12) gleichwertigen Grenzbe- 
dingungen 


Y—0, (23) 
ow 
geniigt. 
Unter Verwendung von (1) erhalt man fiir 40 den Ausdruck 
1 od W 
Aj = a ot Het (25) 


in dem die Laplacesche Ableitung des Temperaturfeldes durch den _partiellen 
Differentialquotienten nach der Zeit sowie die Warmequellenverteilung ausgedriickt 
wird, so daB (22) auch auf die Form 


AAW = 7 > 9 (-2 4 4) 


a oa ria A (26) 
gebracht werden kann. 

Betrachten wir nun eine diinne, senkrecht zu ihrer Mittelebene belastete Platte 
von der gleichen Form wie der Querschnitt des hier untersuchten Korpers, die lings 
ihres Randes unverdrehbar gelagert ist. Nach der Kirchhoffschen Theorie besteht 
zwischen den auf die Breiteneinheit bezogenen Biege- und Torsionsmomenten m;,, 


und den Verschiebungen w(x, y,t) der Mittelfliche die Beziehung 


ew . 9 
Mm; ip N a, 02}, (1 v) mes A WwW 6; kj} o> (27) 
é 5 Bot. i h® 
(i = k... Biegemomente, 7 + k... Torsionsmomente); hierin bedeutet N = | Sey 


die Plattensteifigkeit. Setzen wir die Querzahl des im iibrigen beliebigen Platten- 
materials, das also keineswegs mit dem Material des betrachteten Korpers iiberein- 
stimmen mu8, gleich Null — itber den Wert von N wird noch verfiigt werden —, 
dann geht (27) in den einfachen Ausdruck 


Mik N= % (28) 


tuber. 
Die Verschiebungen sind durch die Plattengleichung 


1 
A Aw = P(e ¥ 1) (29) 


— hiebei entspricht p der je Flicheneinheit wirkenden, zeitlich veranderlich ange- 
nommenen Belastung — sowie die Randbedingungen bestimmt. Zufolge der voraus- 
gesetzten Hinspannung gilt langs des Plattenumfanges s: 


w= 0 (30) 
sowie 
Bie (31) 


12* 
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woraus mit (28) das Verschwinden der tangential an der Randkurve wirkenden 
Biegemomente m,, und der Torsionsmomente m,,, folgt, mithin auch 

m,, (8) = 0 (32). 
und 

Msn (s8) = 9 (33) 
gelten mu. 

Aus dem Vergleich von (22) mit (29) sowie von (23) baw. (24) mit (30) bzw. (31) 
ist nun zu ersehen, daB die Differentialgleichungen der Spannungsfunktion y und 
der Biegeflache w gleichartig aufgebaut sind und die diesen Funktionen vorgeschrie- 
benen Randbedingungen tbereinstimmen. 

Erstreckt sich diese Ubereinstimmung auch auf die rechten Seiten der Differential- 
gleichungen (22) und (29), dann gehen die Funktionen y und w bis auf einen zur Er- 
haltung der Dimensionsrichtigkeit erforderlichen konstanten Faktor C ineinander 
iiber. Setzen wir also mit (22) und (29) 


1+ C 
ae fees a Ad =~ p(x, Yy; t), (34) 
mit anderen Worten: bringen wir auf die Platte eine fiktive, der Laplaceschen Ableitung 
des gegebenen Temperaturfeldes proportionale Belastung auf, dann kann die 
Spannungsfunktion mit 


Y= Cw (35) 


durch die Biegefliche ausgedriickt werden. Da diese Ergebnisse offensichtlich auch 
fiir stationaire Temperaturverteilungen gelten, konnen wir behaupten: 

Steht ein prismenformiger Korper mit einfach zusammenhangendem Querschnitt 
und freier Mantelfliiche unter dem Einflu8 eines in der Langsrichtung unveranderlichen, 
sonst beliebigen stationéren oder instationiren Temperaturfeldes mit oder ohne 
Warmequellen, dann ist der von diesem Feld hervorgerufene, den Randbedingungen 
geniigende Spannungszustand durch die Biegefliche einer gedachten, lings des Quer- 
schnittrandes des Kérpers eingespannten elastischen Platte, auf die eine der Laplace- 
schen Ableitung des Temperaturfeldes proportionale Belastung wirkt, vollstandig 
bestimmt. 

Zur Vereinfachung der praktischen Berechnung setzen wir fiir die Plattensteifig- 
keit N, die die Dimension eines Momentes, also Kraft x Linge hat, unter Einfiihrung 
der Konstanten C,, die dem Kinheitsmoment entsprechen mége, 


I1—y | 
INS Ss Lee (36) 
— unter » ist die Querzahl des behandelten Korpers zu verstehen — und wahlen 


die in (35) auftretende Konstante C gleich der Langeneinheit ; damit erhalt man aus (34) 
die GroBe der aufzubringenden Last dimensionsrichtig, das hei8t als Kraft je Flachen- 
einheit, zu 


pie Ad. (37) 


Mit dieser Belastung also kénnte die Biegeflache w(«, y, t) bestimmt, aus ihr gemaB (35) 
die Spannungsfunktion (a, y,¢) und schlieBlich nach (20) das Spannungsfeld be- 
rechnet werden. | 

Praktisch ware mit diesem Verfahren nichts gewonnen, da das Aufsuchen einer 
die Randbedingungen (30) und (31) befriedigenden Lésung von (29) die gleichen 
Schwierigkeiten mit sich bringt, wie die Bestimmung einer Spannungsfunktion F 
gemiB (13), die die zur Tilgung der Randwerte o,,, (s) und 7,, (s) erforderlichen 
Spannungen o,,,,(s) und 7,» (s) liefert. 
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Abgesehen davon, da8 das hier erérterte Plattengleichnis die experimentelle Ab- 
schatzung derartiger Temperaturspannungszustinde gestattet, besteht der Vorteil 
seiner Anwendung vor allem in der Moglichkeit, die zahlreichen zur Berechnung 
querbelasteter elastischer Platten entwickelten Verfahren? und Hilfsmittel bei der 
Lésung ebener Wirmespannungsprobleme zu verwerten. 

Die fiir einzelne Plattenformen, wie Rechteck, Kreis usw., berechneten EinfluB- 
felder der wichtigsten Schnittgré8en — wir denken hier etwa an das von A. Pucher? 
verfaBte Werk — erlauben die rasche Bestimmung der Momente fiir eine beliebige 
Belastung. Handelt es sich dagegen um besondere Querschnitte, fiir die keine gleich- 
wertige Plattenlosung existiert, dann kann zur Ermittlung der Schnittmomente 
gegebenenfalls das zur Untersuchung von Gewdélbestaumauern entwickelte Versuchs- 
lastverfahren (Trial load method‘) angewandt werden. Die Platte wird hierbei durch 
zwei einander rechtwinkelig kreuzende Scharen liickenlos aneinander liegender, 
beidseitig eingespannter Trager von der Breite 1 ersetzt gedacht und die durch das 
Aufbringen der Belastung auf eine der beiden voneinander unabhangig angenommenen 
Tragerscharen gestérte Kontinuitit durch passend gewahlte Last- und Torsions- 
momentenverteilungen wieder hergestellt. Auch hier laBt sich der Rechenaufwand 
bei ausreichender Genauigkeit etwa durch Verwendung von EinfluBlinien gegenitiber 
der im allgemeinen auBerst miithsamen unmittelbaren Lésung des Scheibenproblems, 
fiir das kaum Hilfsmittel zur Verfiigung stehen, wesentlich verringern. 

Auf die Méglichkeit, ebene Spannungsaufgaben durch Zuriickfiihrung auf Platten- 
probleme zu losen, hat A. Timpe® schon vor mehr als fiinfzig Jahren hingewiesen. 
In den von ihm und anderen Autoren, wie K. Wieghardt®, H. M. Westergaard’, 
H. Cranz® untersuchten Fallen handelt es sich allerdings nicht um Temperatur- 
spannungszustinde, sondern um die Beanspruchung verschiedener unter auBerer 
Belastung stehender Querschnitte. Am belasteten Rand sind hier also die den 
auf die Scheibe wirkenden Normal- und Schubspannungen entsprechenden Ab- 
leitungen der Biegefliiche vorgegeben, wihrend die unbelasteten Bereiche des Platten- 
randes wieder als eingespannt anzusehen sind. _ 

So wurde z. B. der im Mittelquerschnitt des Boulderdams* vom Wasserdruck 
bewirkte Beanspruchungszustand mit Hilfe des Plattengleichnisses unter Anwendung 
des oben beschriebenen Versuchslastverfahrens berechnet. 

Untersuchungen von Warmespannungsfeldern in der hier geschilderten Art diirften 
dagegen bisher noch nicht vorgenommen worden sein’. 


Vgl. K. Girkmann: Flachentragwerke, 4. Aufl., S. 159ff. Wien. 1956. 

A. Pucher: EinfluBfelder elastischer Platten. Wien. 1951. 

Stress studies for Boulder Dam, Boulder Canyon Project, Final reports, Denver, Coll., 
- A, 1939; S. 119 ff. 

A. Timpe: Probleme der Spannungsverteilung in ebenen Systemen, 8S. 6. Leipzig. 1905. 
K. Wieghardt: Uber ein neues Verfahren ... Forschungsarb. Ingenieurwes. H. 49, 8. 15f. 
1908). 

é 2 M. Westergaard: Computation of Stresses... Publ. Roads 1—23 (March 1930). 

8 H. Cranz: Die experimentelle Bestimmung der Airyschen Spannungsfunktion mit Hilfe 
des Plattengleichnisses. Ingenieur-Arch. 10, S. 159ff. (1939). 

9 In einer dem Verfasser erst wibrend der Drucklegung der vorliegenden Arbeit bekannt- 
gewordenen Abhandlung von W. Novacki, Warschau (Abhandl. IVBH, Bd. 16, Zurich 1956), 
wird das durch stationaéire Warmequellen in Scheiben hervorgerufene thermisch-elastische 
Verschiebungspotential mit den Durchbiegungen einer diinnen, frei aufliegenden Platte in 
Beziehung gesetzt [s. unsere Gl. (21)]. Unter der Annahme, daB langs des Randes # = 0 gilt, 
verschwinden die Randnormalspannungen; zur Tilgung der Randschubspannungen Gl. (12) 
muB allerdings ein zusiitzlicher, aus einer Airyschen Spannungsfunktion abgeleiteter Spannungs- 
gustand eingefiihrt werden. Wahrend hier also der endgiiltige, den Randbedingungen genugende 
Spannungszustand aus den Momenten der unverdrehbar gelagerten Platte berechnet wird, ver- 
wendet W. Novacki die Plattenanalogie lediglich zur Bestimmung des thermisch-elastischen 


LBP 


om aR pw wo 
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Mit Hilfe der Plattenanalogie lat sich aber auch das durch Warmeleitungsgleichung 
und Randbedingungen bestimmte, in den Querschnitten eines prismatischen Korpers 
auftretende zweidimensionale Temperaturfeld statisch deuten: Die dem Querschnitt 
entsprechende, lings ihres Randes hier frei verdrehbar gelagerte Platte stellt sich 
unter einer der Laplaceschen Ableitung der Temperaturverteilung proportionalen 
Belastung sowie bestimmter von den Oberflachenbedingungen abhangiger Rand- 
momentenverteilungen so ein, daB das skalare, von der Momentensumme gebildete 
Feld bis auf einen Proportionalitatsfaktor und lineare Glieder mit dem Temperatur- 
feld iibereinstimmt. Fiir die als Momentensumme m= m,,+m,, bezeichnete 
skalare Invariante des Momententensors gilt nach H. Marcus die Differentialgleichung 


Am = — p; (a) 


nach (34) ist aber — p der positiven Laplaceschen Ableitung des Temperaturfeldes 
proportional, so da (a) bis auf einen konstanten Faktor in die Differentialgleichung 
der Wirmeleitung (1) tibergeht; die zur Erfiillung der Randbedingungen (2) erforder- 
lichen Randmomentenverteilungen, in denen wir m anstatt 9 zu schreiben haben, 
lassen die Differentialgleichung unverandert. Differentialgleichung und Randbe- 
dingungen von % und m stimmen daher iiberein, so daB zwischen der Momenten- 
summe und dem Temperaturfeld die Beziehung 


9=Am (b) 


besteht, wobei konstante und lineare Glieder auBer acht bleiben. 

Zur Vermeidung von Unklarheiten bei der Umrechnung der Schnittmomente m, ,, 
My, M,, auf die Scheibenspannungen o,,, 0,,, Tz, Soll nun die aus (35) mit (20) 
und (28) folgende Beziehung zwischen dem Momententensor der Platte und dem 
Spannungstensor der Scheibe aufgestellt werden. 

Wir fithren hierzu (35) in (20) ein und erhalten 


v1, = 260 (se — Awd), (38) 
wahrend fiir die Schnittmomente gemaB (28) 
min = —N ae (28a) 
gilt. Mit Hilfe des «-Tensors der Ebene", der durch 
le Whee 
Gp il fir] 1 £2 (39) 
lees LGn? y= ioe 


definiert ist, koOnnen wir die in (38) und (28a) rechts stehenden Differentialausdriicke 
miteinander verkniipfen, und zwar gilt — wenn wir die Indizes 7 und k in (28a) durch p 
und q ersetzen — 
ew Pw 
OF, Oa. ae Aw bin = Fi» Fea CO, OX,’ (40) 
SO dab die Scheibenspannungen durch die Schnittmomente in der allgemeinen, fiir 
beliebige orthogonale Koordinatensysteme geltenden Form 


2G C 
Cin =p Sivera™ na (41) 


Verschiebungspotentials und erfiillt die Randbedingungen nach den iiblichen Verfahren der 
Scheibentheorie. 

a 0 A. Duschek-A. Hochrainer: Grundziige der Tensorrechnung, II. Bd., 8. 209. Wien. 
950. 
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ausgedriickt werden kénnen. Setzen wir noch fiir NV aus (36) ein und fassen die in (41) 
rechts stehenden Faktoren zu einer neuen Konstanten 


- 2GC C l+y 
La) i 
Rae (Ce Te a) 
von der Dimension [/-?] zusammen, so folgt weiter 
O;,= Ke, p Eka "nq (43) 


und wir gelangen nach Ausfithrung der Tensormultiplikationen zu den speziellen 
Umrechnungsformeln 


Cog ST i Wy (44) 
Cyy = WU es (45) 
hgy So Mae Ke. (46) 


Die Normalspannungen sind demnach proportional den in ihre Richtung fallenden 
Momentenvektoren, wahrend die Schubspannungen durch die gleiche Proportionalitiats- 
konstante mit den negativen Vektoren der Torsionsmomente verkniipft sind. Die 
groBten Randspannungen o,, werden daher den gré8ten Einspannmomenten pro- 
portional sein. 


Betrachten wir noch kurz die durch (37) gegebene Belastungsfunktion. Wird die 
Laplacesche Ableitung des Temperaturfeldes mit (25) durch Zeitableitung und Warme- 
quellenverteilung ausgedriickt, so geht (37) tiber in 


Cin alaou: Ww 4 
ssa ae at A ee) 
Bei den in der Technik untersuchten thermischen Erscheinungen handelt es sich 
nun zumeist um Abklingungsvorginge; das Temperaturfeld ist dann — falls eine 


Produktdarstellung méglich ist — etwa in der Form 
O = O,*(x, y)e~4* + OA*(a, y) e~%! 


gegeben, wobei die beiden Terme #,* und #,* im allgemeinen durch mehrfache Reihen 
ausgedriickt sein werden. Im Falle einer Warmeabgabe an der Oberflaiche wird die 
Zeitableitung des Temperaturfeldes gegen die AuBenflache zu stairker abfallen, woraus 
im allgemeinen auf eine VergroBerung der Einspannmomente und damit auch der 
Randspannungen geschlossen werden kann. Punktformige Warmequellen bzw. 
-senken, also etwa in der Lingsrichtung des prismatischen Korpers angeordnete 
Heiz- oder Kiihlrohre, sind als Einzellasten, die auf das Flachenelement dx- dy der 
Platte wirken, aufzufassen; tiber endliche Bereiche wirkende Quellen oder Senken 
entsprechen demgema’ Lastverteilungen. Handelt es sich, wie z. B. beim Abbinden 
des Betons der annahernd prismatischen Baublocke von Talsperren, um Aufheizungs- 
vorginge zufolge gleichmaBig tiber den Querschnitt verteilter Warmequellen, dann 
wird die Platte neben der dem zeitlichen Temperaturgefalle ee proportionalen Last- 
verteilung durch eine der Quellintensitat je Flacheneinheit entsprechende Belastung 
beansprucht. Bei stationiren Aufheizungsvorgiingen dieser Art ist der Spannungs- 
zustand des Querschnittes daher allein durch die Biegeflache der unter Vollast stehen- 
den Platte bestimmt. 

In diesem Zusammenhang sei noch erwihnt, daB der von E. Melan" allgemein 
bewiesene Satz, nach dem ein ebenes Feld bei stationarer und quellenfreier Temperatur- 
verteilung spannungsfrei bleibt, mit Hilfe der Plattenanalogie fiir einfach zusammen- 


11 —. Melan: Temperaturverteilungen ohne Warmespannungen. Osterr. Ingenieur-Arch. 6, 
1 (1951). 
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hingende Bereiche offensichtlich statisch gedeutet werden kann. Bei stationarem 


Temperaturfeld, das heiBt bei - 0 reduziert sich die durch (37a) gegebene Be- 


t — 
lastung auf p= tga sind keine Quellen vorhanden, dann bleibt die Platte 


unbelastet und der Querschnitt des prismatischen Kérpers somit spannungsfrei. 


(Eingegangen am 2. April 1957) 


Die Kontrolle 
der Druckfestigkeit von Beton durch Stichproben. 1. Teil 


Von W. Eberl und G. SchneeweiB 
Mit 10 Textabbildungen 


Zusammenfassung. In dem vorliegenden 1. Teil der Arbeit wird die an Betonprobekérpern 
beobachtete Streuung der Druckfestigkeit mit Hilfe der Varianzanalyse in ihre auf verschiedene 
Ursachen zuriickgehenden Streuungsanteile aufgespalten. In einem bald folgenden 2. Teil 
werden die hier erzielten Ergebnisse fiir die Entwicklung eines stichprobenweisen Prifverfahrens 
von Beton nutzbar gemacht werden, das sowohl der technischen Forderung nach Sicherheit 
als auch dem wirtschaftlichen Verlangen nach méglichst geringen Herstellungs- und Priifkosten 
in gleicher Weise Rechnung tragt. Die Ergebnisse der Arbeit sollen eine mathematisch fundierte 
Stellungnahme zu den bestehenden deutschen und dsterreichischen Normen erméglichen. 


I. Einleitung 
A. Ubersicht 


In dieser Arbeit wird im Anschlu8 an+ untersucht, auf welche Weise und mit 
welchem Erfolg sich die Methoden der Stochastik (Wahrscheinlichkeitsrechnung und 
mathematische Statistik) auf die seit langem gebrauchliche stichprobenweise Kontrolle 
der Druckfestigkeit von Beton nutzbar machen lassen. Die Verfasser danken dem 
Vorstand der Technischen Versuchs- und Forschungsanstalt der Technischen Hoch- 
schule Wien, Herrn Prof. Dr. Slattenschek, bestens fiir die groBziigige und wohl- 
wollende Foérderung bei der Durchfiihrung dieser Untersuchung. 

Der vorliegende erste Teil behandelt die der Druckfestigkeit von Beton zugrunde 
liegenden Wahrscheinlichkeitsverteilungen. Er zeigt, wie man in Form einer Streuungs- 
zerlegung die an Probekérpern festgestellte Gesamtstreuung der Druckfestigkeit in 
ihre auf verschiedene Ursachen zuriickgehenden Streuungsanteile aufspalten kann. 
Im demniachst erscheinenden zweiten Teil werden diese Ergebnisse fiir eine die 
technischen und wirtschaftlichen Erfordernisse in gleicher Weise beriicksichtigende 
Methode der stichprobenweisen Betonpriifung verwertet werden. 

Vom technischen Standpunkt aus wird im strengsten Fall verlangt, da die Druck- 
festigkeit des Betons an keiner Stelle bestimmte Sollwerte unterschreitet. Das wird 
erreicht, wenn der Mittelwert der Druckfestigkeiten so gro8 ist, da8 trotz der Streuung 
der Druckfestigkeiten die Wahrscheinlichkeit einer Unterschreitung der Sollwerte 
mit hinreichender Sicherheit vernachlassigt werden kann. Einer groBen Streuung 
mu in diesem Fall durch einen entsprechend hohen Mittelwert Rechnung getragen 
werden. 

Aus wirtschaftlichen Griinden wird man hingegen versuchen, den Mittelwert der 
Druckfestigkeiten méglichst niedrig zu wihlen. Dies ist mit den technischen An- 
forderungen nur dann in Hinklang zu bringen, wenn man die Streuung der Druck- 


* E. Tremmel-A. Wogrin: Mathematisch-statistische Auswertung der Giitepriifungen von 
Massenbeton. Bauing. 30, 28—32 (1955). 
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festigkeiten klein halt. Die Grenzen, bis zu denen eine solche Verringerung der 


Streuung getrieben werden soll, bestimmen sich wiederum durch wirtschaftliche 
Erwagungen. 


B. Das Zahlenmaterial 


Die MeBwerte in Tab. 1 sind die Druckfestigkeiten in kg/cm? von Betonwiirfeln 
mit 20cm Kantenlinge. Sie wurden im Rahmen der laufend durchgefiihrten Giite- 
priifungen von Beton ermittelt und sind nach Betonwerken (W, bis W,) und Stich- 
proben geordnet. Eine Stichprobe umfa8t jeweils drei MeBwerte, die an Probekérpern 
aus derselben Mischung gewonnen wurden. Die Numerierung der Stichproben ent- 
spricht der zeitlichen Reihenfolge der Herstellung und Priifung der Mischungen, wobei 
von jedem Betonwerk etwa eine Mischung woéchentlich untersucht wurde. 


Tabelle 1. Druckfestigkeiten in kg/em? von Betonwirfeln mit 20cm Kantenlange 
von 7 Betonwerken 


Stich- 
probe 


f) 


Betonwerk 


Wy 


We 


Ws 


Wa 


Ws 


We 


W, 


OMA DAR WHY KE 


358 392 368 
324 307 308 
235 228 237 
317 346 276 
299 284 293 
330 376 381 
333 389 371 
333 334 364 
443 489 401 
434 354 366 
328 341 374 
279 302 320 
453 458 410 
261 279 244 
353 345 361 
301 402 379 
250 230 278 
335 342 300 
290 352 358 
239 349 315 
359 397 394 
324 336 352 
328 302 316 
285 285 303 
314 318 355 
271 245 209 
246 272 317 
322 386 328 
378 368 353 
419 344 355 
322 347 280 


251 218 246 
265 264 177 
193 206 240 
292 278 276 
281 291 321 
290 279 232 
361 320 398 


235 210 228: 


209 184 177 
167 176 205 
208 163 194 
272 252 232 
206 202 208 
265 279 262 
240 250 249 
231 239 264 
244 238 213 
262 281 273 
222 212 203 
207 167 225 


203 173 218 
334 260 271 
209 215 230 
242 292 252 
176 164 170 
211 198 210 
238 258 267 
239 209 223 
202 190 186 
220 186 228 
203 209 211 
189 174 190 
300 345 292 
255 259 250 
291 298 295 
167 162 163 
243 235 251 
226 215 225 
264 290 279 
206 189 190 
287 295 265 
236 245 239 
204 196 203 


250 256 239 


378 385 404 
389 324 334 
254 290 309 
204 218 201 
252 272 258 
362 351 350 
270 240 245 
322 370 324 
343 294 298 
261 252 245 
243 255 256 
268 271 278 
318 285 288 
306 326 305 
353 350 346 
274 253 277 
264 293 274 
282 277 285 
255 262 266 
233 264 280 
306 295 308 
234 214 224 
264 338 261 
206 239 239 
244 250 244 
250 242 274 
290 301 294 


298 318 316 
381 340 380 
406 400 427 
359 298 339 
399 363 409 
322 335 305 
379 333 322 
376 403 375 
353 307 282 
336 392 391 
361 361 321 
281 315 336 
264 249 256 
281 243 275 
366 317 410 
410 370 313 
384 318 336 
307 328 321 
472 422 466 
316 286 311 
250 234 221 
382 336 330 
322 303 337 
354 362 362 
336 305 321 
297 282 290 
381 361 335 
283 260 295 
353 358 382 


228 229 251 
226 207 213 
405 314 300 
118 105 113 
315 338 358 
215 195 206 
256 247 257 
183 223 264 
256 270 251 
196 239 236 
165 159 158 
262 224 256 
201 200 201 
180 169 164 
134 146 132 
141 164171 
226 216 198 


279 284 278 
268 268 256 
246 278 243 
176 176 178 
383 359 378 
210 229 242 
278 294 284 
246 240 272 
193 204 184 
135 147 152 
169 185 170 
320 358 344 
170 178 164 
199 214 196 
232 215 202 
138 138 146 
244 252 277 
276 304 286 
302 320 285 
316 290 293 
260 256 256 
184 183 165 


Allgemein wird die Druckfestigkeit des k-ten Probewiirfels aus der j-ten Stich- 
probe des i-ten Betonwerkes W; mit 2%;;;, bezeichnet. Stehen die Ergebnisse von 


a Betonwerken W,,.. 


., W, zur Verfiigung, so lauft ¢ von 1 bis a. Da nicht voraus- 


gesetzt wird, daB von jedem Betonwerk gleichviele Stichproben vorliegen, ist 7 = 1,..., b;, 
wo b, die Anzahl der Stichproben von W, bedeutet. Dagegen soll jede Stichprobe 
den gleichen Umfang c haben, also k = 1,...,c. In Tab. 1 ist a= 7, die 6, sind 
der Reihe nach 31, 20, 23, 28, 29, 17 und 22, und c = 3. 


12a 


174 W. Eberl und G. Schneewei8: 


C. Die Wahrscheinlichkeitsverteilungen 


Es wird nun die fiir die-Anwendbarkeit der Stochastik grundlegende Annahme 
gemacht, da die Zahlen von Tab. 1 Beobachtungen von Zufallsgr6Ben sind. Diese 
haben folgende Wahrscheinlichkeitsverteilungen (Wvn): 

1. Die Wvn der Druckfestigkeiten innerhalb der Mischungen, kurz 
Mischungsverteilungen. Fiir eine bestimmte Mischung besteht eine stetige Wv 
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Abb. 1. ——— Wvn der 2;,;, in fiinf Mischungen des z-ten Betonwerkes; ------- Wvn der y,; im 7-ten 
Betonwerk; —----- Wv der pw; 


fiir die Druckfestigkeit eines beliebig herausgegriffenen Probekérpers. Von der 
stochastischen Abhangigkeit der Druckfestigkeiten innerhalb einer Stichprobe wird 
abgesehen, indem die Menge der Mischung als groB angenommen wird. Der Mittel- 
wert dieser Wv heiBt Mischungsmittelwert, die Streuung Mischungsstreuung. Handelt 
es sich um die j-te Mischung des i-ten Betonwerkes W,, so wird dieser Mittelwert mit 
fiz, die Streuung mit o?; bezeichnet. 


2. Die Wvn der Mischungsmittelwerte innerhalb der Baustellen 
(eq Betonwerke). Fir jedes W, besitzen die in 1 eingefiithrten Mischungsmittelwerte 
44; eine stetige Wv mit dem Betonwerksmittelwert uw; und der Streuung o?. 


3. Die Wv der Betonwerksmittelwerte. Die in 2 beschriebenen Beton- 
werksmittelwerte 4; unterliegen einer stetigen Wv mit dem Gesamtmittelwert uw) und 
der Streuung o%. 
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Die Mittelwerte und Streuungen der unter 1, 2 und 3 angefiihrten Wvn sind un- 
bekannte Parameter von stetigen Wvn. Sie kénnen mit Hilfe geeigneter Funktionen 
der x;;, geschitzt, aber nicht berechnet werden und sind von ihren Schatzwerten 
begrifflich scharf zu trennen. 

Abb. 1 stellt einige dieser Wvn dar: Die fiinf durchgezogenen Kurven der oberen 
Reihe geben die Wahrscheinlichkeitsdichten der Druckfestigkeiten x,;,, 7 = 1,..., 5, 
fiir fiinf Mischungen des i-ten Betonwerkes wieder. Die strichlierte Kurve rechts 
deutet die Wahrscheinlichkeitsdichte der Mischungsmittelwerte faz von W, an. 
In entsprechender Weise enthalt die untere Reihe zu den Wahrscheinlichkeitsdichten 
der Mischungsmittelwerte fiir drei beliebige Betonwerke die strichpunktiert wieder- 
gegebene Wahrscheinlichkeitsdichte der Betonwerksmittelwerte «4, mit dem Gesamt- 
mittelwert [o. 

o7;, of und o? sind, roh gesprochen, MafBe fiir die Breite der entsprechenden 
Verteilungskurven, die natiirlich durchaus nicht symmetrisch sein miissen. 

05 besitzt eine geringe technische Bedeutung. 

In o? kommen alle Unterschiede in der Art, Menge und dem Zustand der Beton- 
komponenten, aber auch die von Mischung zu Mischung auftretenden UngleichmaBig- 
keiten in der Herstellung, Lagerung und Priifung der Probekérper zum Ausdruck. 

Durch of werden die Schwankungen der Druckfestigkeit zwischen Probekérpern 
derselben Mischung gekennzeichnet. 07, wird im wesentlichen von der GleichmaBig- 
keit der Durchmischung abhangen. 

Mit Schatzwerten fiir o? und of werden also auch Mafe fiir die Zielsicherheit 
der Betonzusammensetzung und die GleichmiéBigkeit der Betondurchmischung ge- 
wonnen. Die Ermittlung dieser Schatzwerte erfolgt mit Hilfe von Streuungszer- 
legungen (Varianzanalysen). 


D. Unabhingigkeit von der Zeit 


Hiner statistischen Auswertung des Zahlenmaterials von Tab. 1 muB die Unter- 
suchung vorangehen, ob sich die in C 1 und C 2 beschriebenen Wvn nicht etwa im 
Zeitraum, iiber den sich die Messungen erstrecken, in systematischer, das heiBt zeit- 
abhingiger Weise geindert haben. Drei Beispiele sollen zeigen, wie es zu solchen 
Anderungen kommen kann: 

1. Ein Betonwerk W,; hat seine veralteten Mischaggregate durch neue ersetzt 
und dadurch eine wesentliche Verringerung der Mischungsstreuungen 0}, 0%,... 
erzielt. 

2. Infolge mehrmaliger Beanstandung von Probewiirfeln wegen zu geringer Druck- 
festigkeit steigen die Mischungsmittelwerte j,, Mi2,... des Betonwerkes W, in 
mehr als zufalliger Weise an. 

3. Im Zuge eines groBangelegten Planes zur strafferen Standardisierung der Beton- 
herstellung wird die Streuung of zwischen den Betonwerksmittelwerten ,, 
fo,... wesentlich verringert. 

Die Priifung auf Zeitabhangigkeit der Wvn geschieht sehr einfach mit Hilfe von 
Lauftests?: i 

In einer beliebigen Permutation von n, Plus- und n, Minuszeichen, etwa 


re eh ligt sia anc ae aR ae en ean Pak ee ee 
nennt man jede nicht weiter ausdehnbare Aufeinanderfolge gleicher Vorzeichen einen 
Lauf. Seine Lange ist die Anzahl seiner Vorzeichen. Die Laufzahl L ist die Gesamt- 
zahl der Liufe. Im angefiihrten Beispiel ist n; = 6, n. = 10, L =7 und die Lauf- 


2 A. M. Mood: The Distribution Theory of Runs. Ann. Math. Stat. 11, 367—392 (1940). 
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langen sind der Reihe nach 2, 1, 3, 3, 3, 2, 2. Zu einer solchen Vorzeichenfolge 
gelangt man, wenn man in einer Mefreihe 2%, 2, ...,%, einer stetigen ZufallsgréBe z 
jeden Wert durch + oder — ersetzt, je nachdem er oberhalb oder unterhalb des 
Medianwertes Zmeq der MeBwerte® liegt. Ist eine ungerade Zahl, so wird die dem 
Medianwert entsprechende Beobachtung weggelassen, so da man auf jeden Fall 
gleichviele Plus- und Minuszeichen erhalt. Kin aus Plus- (Minus-) Zeichen bestehender 
Lauf heiBt Ober- (Unter-) lauf, weil die entsprechenden z-Werte oberhalb (unterhalb) 
VON Zmea liegen. 

Zu einer weiteren Vorzeichenfolge gelangt man, wenn man fiir die Beobachtungs- 
reihe z,,...,2, der Reihe nach die Vorzeichen der Differenzen z;,, — z, aufeinander 
folgender Werte (¢ = 1,...,2— 1) anschreibt. Nullen, die dabei auftreten, laBt 
man weg. Liiufe aus Plus- (Minus-) Zeichen entsprechen ansteigenden (abnehmenden) 
z-Werten und heiBen daher Aufwérts- (Abwérts-) léufe. 

Man setzt nun voraus, daB alle Beobachtungen unabhangig voneinander aus 
derselben Wv oder aus solchen, die sich nur in zufalliger Weise voneinander unter- 
scheiden, hervorgehen. Dann 148t sich mit Hilfe der Kombinatorik unter anderem 
die Wv der Laufzahl L ermitteln. Als Mittelwert It {ZL} und Streuung © {L} der 
Laufzahlen LZ erhailt man im ersten Fall 


Dy 
M{L,} = =F, 


n—l1 


Sh} =F(l-gaq)=7a 


und im zweiten Fall 


2n—1 
NL} = —, 
16n—29 
i) {L2} ce se 


Fiir » = 20 kann sowohl ZL, als auch LZ, mit guter Naherung als normalverteilt 
betrachtet und dementsprechend auf Signifikanz beurteilt werden. 

Entstammen aber 2, 2:,...,2, der Reihe nach » Wvn, deren Mittelwerte eine 
wachsende (fallende) Folge bilden, so werden die z-Werte im allgemeinen in mehr 
als zufalliger Weise zu-(ab-) nehmen. Hand in Hand damit werden lingere Liufe 
auftreten und die Laufzahl L wird klein ausfallen. Daher darf umgekehrt ein signifikant 
kleines L als Anzeichen einseitiger Verschiebungen der zugrunde liegenden Wv an- 
gesehen werden. 


c 
ay, i D> &;;,/¢ ist der Durchschnitt aus den Druckfestigkeiten der j-ten Stich- 
=1 
probe des 7-ten Betonwerkes, w,;; = max x;;, — min 2;;, (Differenz von groBter und 
k k 


kleinster Druckfestigkeit in einer Stichprobe) die zugehdrige Spannweite. Hine 


systematische Anderung der j,1,..., fay, Oder der oj,..., 03), wiirde eine entspre- 
chende der 2;,..., %ip,, bzw. der wy,..., W;,, zur Folge haben. Mit Hilfe von 
Lauftests der x;;, und w,;, 7 = 1,...,6;, lassen sich daher systematische Anderungen 


der Mischungsmittelwerte und -streuungen (Beispiele 1 und 2) und damit auch die 


* Die MeBreihe sei der GréBe nach geordnet: 24) < 2) <...< 2m)» (Wenn man die Werte 


einer stetigen ZufallsgroBe genau genug miBt, werden die MeBergebnisse im allgemeinen ver- 
schieden ausfallen.) Dann ist 


med = 2(m +1) fir 2 = 2m + 1 
und 


1 
med = Fy (%m) + 2m + 1) fiir n = 2m. 
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Veranderlichkeit der Wvn C 2 und © 1 nachweisen. Fiir die entsprechende Priifung 
der Wv C 3, auf die sich Beispiel 3 bezieht, hatte man eine groBere Anzahl von Beton- 
werken zu untersuchen. 

Tab. 2 gibt die Ergebnisse der Lauftests fiir die vorliegenden 7 Betonwerke 
wieder: nur in einem Fall liegt m knapp unter 20, sonst immer dariiber, also kann 
die Laufzahl durchwegs als normalverteilt angesehen werden. Da alle beobachteten 
Laufzahlen nur um weniger als das Doppelte der Standardabweichung hinter ihrem Mittel- 
wert zuriickbleiben, liefern die Lauftests weder einen Einwand gegen den stationaren 
Charakter der Wvn C 1, noch widerlegen sie das Vorhandensein der Wv C 2. 


Tabelle 2. Priifung der Wvn auf Abhangigkeit von der Zeit. 


Anzahl Ober- und Unterlaufe Aufwarts- und Abwartsliufe 
Beton- ” Get Beobachtete | ia Beobachtete | 
werk eee Laufzahl L, von _‘Mittelwert Asteria Betr L,von _ | Mittelwert aaeoet 
b; z., w.. | ® {La} V6 {I} & w m {Za} Vs 
is. ig 1 ij. | S {Ls} 
W, 31 12 16 le 2 16° 2:7 19 17 20°3 2°3 
W, 20 10 7 11:0 2°2 11 11 jee E320) 18 
Ws 23 13 8 12°5 2°3 18 16 few L5sO 19 
W, 28 16 15 15:0 2°6 15 16 | 18°3 2°2 
W; 29 20 17 15°5 2°6 19 20 | 19-0 2°2 
W, 17 8 13 9°5 2°0 13 15 Niner 8 BO, 16 
W, 22 10 16 12°0 2°3 11 1D, Nee 1-9 


E. Ubergang zu Logarithmen 


Infolge der grundlegenden Bedeutung, die der Normal- oder Gauf-Verteilung 
fiir einen GroBteil der heute iiblichen statistischen Priifverfahren zukommt, soll nun 
untersucht werden, ob die in C1, C 2, C3 
angeftthrten Wvn als Normalverteilungen 
angesehen werden kénnen oder sich auf ein- 
fache Weise in solche transformieren lassen. 
Linksschiefe Wvn, die im Versuchswesen 
haufig auftreten, konnen z. B. oft durch 
Logarithmieren der Zufallsgr6Be in Normal- 
verteilungen iibergefiihrt werden. 


A Modelwert 
4 Mittelwert 
v Medanwert 


A vl 

a) Form der Wvn. Abb. 2 zeigt sche- pb. 2. Beispiel einer linksschiefen Wv 
matisch eine (eingipfelige) linksschiefe Wv 
einer ZufallsgréBe z mit der Wahrscheinlichkeitsdichte f(z). Im Gegensatz zu den 
symmetrischen Wvn fallen Mittelwert ~, Medianwert » und Modalwert A nicht 


+co 


zusammen. Das dritte Moment \(z — u)* f(z) dz ist zum Unterschied von symmetri- 


schen Wvn von Null verschieden. 

Eine Linksschiefe der Wvn C2 kann mit Hilfe vereinfachender Annahmen 
plausibel gemacht werden: 

Beriicksichtigt man nur den Kinflu8 des Zement-Wasser-Wertes z/w, so bestimmt 
sich die Druckfestigkeit x nach der Formel von Bolomey 

eek! ky 2[w, 

wobei k, und k, Konstante sind. Die Wv von x wird also im wesentlichen durch die 
des Quotienten z/w bestimmt. 
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z und w mdgen unabhangig voneinander nach G(u,, o2) bzw. G(Ww, 0%) verteilt 
sein‘. In der Praxis sind o, und o,, klein gegentiber uw, bzw. w. Bezeichnet man daher 
die auf die Mittelwerte , und w,, bezogenen ZufallsgroBen mit Az und Aw, also 


Az=2— pf, AW=W— fy, 
so werden Az/u, und Awju, im allgemeinen klein ausfallen. 


Fiir den Mittelwert von z/w gilt: 


ed ee bd LT 


Reihenentwicklung in der geschwungenen Klammer und Veérnachlassigung aller 
hoheren Potenzen von Aw/u, als der zweiten liefert 


mts} = se (1 +e) 


SI 


AS 
SSs 


YM 
S 


Summenhautigkelispozelt 


Abb. 3. Summenhiufigkeiteh der Druckfestigkeiten x,;, der einzelnen Betonwerke im Wahrschein- 
lichkeitsnetz 


Durch eine einfache Transformation eines Doppelintegrals laBt sich zeigen, daB u,/u, der 
Medianwert von z/w ist. Somit ist fiir z/w der Medianwert kleiner als der Mittelwert 
und dasselbe gilt fiir x Ist also o2/u2, groB genug, so erhdlt man eine merkliche 
Linksschiefe der Wv von 2. 


Abb. 3 zeigt auf Wahrscheinlichkeitspapier und in Prozenten die aufsummierten 
Haufigkeiten der Druckfestigkeiten 2,;, ftir die einzelnen Betonwerke W, bis W,. 
Verlaufen die Streckenziige einigermaBen geradlinig, so konnen die a;,,, fiir jedes W; 
als Beobachtungen einer normalverteilten Zufallsgr6Be angesehen werden. Wie 
Abb. 3 zeigt, sind die Streckenziige zum gréBeren Teil nach oben durchgebogen, 


Eine Gau8-Verteilung mit dem Mittelwert « und der Streuung o? wird kurz mit G@(u, o?) 
bezeichnet. Genau genommen handelt es sich hier um beiderseitig verkiirzte GauSverteilungen, 
da die Wahrscheinlichkeit z. B. eines w-Wertes auBerhalb von [e, 2 uy —e] exakt 0 ist, wobei ¢ 


eine kleine positive Gré8e ist. Dadurch rechtfertigt sich auch die folgende gliedweise Integration 
einer Reihe. 
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lassen also auf schwach linksschiefe Wvn schlieBen. Es liegt daher die Frage nahe, 
ob nicht die y;;, = log #,;;, mit mehr Berechtigung als normalverteilt betrachtet 


werden k6nnen. 

In Tab. 3 sind die Logarithmen y;,,, der Druckfestigkeiten von Tab. 1 zusammen- 
gestellt. 

Abb. 4 stellt die aufsummierten Haufigkeiten der y,;, fiir die einzelnen Betonwerke 
dar. Ein Vergleich der Abb. 3 und 4 zeigt, da die Streckenziige in Abb. 4 zwar 
etwas geradliniger verlaufen, die Unterschiede jedoch nicht so groB sind, daB sie die 
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Abb. 4. Summenhiaufigkeiten der logarithmierten Druckfestigkeiten y,;, der einzelnen Betonwerke 
im Wahrscheinlichkeitsnetz 


Entscheidung zwischen der Verwendung der Urwerte w;;;, oder deren Logarithmen y;;;, 
herbeifiihren. Im folgenden wird jedoch eine Beziehung gezeigt, die diese Frage 
zugunsten der y;;;, entscheidet. 


b) Vorbereitungen. Zur Vereinfachung werden Summen der 2;;, bzw. y;;, mit 
den entsprechenden GrofSbuchstaben bezeichnet, und Zeiger, iiber die summiert 
wird, durch Punkte ersetzt. Fiir die zugehérigen Durchschnitte werden Kleinbuch- 
staben mit den Zeigern der Summen verwendet. Also: 


c 
BG 
Xi5, = LQ, Vise C45, = = | 
b5 ¢ 
SN aye Xi. 
tn eee Liste Li, = = (1) 
Seen qe 
a b; ¢c 
D.C) Sana x gh eee 
. —_ aok ane B 
¢=1j=1h=1 1¢ | 


a 
—- S 0 
wo B, = 2b: In genau entsprechender Weise gewinnt man ausgehend von den y;;, 


die Summen Y,,;,, Y;,., Y.., und die Durchschnitte y,, , Yio Yous 
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Im weiteren Verlauf wird sehr oft von folgenden Tatsachen Gebrauch gemacht>: 
OG) rig sca +n seien ” voneinander unabhingige ZufallsgroBen mit WM {z,;} = uw, 
Che 8d bes hae, ytd ip Og Belenian beliebige Konstante. Dann gilt 


OTe agent = ai, eS Ea, | (2) 
© {a 2, +... +,2,) =atot +... +a2 0%. j 


B) Durch die Spezialisierung uw, =...=u, =m, 0? =... =02 =o? und a, = = 
= a, = 1/n ergibt sich, wenn man (z, +... + 2,)/n =z setzt: 
MZB=n, SZ=o%n. (3) 


In Worten: Der Durchschnitt von n voneinander unabhangigen, mit dem Mittelwert 
und der Streuung o* verteilten ZufallsgréBen hat ebenfalls den Mittelwert mu, aber 
die kleinere Streuung o?/n. 

y) Sind insbesondere z,,...,z, unabhiingig voneinander nach G(u, 02) verteilt, 
so ist z nach G(u, o?/n) verteilt. 

In * wird gezeigt: 


6) Sind z,,...,2, ” voneinander unabhangige Beobachtungen einer nach G(w, 0?) 
verteilten Zufalisgr6Be z, so sind Durchschnitt z und Spannweite w (= max z, — 
— min z,) dieser » Werte unabhingig voneinander verteilt. be 

k 


c) Abhangigkeit von Mittelwerten und Streuungen. Fiir die im Abschnitt II folgende 
Streuungszerlegung ware es sehr vorteilhaft, wenn die Verteilungen C 1 als GauB- 
Verteilungen gleicher Streuung angenommen werden kénnten. Die Erfahrung lehrt 
jedoch, dafS die Druckfestigkeiten innerhalb einer Mischung im allgemeinen um so 
stirker streuen, je groBer ihr Mittelwert ist, und diese Tatsache kommt auch im 
vorliegenden Zahlenmaterial zum Ausdruck. Hingegen weisen die logarithmierten 
Druckfestigkeiten y,;;;,, GauB-Verteilungen auf, deren Streuungen fiir alle Mischungen 
als gleich angesehen werden diirfen. Das soll nun im einzelnen gezeigt werden. 

Berechnet man fiir jede Stichprobe des i-ten Betonwerkes W; sowohl den Durch- 
schnitt 2,; als auch die Spannweite w,;,7 = 1,...,0;, so erhalt man eine Folge von 
b; Zahlenpaaren | 

(i, War), - > <5 (Xeo,.. Wio,)s (4) 


die sich auf verschiedene Mischungen beziehen. Sie sind daher voneinander stochastisch 
unabhangig. Dies hat zur Folge, da die Wahrscheinlichkeitsdichte fiir das gleichzeitige 
Auftreten der Zahlenpaare (4) gleich dem Produkt der Wahrscheinlichkeitsdichten 
der einzelnen Zahlenpaare ist. 


Man legt nun die zu priifende (im vorliegenden Fall: zu widerlegende) Annahme 
oder Nullhypothese zugrunde, da& die x,,, innerhalb einer jeden Mischung von W;, 
mit derselben Streuung 07, = of = ...= 0%, = 07” normalverteilt sind. Da nach 


b 6 dann x,; und w,; unabhangig voneinander verteilt, also auch unkorreliert sind,’ 
so ist die Wahrscheinlichkeitsdichte eines jeden Zahlenpaares (%;;,, wW;;), 7 = 1,..., 0:, 
ihrerseits gleich dem Produkt der Wahrscheinlichkeitsdichten von 2,;, und w,;. 


Trifft die Nullhypothese zu, so hat jede Folge 
(Li ’ wi), Cera ti) (in,.. Win,» (5) 
die aus (4) durch Permutation der w,,,..., Wi», hervorgeht, dieselbe Wahrscheinlich- 
5 A. Hald: Statistical Theory with Engineering Applications & Statistical Tables and 
Formulas. New York: Wiley. 1952. 
6 J. F. Daly: On the Use of Sample Range in an Analogue of Student’s t-test. Ann. Math. 


Stat. 17, 71—74 (1946). 
Ingenieur-Archiv XI/3 13 
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keitsdichte wie die Folge (4), da die Spannweite infolge der iibereinstimmenden 

Streuungen in allen Mischungen dieselbe Dichte besitzt. Die Nullhypothese wird nun 

in einer ihrer Auswirkungen, naimlich der Unkorreliertheit der x;;, und w,; getestet. 
Zu jeder Folge (4) gehdrt ein Wert r des x-Korrelationskoeffizventen 


by 
AT a Ee ) (w; j see Cg ) 
j= 
fo aimee 5 7° 
oy (ig, 0) a we 
j=l j=1 
Durch Permutation der Spannweiten bekommt man insgesamt b;! = 6; (b; — 1)... 


3.2.1 gleichwahrscheinliche, im allgemeinen verschiedene Wecte von r. Zu einem 
entsprechend klein vorgegebenen « (etwa 0:05) kann man dann eine nichtnegative 
ganze Zahl M bestimmen, so da die Wahrscheinlichkeit, einen der M groBten oder 
der M kleinsten r-Werte zu erhalten, unter den getroffenen Annahmen héchstens 
den Wert « bekommt. Ist der (4) entsprechende r-Wert einer dieser extremen Werte, 
so hat man einen signifikant groBen bzw. kleinen Wert erhalten und mu8 ihn als 
statistischen Einwand gegen die Annahme unkorrelierter x-Durchschnitte und 
x-Spannweiten und damit als statistische Widerlegung der Nullhypothese werten. 

In’, II, S. 131f., wird gezeigt, daB dieser Test mit guter Naherung durch den 
allgemein iiblichen ¢-Test ersetzt werden darf. Danach wird als PriifgroBe der Wert 


He 
se V1i—r i oe 
mit b; — 2 Freiheitsgraden verwendet. 

Die Gleichheit von Streuungen lieBe sich auch mit Hilfe eines Bartlett-Tests 
priifen, jedoch ist dieser in diesem besonderen Fall zu unempfindlich. Statt mit den 
Spannweiten hatte der Test auch mit den Streuungsschitzwerten durchgefiihrt werden 
konnen, es wurden aber im Hinblick auf die einfachere Berechnung die Spannweiten 
gewahlt. Der dadurch verlorengehende Teil der Information ist nicht groB. 

Tab. 4 enthalt die Korrelationskoeffizienten und Priifgré8en fiir die 2;;, und y;;,. 
Die durchwegs positiven Vorzeichen der x-Korrelationskoeffizienten (Spalte 2) und 
die beiden héchst signifikanten t-Werte (Spalte 4) deuten auf eine positive Korrelation 
zwischen den aus den «,;, abgeleiteten Durchschnitten und Spannweiten hin. Hin- 
gegen zeigt das Abwechseln der Vorzeichen bei den y-Korrelationskoeffizienten 
(Spalte 3) zusammen mit dem Fehlen signifikanter t-Werte (Spalte 5), daB man durch 
Ubergang zu den Logarithmen zu unkorrelierten Durchschnitten und Spannweiten 
gelangt. 


t 


Tabelle 4. Korrelationskoeffizienten von Stichprobendurchschnitt und -spann- 
weite fur die Urwerte (#;;,) und ihre Logarithmen (Yian) 


Korrelationskoeffizient r der Priifgr6Be t von 7 der Zahl der 
Betonwerk Basy Vein tie age Freiheitsgrade 
WwW, + 0°254 — 0°014 + 1°41 — 0:08 29 
W, + 0-211 — 0°203 + 0°92 — 0°88 18 
W, + 0°485 + 0°256 + 2°54** + 1:21 21 
W, + 0161 — 0:088 + 0°83 — 019 26 
W; + 0°294 + 0071 + 1:60 + 0°37 27 
W, + 0°595 + 0°249 + 2°87** -+- 0:99 15 
W, + 0°387 — 0:078 + 1:88 — 0°35 20 


7M. G. Kendall: The Advanced Theory of Statistics, 2 Bde. London: Griffin. 1948. 
** Hochst signifikant. (« = 0°01). 
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Diese Ergebnisse bestiitigen die bereits erwihnte Tatsache, daB Mischungen mit 
héherer mittlerer Druckfestigkeit im allgemeinen auch eine groBere Streuung der 
Druckfestigkeit aufweisen. Fiir die y,,;, ist dies nicht der Fall. Der Streuungszerlegung 
werden daher die Logarithmen der Druckfestigkeiten zugrunde gelegt. Da diese 
nach Abb. 4 auch eine geringere Abweichung von der Normalverteilung als die Ur- 
werte zeigen, ist ein weiterer Vorteil. 


II. Streuungszerlegung 


Die folgende Streuungszerlegung hat den Zweck, die Gesamtstreuung 6? der 
Yisx = log %,;, in drei Anteile aufzuspalten: 

1. Streuung o? der y,;, innerhalb der Mischungen; 

2. Streuung of der y-Mischungsmittelwerte j,; desselben Betonwerkes; 

3. Streuung o2 der y-Betonwerksmittelwerte j,. 
Die y-Wvn entsprechen den in I C angefiihrten x-Wvn. 


Die y;;, sind in Tab. 3 in gleicher Weise wie die MeBwerte x,;, in Tab. 1 zu- 
sammengestellt. 
A. Modell 


Das der Auswertung zugrunde liegende Modell besteht in der Annahme, da sich 

die y,;, in der Form 

Yair = + Oe + Biz + Yizn (6) 
darstellen lassen, wobei w ein unbekannter Parameter ist, und die anderen drei 
Summanden auf der rechten Seite von (6) unabhangig voneinander nach G(0, o2), 
G(0, of) bzw. G(0, o?) verteilt sind. 

Die Bedeutung dieser Zerlegung ergibt sich aus folgenden Feststellungen: 

Die zu einer bestimmten Mischung M, gehorigen y;,,;, besitzen eine Wv mit dem 
Mittelwert “,,; =" + «;, + 8, Die Wv dieser y-Mischungsmittelwerte w,; fiir ein 
Betonwerk W, hat den Mittelwert u;= + «,. Der Parameter uw ist der Mittelwert 
der Wv der y-Betonwerksmittelwerte (;. 

Man erhalt daher formelmiaBig die y,;,, 6,; und «;, ndem man die y;;,, “;; und pu; 
auf ihre Mittelwerte bezieht, das heiBt von ihrem Mittelwert aus miBt: 


Vise = Yin — Mas 
Bis = Bas — bi 
Os = hi — Uh 
Da die «,, B,; und y,, als voneinander unabhangig vorausgesetzt werden, gilt 
fiir die Streuung (engl. variance) o? der y;;, 
Go? = o2 + a} +8. 


Daher heiBen o%, of und o% Varianzkomponenten. 


B. Voraussetzungen 


Die Annahme, da sich die y,;, als Summen eines Parameters und dreier unab- 
hingig voneinander normalverteilter Zufallsgro8en darstellen lassen, bedarf einer 
Priifung. Diese erfolgt schrittweise und wird eine Aufspaltung des Zahlenmaterials 
notwendig machen. 

Zunaichst ist es unwabrscheinlich, daB die y;; unabhingig vom Betonwerk die 
gleiche Streuung o% aufweisen, da nicht alle Betonwerke W, iiber gleiche Mischanlagen 
verfiigen. Man nimmt daher, um zu einer groben Einteilung der Betonwerke in dieser 
Hinsicht zu gelangen, zunichst einmal probeweise an, daB die y;;, innerhalb der 

13* 
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Mischungen aller W, mit der gleichen Streuung o? normalverteilt sind, und berechnet 
fiir jedes W, mit Hilfe der Stichproben einen Schatawert 62 von 03: 


b; c 
pe Cy aaa Yas.) 
a isa (=a! 
iy b (e—1) 


Je zwei dieser. Schiitzwerte werden hierauf einem F-Test unterworfen. Das Er- 
gebnis zeigt Tab. 5: Die ersten drei Zeilen bzw. Spalten geben der Reihe nach das 
Betonwerk, den mit 10% multiplizierten Streuungsschaétzwert und die Zahl seiner 
Freiheitsgrade (F. g.) an. Die Felder der Tabelle enthalten die Quotienten der dariiber 
und der linksstehenden Streuungsschitzwerte, also die zur Priifung erforderlichen 
F-Werte. Signifikante F-Werte (« = 0:05) sind mit einem, hochst signifikante (« = 0°01) 
mit zwei Sternen versehen. 


Tabelle 5. F-Werte zu den Schatzwerten fiir die Streuung der y;;, innerhalb der 
Mischungen 


nach abnehmendem s? geordnet 


Betonwerk 
We Wi We Ws Ws Wa WwW; 
10°, s? 1769 1469 _ 1295 972 738 676 481 
F. g. 40 62 34 58 46 56 44 
Mos ae 481 44 | 3°68** | 3:05** | 2°69** | 2-02** 1:53 1:41 
fo W, 676 66 | 2°62** | 2 17** 7] 1-go% 1:44 1:09 
gs W, 738 46 | 2-40** | 1-99** | 1°75* 1°32 
a0 W,; 972 58 | 1:82* 1°51 1330 
Sa Wee) 1296 SA (137 1:13 
rs a 1469 62 | 1:20 
oo Ws 1769 40 


Kine weitere, leicht tiberprifbare Folge des in II A aufgestellten Modells betrifft 
die y;;,, deren Wv fiir Mischungen desselben Betonwerks auf Grund von I E by normal 
sein und — bei Zutreffen der Voraussetzungen — die von Betonwerk zu Betonwerk 
gleiche Streuung of + 07/3 besitzen muB. Die Vertraglichkeit der entsprechenden 
Streuungsschatzwerte wird in Tab. 6 gepriift, die fiir jedes W,; den mit Hilfe der y;; 
berechneten Schatzwert by 
a (Vis. SSP 
iBy b; —1 


von of + 02/3 enthalt und im tibrigen wie Tab. 5 angelegt ist. 


Tabelle 6. F-Werte zu’den Schaétzwerten fiir die Streuung der y;; zwischen den 


Mischungen 
eae nach abnehmendem s2 geordnet 
We W, Ws; W. Wi Wa Ws 
10°. s? 157090 138500 | 60600 55801 45632 42832 38153 
F.g, 16 21 22 19 30 27 28 
g 8 W; 38153 28 4°12** 3°63** 1°59 1°46 1:20 1°12 
As W, 42832 27 3°67** 3230" 1°41 130 1:07 
2 8 W, 45 632 30 3°44** 3°04** 1°33 1°22 
gw W, 55801 19 2°82* 2°48* 1:09 
=% W, | 60600 | 22 | 2:59* 2:29 
‘Sq W,. | 188500;| 21 | 1:13 
aS W, | 157090] 16 
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Tab. 7 faBt die Signifikanzaussagen der Tab. 5 und 6 zusammen und gestattet 
eine Gruppenbildung unter den Betonwerken, so da die Streuungsschitzwerte von 
Betonwerken derselben Gruppe miteinander vertraglich sind, das hei8t keine signi- 
fikanten F-Werte liefern. Zwar ist diese Einteilung in Gruppen nach Tab. 7 nicht 
in eindeutiger Weise bestimmt, doch wurde aus einer teilweisen Kenntnis der Misch- 
aggregate heraus die folgende Gruppierung gewihlt: 

Gruppe I: W,, W,; 
Gruppe II: W,, W,, Ws; 
Gruppe III: W,; 

Gruppe IV: W,. 


Tabelle 7. Zusammenfassung der Ergebnisse von Tabelle 5 und 6 


Nach Betonwerk 
| Tabelle | 7 |. fo ; Nl 
Ws Ws, Ws We | W, 


| z zi. | tS * OK | oe 5 ab. "Xx 
| W, 6 aw | yo ee BP pe *K eK 
| 
| 5 
W, 6 
5 
W, f 6 
5 5 
q W 
e ‘ § 
a = 
5 
5 
W, 6 
5 
W, 6 


** Hochst signifikant. 
* Signifikant. 
— Vertraglich. 


Die Gruppenbildung ist in Tab. 7 durch stark umrandete Quadrate angedeutet. 
Die bisher verwendeten Bezeichnungen werden auch fiir die Gruppen beibehalten, 
jedoch der Index der Gruppe (/ = I, II, III, IV) angefiigt. Innerhalb einer Gruppe 
sind daher a, Betonwerke W,, mit je b,;, Mischungen zusammengefaBt. Ob inner- 
halb der einzelnen Gruppen I bis IV die fiir die Zulassigkeit des Modells notwendigen 
Voraussetzungen einer statistischen Priifung standhalten, kann in zwei Schritten 
untersucht werden: 

1. Priifung, ob die Wvn der logarithmischen Druckfestigkeiten y;;, innerhalb 
der einzelnen Mischungen (Stichproben) normal sind und innerhalb jeder Gruppe 
die gleiche Streuung o?,, 1 = I,..., IV, aufweisen. Dieser Test zerfallt in die Unter- 
suchung der 

ax) Stichprobenstreuungen, 

B) Stichprobenschiefen. 

2. Nachweis, daB die Giiltigkeit der anderen Voraussetzungen, auf denen das 
Modell beruht, aus in 1 unwiderlegt gebliebenen Annahmen folgt, 


186 W. Eberl und G. SchneeweiB: 


Durchfiihrung: 1 «) Wenn die y-Wvn innerhalb der Mischungen der J/-ten Gruppe 
alle die gleiche Streuung o?, J =I,...,1V, haben (1. Voraussetzung), so ist, s0- 
3 


bald W, zur I-ten Gruppe gehoct, 8% = D'(yixn — Yu.)?/2 ein Schatzwert von 0}, 


k=1 


2 


mit 2 F.g. Nimmt man ferner an (2. Voraussetzung), daB alle Mischungsverteilungen 
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Abb. 5. Die Summenhiufigkeiten von x? mit 2 F. g. im Wahrscheinlichkeitsnetz 
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Abb. 6. Die Verteilung der Stichprobenstreuungen von 
Gruppe I (Betonwerke W 1 und W 2) 


nur nadherungsweise erfiillten ersten Voraussetzung 
zufalliger Weise von der Geraden abweichen, 


normal sind, so ist 2 s3/o?, wie 
77 mit 2F.g. (s. Abb. 5) ver- 
teilt. 


Die Richtigkeit der beiden 
Voraussetzungen wird nun 
nach 5, §. 293, gepriift. Man 
hat dazu fiir jede Gruppe die 
s?, in ansteigender Reihenfolge 
als Abszissen und die zu den 
aufsummierten relativen Hau- 
figkeiten der s}, gehérigen Per- 
zentile von y? mit 2F.g. als 
Ordinaten aufzutragen. Die so 
entstehende Punktreihe wird 
sich bei Zutreffen der Annah- 
men entlang einer Geraden 
hinziehen, die mit der Steigung 
2/o?, durch den Ursprung 
geht. Infolge der Stichproben- 
schwankungen der s?, und der 
werden die Punkte jedoch in 
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Die Abb. 6 bis 8 zeigen die den Gruppen I bis III entsprechenden Punktreihen, 
die mit der Annahme von Gau8-Verteilungen gleicher Streuung vertraglich sind. 
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Abb. 7. Die Verteilung der Stichprobenstreuungen von Gruppe II (Betonwerke W 3, W 4 und W5) 


Der Verlauf der Punktreihe von Gruppe IV in Abb. 9 laBt dagegen die Voraussetzungen 
als nicht zutreffend erscheinen. Von einer Sonderuntersuchung dieser Gruppe wurde 
jedoch Abstand genommen. 

16) In IEKa wurde festgestellt, dak 
das dritte Moment von symmetrischen yy+ 


Wvn verschwindet, wihrend es bei schiefen 
Wvn von Null verschieden ist. An Stelle 9 
dieses dritten Momentes werden fiir die 
Stichproben die ZufallsgroBen 8 
c 
Ky=c > (Yisr = Fi. FNC Ay (6 = 2) 7 
k=1 “i 
Q 
8.10 + N 6 
& 
7 N 
x f 
5 Ny 
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S ty 4 
‘ 5 
EB 4 3 
SS 
7 
7 


3 4 Io g 
Prityrobe y? Prifygribe yy? 


Abb. 8. Die Verteilung der Stichprobenstreu- Abb. 9. Die Verteilung der Stichprobenstreu- 
ungen von Gruppe III (Betonwerk W 6) ungen von Gruppe IV (Betonwerk W 7) 


ermittelt. Wenn die Mischungsverteilungen normal sind und innerhalb jeder Gruppe 
die gemeinsame Streuung o2, aufweisen, f= 1... AV. dant. ist nach * der 
Mittelwert IN {K,;} = 0 und die Streuung 


6c 
S{Ka} = (¢— 1) 2) 


8 G. A. Bennett-N. L. Franklin: Statistical Analysis in Chemistry and the Chemical 
Industry. New York: Wiley. 1954, 


a. (7) 
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Nach dem zentralen Grenzwertsatz ist der Durchschnitt K, aller K,; der /-ten 
Gruppe annihernd nach G(0,S{K,,}/g,) verteilt, wenn g, die Gesamtzahl der 
Mischungen in der J-ten Gruppe ist. Die nach G(0, 1) verteilte GréBe 
Uy= K IVS {K;;}/g, kann daher zur Priifung auf Schiefe herangezogen werden und 
ist fiir alle 4 Gruppen in Tab. 8 angegeben. Die erhaltenen u-Werte ergeben fiir die 
Yijx keinen Einwand gegen die Voraussetzung normaler Mischungsverteilungen. 


Tabelle 8. Priifung der durchschnittlichen Stichprobenschiefen fur die einzelnen 


Betonwerke 
Betonwerk 
etn: W Ws Wa W; We W, 

K 3480 | —13980| 46837 | +9413 | —4675.| 415394 | 41136 

Urwerte |//S {K}| 12971 ‘7460; 2400 3012 8121 6699 1203 
iF — 0°65 —1°87 | 2°38 BERN fe Mn een a2 42:30 | +0°94 

Logarith-| K | —39700 |—114988] +8211 | + 20332 | —18455 | +3038| + 4395 
mierte | |/S {K} 30335 49909} 12544 9957 16880 33907 6746 
MeBwerte U —1°31 — 2°30 + 0°65 + 2°04 —1:09 + 0:09 + 0°65 


Auf Grund von 1 « und 1 f wird die Annahme, daf die Mischungsverteilungen 
der y,;, normal sind und innerhalb einer jeden der 3 ersten Gruppen die gleiche 
Streuung aufweisen, den weiteren Betrachtungen zugrunde gelegt. 


2. MiBt man in jeder Mischung die y;;, von ihrem Mittelwert y,; aus, mit anderen 
Worten: geht man von den y;;;, durch 


Vise = Yerr — bes 
zu neuen ZufallsgroBen y,;, tiber, so mtissen diese fiir jede Mischung nach G(0, o?,,,) 


verteilt sein, wenn W, zur [-ten Gruppe gehort. Daher ist y,;, unabhingig von ,; 
verteilt. 


Ein bekannter Satz von Cramér® besagt: Wenn die Summe von zwei unabhingig 
voneinander verteilten Zufallsgr6Ben normalverteilt ist, dann sind es auch die 
Summanden. 


Wie Abb. 4 zeigt, kann die Verteilung der Summenwerte y,;, = 4;; + yi;, als 
Normalverteilung angesehen werden. Nach dem angefiihrten Satz kénnen infolge 
der Unabhingigkeit der Summanden y;; und y,;, diese ebenfalls als normalverteilt 
angenommen werden. Daf die Verteilung der Werte y,;, als normal angenommen 
werden darf, wurde bereits in den Abschnitten 1 « und 1 # gezeigt. Neu ist das 
Ergebnis, daB auch die w,; normalverteilt sind. 


In ahnlicher Weise kénnte auch die Unabhingigkeit und Normalverteilung der 
Bij = Mig — Mm; Und a, =p; —f innerhalb einer jeden Gruppe erschlossen werden, 
doch ist dieser Nachweis wegen der geringen Anzahl von Betonwerken je Gruppe 
(9, = 2, 3, 1,1) kaum notwendig. 


So wie Abb. 1 die x-Wvn, stellt Abb. 10 die y-Wvn dar. Wahrend Abb. 1 in 
der oberen Reihe lauter schwach linksschiefe eingipflige Kurven zeigt, die wegen 
der signifikant positiven Korrelation von 2,;, und w,; um so breiter sind, je hoher sie 
i * H. Cramér:; Uber eine Eigenschaft der normalen Verteilungsfunktion, Math, Z. 41, 405—414 

936). 
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Abb. 10, ——— Wyn der y;;, in fiinf Mischungen des 7-ten Betonwerkes; ---- Wvn der j4,; im 


z-ten Betonwerk; -.-—.- Wv der ji; 


liegen, sind in Abb. 10 simtliche Glockenkurven symmetrisch und innerhalb jeder 
Reihe kongruent. Die untere Reihe von Abb. 10 enthalt natiirlich nur Wvn von 


Betonwerken derselben Gruppe. 
C. Formeln 


Fiir jede der vier Gruppen von Betonwerken 1a%t sich nun eine Streuungszerlegung 
durchfiihren, wobei man von der Identitat 


= (Yer yy =e (Yas 4a} Fee (yy sty ey ae oa Ye Y..P (8) 
ausgeht. Bezeichnet man Durchschnitte der «,, B;;, yi;, ebenso wie die der x,;;, und 
Yesr, 80 gilt: Yin = b+ OH + Bis + Vise 

Yir = Pe %e + Ba + Yu. 

Yi. = fe + oe + Bi +. 

A fs ee bh +e b; o/B, +B. +y., 


(9) 
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mit B, Sb. Der Ausdruck Pane: «,/B, ist der mit den b; gewogene ee eect der 
«, und nur im Falle b, = b, =... = 6, gleich dem Durchschnitt «, - Dx o,/a7°. 
Sind 2,,...,%, ” voneinander unabhingige Beobachtungen einer nace G(u, 0?) 


verteilten ZufallsgréBe z, so ist bekanntlich ° 


wt | Se, — 3H = (nm — 1) 
¢=1 


nv 
wenn z = .’z,/n der Durchschnitt der Beobachtungen ist. Setzt man nun die rechten 


i=1 F 
Seiten von (9) in die Quadratsummen auf der rechten Seite von (8) ein, so gewinnt 
man daher als Mittelwerte dieser Quadratsummen der Reihe nach: 


I. mI = (Yaar — ys) | = Po io (Vis — vas) = B(¢— 1), 


YY 


2. MCS” (vs, — yen) = 08’ MLS (Bis — Be) + (ren — 74.) Ph = 


= 5 [RIE (Bos — Behl + MLS Cs. — 7.3, 
da der Mittelwert des doppelten Produktes wegen der Unabhangigkeit der 6 und y 
und wegen N{6;; — B:.} = M{y.;. — y:..} = 90 verschwindet, 
M Jed’ Ys. — y..)*\— oS) [(b: — 1) 69 + (bi — 1) 2 Jel, 
da y,;, als Darcie aus ¢ Beobachtungen mit der Streuung o? selbst die Streuung 


oz/c besitzt. 
Man erhalt also endgiiltig 


mM {ex (43, = ¥.) = (B, — a) (ec oF ae 02). 
3. M fe Ye. = y...\ =e Ms 2 ( = 2"; a/By) + (B:. — B.) + 
+ (ys. — yaa] | = 2.8m { (ae — Xb; 0)/ By) + M {Be = BP} + ML. 7.9]. 


a 


aj 


Die drei Mittelwerte in der eckigen Klammer sind der Reihe nach, wenn man 
a 
) 0? = B, setzt: 


ji 


M (a? — 2 0; 3b; 05/By + (2 b,j [Bil = 08 — 2b, 03/B, + By o2/By, 
i] j J 


Me { (Ba. shy eet ao 1) of/a b,, 
Me {(y A 1) o2/a b,c. 
Man hat daher schlieBlich 


Me ex i y..)| = Opn: (o% — 2b, o2/B, + By 0% Bi) + ¢ (a — 1) oF a3 
Yy a 

2 Bi— B, 

+(a—1)e8=c > — o2 + (a — 1) of + (a— 1) 02. 

Diese Ergebnisse sind in Tab. 9 festgehalten. Als mittlere Quadrate werden dabei 


die tn cae von Quadratsumme und zugehériger Anzahl der Freiheitsgrade be- 
zeichnet 


** An Stelle von a,; und 6;, wird kurz a; und b; geschrieben, 
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Tabelle 9. Das Schema der Varianzanalyse einer Gruppe mit mehreren Beton- 


werken 
A Mittlere Mittelwerte 
eee | a ee ae Quadrate der mittleren Quadrate 
Site 5(4%;—y.?= oj = oF + ¢0% + 
Zwischen den a me ee Se 
Werken Caan! = i ea el jf | (B? — B,)c 
oe B, eT ae a = 3 
2 petty pee Nae 
Zwischen carey x (Yas. — Yi.) es 
den Stich- 2) ft. = B,—a ieee 6 = o%+co? 
1 1 y2 2 1 2 2 y B 
proben eS 2 Ve ae ties 1p 
aj 
2 Yrsn—Y;;, P= 
Tnnerhalb S3 ow mo Yis.) 2 
der Stich- Mg fg = (€— 1) B,| 3 = B= 02 
proben = =» 54 — 2aDn ya, fs us 
ijk Toe 
Summe S2= > (Yisn— Y...)? f= Bye 1 
ijk 


Diese mittleren Quadrate sind Schatzwerte ihrer Mittelwerte o?, i = 1, 2,3. 
Durch sie gelangt man zu Schatzwerten der Varianzkomponenten o?, 63, o?. Das 
Zeichen ~ bedeutet im folgenden immer ,,wird geschatzt durch“ bzw. ,,ist ein Schatz- 
wert von‘, je nachdem der unbekannte Parameter, z. B. die Varianzkomponente, 
auf der linken oder rechten Seite der Gleichung steht. 


Ks gilt also: 


. B?— B,) 

2 (Bi 2 2 2 oe 

oo Gis oO, + COR + OF & Si, 

a= cof + o2 w 88, (10) 
oo = of & 83. 


2 = 2 2 
~7 8 85) =S 
Ce (Bi Be) e Sige) 
82 — 82 
2 2 3 ery 
OB C a SB; 
y WY OZ y* 


Die Streuung der fiir die Varianzkomponenten ermittelten Schatzwerte 1aBt sich 
auf dem Weg iiber die Streuung der s?, 7 = 1, 2, 3, bestimmen. Bildet man mit Hilfe 
von Tab. 9 die Si/0%, so sind diese GroBen wie Chiquadtat mit f; F. g. verteilt. Die 
Streuung eines 7? mit fF. g. ist bekanntlich > 2/, daher gilt fiir die Streuung von 
S?2/o?: GS {S?/o?} = 2f, und fir die Streuung von s} = S?/f;: S {s}} = 2 of/f,. 
Die s? sind aus voneinander unabhingigen Quadratsummen hervorgegangen und 
daher selbst voneinander unabhangig. Leider hingen ihre Streuungen von den un- 
bekannten o? ab, fiir die wir nur Schatzwerte haben. Erwartungstreue Schatzwerte 
der 264/f; bekommt man so: Es ist I {s?} = of und © {s?} = ne Nach dem 
Verschiebungssatz ist M {sf} = S {s?} + [M {s}7 = = 2 of/f, + oF OEY) ott 
Daher ist It {2 s#/(f; + 2)} = 2o7/f; und somit 2s; tl(f; + 2) ein Ser aca 
Schaitzwert von © {s?}. 
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Infolge der Unabhangigkeit der s? erhalt man daher naherungsweise als Streuungen 
der Schatzwerte fiir die Varianzkomponenten: 


9 2 

Sie 
2 

S {8} == 


S {sy} = 


a — 1)? B3 


83 


( 


1— B,) 


f+ 2 


4 
2 83 


fy + 2° 


fe + 


83 
9 bs 


st 
nee) 


Im Falle o2 = 0 ist s?/s} wie F mit f, und f, F.g., im Falle of = 0 ist s3/s3 wie F 
mit f, und f, F. g. verteilt. Ein zu groBer Wert eines der beiden Quotienten zeigt an, 


daB o2 bzw. o} signifikant von 0 verschieden ist. 


D. Durchftthrung 


Die Tab. 10 bis 13 enthalten nach dem Muster von Tab. 9 die Varianzanalysen 
fiir die Betonwerksgruppen I bis IV. Da die Gruppen IIT und IV jeweils nur aus einem 
Werk bestehen, entfallt in den Tab. 12 und 13 die Zeile ,,zwischen den Werken“. 


Tabelle 10. Die Varianzanalyse der Gruppe I 


Zwischen den Werken..... 


Zwischen den Proben 


Innerhalb der Proben 


She, sek 


eee | 1-612775 


S2 


0°722150 
0°728749 
0°161876 


0°722150 
0°0148724 
0:00158702 o. 


o2 


+3 of + 72°9412 of 
oy 3 oF 


Y 


Tabelle 11. Die Varianzanalyse der Gruppe II 


Zwischen den Werken..... 


Zwischen den Proben 


Innerhalb der Proben 


ey: | 2-261643 


S? 


1:066093 
1067386 
0°128164 


0°533047 
0°01386216 
0°000801025 o 


o2 


a + 3 of + 79°6125 of, 
o, eS OB 


Y 


Tabelle 12. Die Varianzanalyse der: Gruppe III 


Zwischen den Proben 
Innerhalb der Proben 


S2 


0°760708 
0:044039 


0°804747 


F, g. 


50 


3? o 
0°0475443 oy + 3 of 
0°0012953 oy 


Tabelle 13. Die Varianzanalyse der Gruppe IV 


Zwischen den Proben 
Innerhalb der Proben 


32 o2 


a | eee ee eee 


@) Ss ole © ete 


9X0) (0). 916) 6. (w 


S? F. g. 
0°872525 21 
0°021161 44 
0°893686 65 


0:0415488 oy +3 oR 
0:00048093 o 
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Tabelle 14. Ergebnisse der Varianzanalysen 
SSS SEES EEE a ee ee eee ee 
Schitzwerte der Varianzkomponenten 


i s2- 10° s2- 108 8,, ° 10 Lie pd 
a B oy 
I 9697 + 8089 4428 + 985 1587 + 220 0°015712 
IT 6521 + 4735 4354+ 736 801 + 89 0°011676 
TIT — | 15416 + 5284 1295 + 305 0°016711 
BV; — 13689 + 4089 481 + 100 0°014170 


Tab. 14 enthalt die Schatzwerte der Varianzkomponenten mit ihren Standard- 
abweichungen (Wurzeln aus den Streuungen.) 


E. Ergebnisse 


Fir die Untersuchung der Streuungen wurden nicht die Druckfestigkeiten selbst, 
sondern deren Logarithmen herangezogen, weil diese eine Streuung zeigen, die nicht 
mehr von der Hohe der Druckfestigkeit abhingig ist (vgl. I E c). 


In Tab. 14 sind die Schitzwerte der Varianzkomponenten und als deren Summe 
der Schaitzwert der Streuung der y;;, innerhalb der Gruppen angegeben. 


Die Streuung o2 ist ein Ma fiir die Unterschiede der Betonwerksmittelwerte. 
Diese Streuung besitzt im vorliegenden Fall insofern eine geringere Bedeutung, als 
die Unterschiede der Mittelwerte der Betonwerke nicht als rein zufallsbedingt anzu- 
sprechen sind. Die GréBe von o2 legt die Vermutung nahe, daB jedes Betonwerk 
seiner Betonherstellung einen seinem Sicherheitsbediirfnis entsprechenden und gegen- 
tiber der Norm erhohten ,,Sollwert®‘ zugrunde legt. 


Wie bereits eingangs erwahnt, gibt die Streuung of ein MaB fiir die GleichmaBig- 
keit der Zusammenmischung, wahrend o? ein Maf fiir die Gleichmafigkeit der Durch- 
mischung des Betons ist. Ein Vergleich der Streuungen of und o, zeigt, daB die 
groBen Streuungen der Druckfestigkeiten im wesentlichen auf einer mangelhaften 
Kontroile der Betonzusammenmischung beruhen. 


Eine gute Durchmischung bedeutet eben noch keine gute Zusammenmischung. 
Die Streuungen der Gruppe IV zeigen dies deutlich. 


Die in Tab. 14 angefiihrten Schaitzwerte der Varianzkomponenten k6énnen zur 
Uberpriifung der Arbeitsweise der bestehenden Betonwerke herangezogen werden. 


Dies kann iiberschlaglich in der Weise erfolgen, daB man, ausgehend von der Beziehung 


mM {ox (Gi y..) = (b; — 1) (c of + of) 


den Schitzwert 
& (Yi4. — YixY 
J 
eee OF 


unter Vernachlissigung von o2 errechnet und mit den Werten sj in Tab. 14 ver- 
gleicht. Sicherheitshalber kann nachgepriift werden, ob 
PS (Yisn — Ye3-)° 
Ae See ye 
b, (ec — 1) es 


tatsichlich gegentiber of vernachlassigt werden kann. 


Das der Untersuchung zugrunde gelegte Varianzkomponentenmodell bringt die 
zusammengesetzte Natur der Gesamtstreuung klar zum Ausdruck und wird die 
Grundlage fiir ein stochastisch fundiertes Priifverfahren liefern, das auf die Uber- 


194 G. Schumpich : 


wachung sowohl des Zusammen- als auch des Durchmischens gerichtet sein wird. 
Die bestehenden Betonnormen tragen diesem Gesichtspunkt nicht Rechnung und 
sind in dieser Richtung zu erweitern”. (Bingegangen am 4. April 1957) 


11 Die Verfasser sind gerne bereit, Betonwerken, in denen soleche Kontrollen von Hilfskraften 
durchgefiihrt werden, entsprechende Vordrucke zu senden. : 

Die folgenden beiden Biicher beschreiben Methoden, die in dieser Arbeit Verwendung finden, 
wiahrend sich die beiden Abhandlungen auf denselben Gegenstand beziehen: 

O. Kempthorne: The Design and Analysis of Experiments. New York: Wiley. 1952. 

A. Linder: Planen und Auswerten von Versuchen. Basel-Stuttgart: Birkhauser. 1953. 

H. Risch: Betrachtungen zur Priifung der Betonfestigkeit. Beton- u. Stahlbetonbau 51, 
135-138 (1956). 

W. Arneth: Auswertung der Giitepriifungen an Betonprobekérpern und Bohrkernen. 
StraBe u. Autobahn 7, 424—427 (1956). 


Beitrag zur Kinetik und Statik 
ebener Stabwerke mit gekriimmten Staben** 


Von G. Schumpich, Hannover 
Mit 12 Textabbildungen 


Zusammenfassung. Ausgehend von Differentialgleichungssystemen fiir freie Schwingungen 
eben gekriimmter Stabe wird eine Matrizenmethode, das sog. Reduktionsverfahren, zur Berechnung 
der freien Schwingungen ebener Stabwerke mit gekriimmten Staben entwickelt. Das Verfahren 
wird auf erzwungene Schwingungen und die Statik solcher Stabwerke ausgedehnt. 


I. Einleitung 


Die auf dem Gebiet der Dynamik gekriimmter Stabe bis jetzt erzielten Fortschritte 
findet man in dem Buch von K. Federhofer ,,Dynamik des Bogentragers und 
Kreisringes‘‘ [1]? zusammenfassend dargestellt. Da die exakte Losung schon ver- 
haltnismaBig einfacher Probleme nur mit einem sehr groBen Rechenaufwand erhalten 
werden kann, ist man gezwungen, Naherungsverfahren anzusetzen. So berechnen 
verschiedene Verf. Eigenfrequenzen und Eigenschwingungsformen mit Hilfe des 
Ritzschen Verfahrens. 

In der letzten Zeit ist eine Reihe von Arbeiten erschienen, unter anderen die 
Arbeiten [5] bis [11], in denen verschiedene technische Probleme nach einer Matrizen- 
methode gelost werden, welche jetzt als Reduktionsverfahren bezeichnet wird [6]. 
Grundsatzlich laBt sich das Reduktionsverfahren auf alle Probleme anwenden, die 
durch gewohnliche lineare Differentialgleichungen beschrieben werden kénnen. In 
der vorliegenden Arbeit wird das Reduktionsverfahren fiir beliebige ebene Stabwerke 
mit gekriimmten Staben entwickelt. 

Zur Herleitung der sog. Ubertragungsmatrizen des Reduktionsverfahrens wird 
zuerst eine allgemeine Methode vorgeschlagen, bei der man von einem System 
von  linearen homogenen Differentialgleichungen 1. Ordnung mit konstanten 
Koeffizienten ausgeht. Dabei hat man die Vorteile, da man sich von vornherein 
des Matrizenkalkiils bedienen kann, und da® daritber hinaus die etwa notwendigen 
Grenziiberginge bei der Herleitung der Ubertragungsmatrizen fiir Sonderfille 


praktisch schon an der Koeffizientenmatrix des Differentialgleichungssystems durch- 
gefiihrt werden kénnen. 


* In Verehrung Herrn Professor Dr. K. Federhofer gewidmet 
1 Gektrzte Fassung der Dissertation des Verf., Hannover 1957; Referent: Prof. Dr.-Ing. 
E. Pestel, Korreferent: Prof. Dr.-Ing. K. Jaeckel. 


® Die Zahlen in eckigen Klammern verweisen auf das Schrifttumsverzeichnis am Ende der 
Arbeit. 


Beitrag zur Kinetik und Statik ebener Stabwerke mit gekriimmten Stében 195 


II. Mathematische Grundlagen fiir die Herleitung der Ubertragungsmatrizen 


Wir betrachten das System von n linearen homogenen Differentialgleichungen 
1. Ordnung mit konstanten Koeffizienten 


dz, 
Tage = C1 aa r Ayn &n 
S (neice) hee ersel allele cel hue seria, \e. p—sle © 9 sie’ 1 a 
din oe 
Fe On 1 + Ann Sn, 
welches in Matrizenschreibweise wie folgt geschrieben werden kann: 
dz, 
ee (1b) 
3 d F 
mit den Spaltenvektoren 3(s) und ae und der quadratischen n-reihigen Koeffizienten- . 
matrix 1. 
Die Lésung des Systems (1) bei den Anfangsbedingungen z, (so) = 2,°,.... 2n (So) = 2m 
wobei wir s, = 0 setzen wollen, ist nach [3], S. 120, gegeben durch 


3(s) = Us) - 30 (2) 
mit dem Anfangsvektor 3, = 3(0) und 


U(s) = eds, (3): 


Dabei soll Ul eine quadratische n-reihige Matrix bedeuten. Mit Hilfe der Matrix 11 
kann aus dem Vektor 3, an der Stelle s, = 0 nach (2) der Vektor 3(s) an einer be- 
liebigen Stelle s in einfacher Weise berechnet werden. Die Matrix ll nennen wir 
, Ubertragungsmatrix“’ und den Vektor 3(s) ,,Zustandsvektor, da er den Zustand 
an der Stelle s vollstandig kennzeichnet. 

Die durch (3) eingefiihrte Ubertragungsmatrix 11 kann aus % mit Hilfe der 
Definitionsgleichung 


le E+ As +s: (Use +a, Ashe+..., (4) 
also durch Aufsummierung einer unendlichen Reihe, berechnet werden. Jedoch 


14Bt sich e%* auch durch einen endlichen Polynomausdruck, das sog. Ersatzpolynom 
(vgl. [4], S. 227) von der Gestalt 


els —¢,E+e,Ast... +n_, (As)! = g(As) (5) 


darstellen, so daB die Matrix WU schon durch Aufsummierung von endlich vielen 
Matrizen gewonnen werden kann. 


Die Koeffizienten c, in (5) kénnen nach [4], 8. 227ff., bestimmt werden. Die - 


dort beschriebene Methode zur Berechnung der c¢, sttitzt sich auf den folgenden 


Satz: Die charakteristischen Zahlen x; eines Matrizenpolynoms oder einer durch 
ein Polynom darstellbaren Matrizenfunktion © = /(¥) sind mit den charakteristischen 
Zahlen A; der Matrix 8 verkniipft durch die Gleichung 


= f(A;). 3 (6) 
Da die gegebene Funktion /(8) und das Ersatzpolynom g(%) eine und dieselbe Matrix © 
darstellen, so haben beide dieselben charakteristischen Zahlen x,, und nach dem 


obigen Satz folgt 
g(Ai) = f(Ai). (7) 


3 In dem uns hier speziell interessierenden Fall ist 8 = %-s, f(8) = els, = > Zu,den | 


nachfolgenden Ausfiihrungen vgl. [4], S. 227—236, insbesondere auBer dem hier zitierten Satz 1 
auch den Satz 2. 
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Wir wollen hier noch den folgenden Zusatz beweisen, mit dessen Hilfe die Berechnung 
der Koeffizienten c, des Ersatzpolynoms g(%) auch im Falle mehrfacher charakteri- 
stischer Zahlen von % sich einfach gestaltet. 
Zusatz: Ist A, eine p-fache charakteristische Zahl von %, so gilt auBer 
gy) =f) | 
auch | g'(Ax) = f(r), 
g'"(Ay) = F'(Ay); 


yfo-0(2,) ca fe Aa) 


wobei Ableitungen nach dem Argument durch Striche gekennzeichnet sind. Diesen 
Zusatz beweisen wir mit Hilfe von Grenziibergingen, die den Grenzibergaingen bei 
Zurmiihl [4] entsprechen. Dabei gehen wir zweckmaBig von der Darstellung des 
Ersatzpolynoms in der Gestalt aus: 


9(B) = 6, € +b, (B — A, E) +... + bps (B — A, E) (BW — ALE)... (W— An). (9a) 
Dem Ersatzpolynom entspricht dann das Interpolationspolynom 
g(a) = by + 6, (@ — Ay) +. ~~ +. Ona (@ — Ay) (% — Ag)... (@ — Ana). (10a) 


Es mégen zunachst simtliche charakteristischen Zahlen 4; von 8 voneinander ver- 
schieden sein. Dann hat man nach (7) mit g(%) nach (9a) in 


f(A1) a bo, ) 
f(A2) oo by ale bi(Az ag Ax); 


f(As) or by aie b,(A3 | Ay) ie bo(A3 am, A) (As we Ag); 
: (1la) 


ein Gleichungssystem von 7 linearen Gleichungen vorliegen, aus dem sich die  un- 
bekannten Koeffizienten 6b, rekursiv bestimmen lassen. 
Aus der ersten Gleichung von (lla) folgt 


bo = f(A,). 
Die zweite Gleichung bringen wir auf die Gestalt 


As) — f(A 
b, = Ae (12a) 


Macht man nun den Grenziibergang 1, — /,, so wird der Koeffizient b, = f(A,) 
durch diesen Grenziibergang nicht beeinfluBt, b, geht nach (12a) iiber in 


b,! = lim b, = f(A) (13a) 


Ag > Ay 
und die tibrigen Koeffizienten b, in b,1. Man erhalt nach dem Grenziibergang statt (9a) 
g(B) = f(Ay)- € + f(A) (BW — A, E) + by (B— A, EP +... 4 
Tt Oya (B — 2, EP» (B— AE)... (B— 1,1: ©, (9b) 


statt (10a) 
g(x) = f(Ay) + P(Ay) (@ — Ay) + by (w — aR... + 
+ hg ag (@ — A) (ge een) (10b) 
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und statt (lla) 


f(A) = bo + by (As — A) + 851 (Ag — Ay), 
(11b) 


f(An) = bo + 64" (An — Ay) + 83! (An — Ay) cet 


oh Lanner (An Ty Ay)? On ae As) . ae at See 
Aus (10b) folgt nun ‘ 
g(Ay) = f(y), 
Ay) = f'(Ax) 
Die Behauptung ist damit fiir eine zweifache charakteristische Zahl richtig. Um 
den Satz allgemein fiir eine p-fache charakteristische Zahl zu beweisen, wenden wir 


den Induktionsschlu8 an. Der Satz sei richtig fiir eine (p — 1)-fache charakteristische 


Zahl 2,. Das Gleichungssystem fiir die Koeffizienten des Ersatzpolynoms lautet in 
diesem Fall 


f(A1) = 5o, | 
f'(ay) = by, 
f(A) = 2! b.”, 


f(A ) sa bo, 
f(A) = by, 


: (11c) 
f(A.) = (pp 2)1 088, 


Ag) 06 Or ap a ke enn (Ap — A)? On (Ay AP 


Og Dee Aa se ag Ag A)? os (An Aga) 
LaBt man nun auch A, — 4, gehen, so erhalt man aus der p-ten Gleichung von (11c) 


= OA ali i a 


Ap >A 
(p—2)(4 st 
gg MANN f 1A) Op AY 2 Fst ar? 
oe ey, (Ap ss Aye e 
(p—1) 
is pane ae 
Die Koeffizienten by, b;,..., 02-3 werden durch diesen Grenzproze8 nicht be- 


einfluBt, die ibrigen Rooffizienten 1 p-2 gehen in b,?-1 iiber. Nach dem Grenz- 
iibergang erhilt man ein dem System (11 ¢) entsprechendes System, welches aus diesem 
durch Verschiebung des Index p um eins entsteht, und somit das Interpolations- 
polynom 

4 Das Bildungsgesetz des Gleichungssystems (llc) erkennt man leicht, wenn man noch 
in (11b) den Grenziibergang A; —> A, durchfihrt. 
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fO-Y(Ay) p-1 ‘ 
gst) = faa) + fA) (@ — Aa) +. + it ( — Anas) + 8B (@ — AP + 
+ OPS (@ — Ay)? (3 — Aga) +... + BRIT (a — Ay)? (@ — Apia). ++ (%— dna). (10d) 
Aus (10d) folgt die Behauptung 
(Ax) = f(Ar), 
g' (Ay) = f'(Ad); 


gay) = fO-M(A,). 
Wir bemerken noch, da8 die Koeffizientenmatrix des Gleichungssystems (11) 
eine untere Dreiecksmatrix mit von Null verschiedenen Elementen in der Diagonalen 
ist, so daB die Koeffizientendeterminante nicht verschwindet und das Gleichungs- 


system fiir die b, eindeutig losbar ist. 
Die Koeffizienten c, des Ersatzpolynoms in der Gestalt 


g(B) =¢,E6+¢,8+...+¢,_, 8! 
bestimmen sich nach dem vorangehenden Satz durch Lésung der Interpolations- 
aufgabe, dasjenige Polynom g(x) vom niedrigsten Grade zu bestimmen, welches 
an den Stiitzstellen x; = A; die vorgegebenen Funktionswerte g(/;) = f(A;) = f; und 
Ableitungen 9(A,;) = f(A,;) =f (7 = 1, 2,3,...,p;— 1) annimmt, wobei p; die 
Vielfachheit der charakteristischen Zahl A, ist. ; 


Ill. Differentialgleichungen fiir die Zentrallinie gekriimmter Stabe. 
Bewegungsgleichungen fiir freie Schwingungen 
Die Zentrallinie eines im Sinne der mathematischen Elastizitatstheorie diinnen, 
raumlich gekriimmten Stabes sei durch 


t=r(s), s — Bogenlange, 


gegeben. Jedem Punkt der Raumkurve r ordnen wir ein rechtwinkliges 2, y, z- 
Koordinatensystem zu, dessen x-Achse in Richtung des Tangentenvektors t der 
Raumkurve verlauft und dessen y- und z-Achsen 

te in der Normalebene liegen und mit den Haupt- 

; Me tragheitsachsen des Stabquerschnittes zusam- 


ty 


fk 


My 
“Se 


Abb. 1 Abb. 2 


menfallen. Die Kinheitsvektoren i,j, in Richtung der Achsen bilden ein Rechts- 
system. Den Winkel zwischen der y-Achse und der Hauptnormalen n bezeichnen 
wir mit y. Fir y¥=0 geht das Hauptachsensystem i, j, £ in das Hauptsystem 
(begleitendes Dreibein) t, n, 6 iiber. ; 
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Die Verschiebungs- und SchnittgréBen definieren wir dann als positiv, wenn ihre 


Vektoren am positiven Schnittufer in die Richtung positiver Koordinatenachsen 
weisen (Abb. 1 und 2). 


Aus den bekannten Gleichgewichtsbedingungen und den Elastizitatsgesetzen fiir 
Biegung, Dehnung und Torsion (vgl. [2] und [1]) erhalt man unter Benutzung des 
oben eingefiihrten x, y, z-Koordinatensystems das folgende Differentialgleichungs- 
system fiir die verformte Zentrallinie®: 


: 1 
Gu ER Ve k,v ky w, 


Qe = k,Q, — ky Q, — Pay 
0 =cM,—k,p +k, y, 
M, =—k,M,+k,M, — dy; 
v=—ku+tot+rty, 
yp =k,d+aM,—rty, { (13) 
M, =k,M,—Q,—tM, =O, 
Q,y = —k,Q.+7Q, — Py» 


w =k,u—tv—y, 


yp =—k,6+179+bNM,, 


Meh eM, Oa 
emer s aes d 
Q. = ko, — tQ, — Pp: ( ==). J 
Dabei gelten folgende Bezeichnungen: 
U,V, W Verschiebungen in der Richtung der z-, y- und z-Achse. 
6, Y, Y Winkel der Drehung des Stabelementes um die z-, y- und 
z-Achse. 
i Torsionsmoment. 
M,, M, Biegemomente. 
OQ. Normalkraft. 
Oe.0, Querkrafte. 
F Querschnittsfliche des Stabes. 
=GIp Torsionssteifigkeit. 
+ = HT, = = EI, Biegesteifigkeiten. 
= Windung der Zentrallinie. 
d 1 
Baa te 
= k Kriimmung der Zentrallinie. 
| = cos ¥, ky = = sin yx Kriimmungskomponenten. 
Dar Py> Dz Komponenten der duBeren Krafte (bzw. Massenkrafte) je 
Langeneinheit der Zentrallinie. - 
Gy dysndy Komponenten der auBeren Momente (bzw. Massenmomente) je 


Langeneinheit der Zentrallinie. 


5 Im Differentialgleichungssystem (13) ist der Hinflu8 der Querschnittsverwélbung infolge 
der Torsion nicht beriicksichtigt; Naheres dariiber s. in [1]. 
14* 
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Die GréBen k,, k, und 7 beziehen sich auf den unverformten Zustand der Zentral- 
linie. 

Wirken auf den Stab wihrend der Bewegung keine auBeren Krafte, so wird seine 
Bewegung durch das Differentialgleichungssystem beschrieben, welches aus (13) folgt, 
wenn man in diesem fiir p,, Py, Dz» dy, dy, d, die Tragheitswiderstande einsetzt, also 


Po= — "4, iS Py = = ph, Pei kD; 

dy, = — Ut," 0, dy = — wt, 9, d, = — wie 9. (14) 
bs Masse des Stabes je Langeneinheit der Zentrallinie. 
tay ty, tz Tragheitsradien beziiglich der 2-, y-, z-Achsen. 
t Zeit. Ableitungen nach der Zeit werden durch Punkte gekennzeichnet. 


Mit dem Ansatz fiir Higenschwingungen 
u(s, ¢) = u(s)e>t, 
Q2(S; t) = Qa(s) e%, 
6(s, t) = 6(s) e?*, 


rigs Tera bs (15) 


ee 


Q.(8, t) = Q,(s) eet J 
geht (13) mit p,, p,,... nach (14) in ein System von zwolf gewohnlichen linearen 
Differentialgleichungen 1. Ordnung fiir die zwolf unbekannten Funktionen w(s), 


Q.(s), 6(s),...,Q,(s) tiber*, welches bei Verwendung der Matrizenschreibweise in 
der Form fhe g's (16) 
ds 
geschrieben werden kann mit dem Zustandsvektor’ 
3= {u,Q2, 6, Mz, v, p, M,, —Qy — wv, y, Uy Q} (17) 
0 a 0 0 k, 0 0 0 ky 0 OF p28 ; 
Swot o0e| 0 Oo Oso Gf ES 0 aot eeienees 
ee ee on eee oe ee meee 
0 | 0 Lpeatiz| 0 Ee Wo soe che eee 
ae Oe eo 0 “Qe 0 le oe Gal eee 
es ee a0 10 eae as) 
VS ORR 0 It ae 6.) tera 
Oke e 074.00 0.1 0-1 O. Ol ones 
hy | 0.) 0. |) 0.) yO) (Oy ental Om eae 
a ee kegs 
Per ey Pn zs | 0 [0 Epo 0} 1 
Oke oe ee RG aed Fekete OAC FT & 


8 Fir die Funktionen u(s, t), Q,(8, t), d(s,t),... und die reinen Ortsfunktionen u(s), Q (8), 
6(s),... verwenden wir hier dieselben Funktionssymbole, was jedoch nicht zu Verwechslungen 
fiihren dirfte. 

” Spaltenvektoren werden aus Griinden der Raumersparnis in Zeile geschrieben und durch 
geschweifte Klammern gekennzeichnet. 
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und der quadratischen Matrix % (18). Die stark umrahmten Teilmatrizen in der 
Hauptdiagonalen der Matrix % bedeuten der Reihe nach von links oben nach rechts 
unten die Anteile, die zu Normalschwingungen, Torsionsschwingungen, Biege- 
schwingungen in der y-Richtung und Biegeschwingungen in der z-Richtung gehéren. 
Die Matrizenelemente auBerhalb der eingerahmten Teilmatrizen haben die Be- 
deutung von Kopplungsfaktoren der verschiedenen Schwingungsarten. 


Aus dem allgemeinen System (16) mit der Koeffizientenmatrix (18) lassen sich 
leicht Differentialgleichungssysteme fiir verschiedene Sonderfalle gewinnen. Fiir den 
Sonderfall der freien Schwingungen eines eben gekriimmten unverwundenen Stabes, 


dessen Zentrallinie in der «,y-Ebene liegt (y = 0, = = Oct 0, ho i), 
zerfallt das allgemeine Differentialgleichungssystem in zwei voneinander unabhingige 
Systeme: 


A. Schwingungen in der Ebene der Zentrallinie. Kopplung der Normalschwin- 
gungen mit den Biegeschwingungen in der y-Richtung. 


u U 
Qe Qn 
ove 0 0 (19) 
y 0 0 , 
MU,’ Oe O. APO etFANO saul M, 
=O 0 Py Wie WA!) 0 0 =), 


B. Schwingungen senkrecht zur Ebene der Zentrallinie. Kopplung der Torsions- 
schwingungen mit den Biegeschwingungen in der z-Richtung. 


é | 0 Gus LaeOeh ahohane é 

Me —porie 0 0 0 0 M, 

4D Wee 0 0 0 1 0 O21 (=o (20) 
y ay 0 i) 0 b 0 y 

Bd 0 =k] 0=pw*1,2 0 1 M, 

gi 0 0 juw 0 0 0 Q, 


Diese beiden Schwingungsarten eines eben gekriimmten Stabes sollen nun im folgen- 
den Abschnitt behandelt werden. 


IV. Freie Schwingungen eben gekriimmter Stabe 
A. Herleitung von Ubertragungsmatrizen 


Fall A: Schwingungen in der Ebene der Zentrallinie. 


1 
Fiir den Fall eines homogenen Kreisbogens (k = konst., 1 = konst., HL, = 7 = 


= konst.) sind die Elemente der Koeffizientenmatrix des Differentialgleichungs- 
systems (19) konstant. Zur Aufstellung der Ubertragungsmatrix fir ein solches 
homogenes Kreisbogenstiick mit dem Zentriwinkel ® und dem Radius @ (Abb. 3) 
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kénnen wir demnach die im Abschnitt II beschriebene Methode benutzen. Setzt 
man fiir die Bogenlinge des Kreisbogens 


a= = == (P.O (21) 
und fihrt die Abkiirzungen ein 
EJ, 
@ or 0" — ae gpooh, (22) 


so erhalt man fiir die Ubertragungsmatrix zunachst die Dar- 
stellung 


U = exp Us = exp -@ (23) 


Wir bezeichnen die charakteristischen Zahlen der Matrix Us = %X%o® mit A; 
und die der Matrix 2% oe mit A;. Es gilt 


Die charakteristische Gleichung der Matrix 0 lautet: 


det (EA—%e) = 4° + A244 +44] 4 


2 it itera All cae op Sea ae 
+@{1-A +24 ALLA aa - 41+ -—AS_—A 7 |=0. (25) 
Bei Vernachlassigung des EKinflusses der Rotationstragheit der Stabelemente 

und der Voraussetzung dehnungsloser Deformation erhalt man aus (25) mit 7, = 0, 
f =O 
A’ +2A4+ A27(1—A)+A=O. (26) 


Schon die Diskussion dieser einfacheren Gleichung bereitet Schwierigkeiten. 
Die Koeffizienten und damit auch die Wurzeln der in A? kubischen Gl. (26) sind 
frequenzabhangig. Man muB daher verschiedene Frequenzbereiche unterscheiden 
und fiir jeden dieser Bereiche gesondert eine Ubertragungsmatrix aufstellen. Dariiber 
hinaus wiirde das Arbeiten mit diesen Ubertragungsmatrizen unhandlich sein, weil 
ihre Elemente einen komplizierten Aufbau haben. Aus diesen Griinden wollen wir 
auf die Aufstellung der Ubertragungsmatrizen fiir einen homogenen Kreisbogen in 
dieser Allgemeinheit verzichten® und beschrinken uns auf die Herleitung von Uber- 
tragungsmatrizen fiir folgende, bei den praktischen Anwendungen wichtigen Sonder- 
falle: 

a) Masseloser Kreisbogen mit dem Zentriwinkel ® und dem Radius 0. (u = 0). 

b) Starrer Kreisbogen mit dem Zentriwinkel ® und dem Radius o. (a = 0). 
c) Einzelmasse m. 
d) EHinzelelastizitat e. 


8 Bendtigt man eine solche Ubertragungsmatrix in einem konkreten Fall, so ist es bequemer, 
YX o® < 


sie numerisch fiir eine feste Frequenz w direkt aus der Reihe U = ei??? — G+ 1! 


2 
He sea! +... zu berechnen. 
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Diese gewinnen wir durch Grenziiberginge aus der Darstellung (23) der Uber- 
tragungsmatrix fiir einen homogenen Kreisbogen. Wegen der Stetigkeit der 
Exponentialfunktion gilt 


lim exp Us = exp lim As = exp B. (27) 
Nachstehend sind fiir die oben genannten Sonderfalle die Exponentenmatrizen 8, 


ihre charakteristischen Gleichungen und charakteristischen Zahlen zusammen- 
gestellt. 


a) 


3, = lim We = (28a) 
ud 
Boia Os A= 1, =.0; 
b) 
GS, =m Ws (28b) 
a—-0 


A2(1 + A= 0; 4, = 
a) 
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In den Sonderfallen c und d lautet das Ersatzpolynom fiir e® 
1 1 if 
g(B) = C+ 894+ 845, B+ _ Bt aw (29) 
Die Koeffizienten c, in den Sonderfillen a und b im Ersatzpolynom 


g(B) = eo E +0, B + c, B? + cz BF + c, Bt + ¢; VB (30) 


bestimmen sich aus dem folgenden Gleichungssystem (3 zweifache charakteristische 
Zahien). . 
g(0) = 1, g'(0) = 1, 
g(t PD) = e?, Feo) =e, (31) 
g(— +®) =e '?, C\ ee 1@) = e-*®, 

mit 

g(x) = cy) +, 2 + 6,27 +¢,2° +c, 24 +6; 2 

g'(z) =¢c, + 2c,% 4+ 3c, 22 + 4c, 2° + 5c; 24 

Aus (31) erhalt man 


C= c= 1, 


} 
nsapsate ee ro [4 (cos 6 — 1) + Osin GO}, e 
Cs sar [P (cos ® — 1) — 5 (sind — G)}, | (32) 
= — — [2 (cos 6 — 1) + @sin PD], | 
C= zor lO (cos S — 1) — 3(sin d— 6)}. | 


Durch Einsetzen der entsprechenden Exponentenmatrizen © in (29) und (30) 
mit den Koeffizienten c, aus (32) erhalt man die Ubertragungsmatrizen (34a), (34b), 
(34c) und (34d). In (34e) ist die Ubertragungsmatrix fiir einen masselosen starren 
Kreisbogen angegeben; in (34f) die Ubertragungsmatrix (Eckmatrix) fiir den Fall, 
dai zwei Stabe unter einem Winkel @ eine Kcke bilden (Abb. 4), die etwa durch 
den Grenziibergang o —> 0 aus der Ubertragungsmatrix (34e) gewonnen werden kann. 
Die in (34a) und (34b) vorkommenden Funktionen F,(®) und ihre Reihenentwick- 
lungen sind: 


FD) == (®cos + sin) — 6-2 4 3 4 or 4 

F,(®) = Osin® + cos@—1= 5,8 — 4 O42 os 7 os 4 

F,(®) = = (sin ® — 6 cos) = 4, Ps = ae : = 01 P+... 

F,(®) = 2sin  — Gcos  — == Be 3 5p +9 or e+... = 
F,(®) = (2 — 2 cos ® — @ sin®) = fe ee a Ovi... “#59 
F,(®) = + (26 + Ocos @ — 3sin ) = ere ae eT ee), J 
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Ubertragungsmatrix U,. Masseloser Kreisbogen. 


af o F,(®) | eee. Gai 
cos ® ; u sin® |o0(1—cos®) |ag(@—sind) 2 
+a F(®) | tae FG) 
0 0 —sin B® 
Cl) Vr 
o sin ® a 0? (1— cos ®) Ree) (34a) 
+aoF,(®) ‘F(®) 
ao? (®— sin D) 0 ee a — ao® ao? (1 — cos D) 
o (1 — cos ®) Ont} 0 / 1 osin ® 
sin @ Oa: Ona a 0 cos ® 
Ubertragungsmatrix U,. Starrer Kreisbogen. 
cos @ 0 sin @ o(1—cos®) | 0 0 
—pw*? oe Pcos® pe ee ate) 0|—sin® 
—sin ® 0 cos B o sin ® 0 0 
0 0 0 1 | 0 (34b) 
2 2 | @ Pend) [oo a pam e 
—2uo oe F(®) e (1 — cos ) Uo e o( ) + uo oe F,() esi 
poo Psin D snd | wooDcos® | pore2F(®) | 0) cos® 
Ubertragungsmatrix U,,. Einzelmasse. 
(34c¢) 
Ubertragungsmatrix U,. Hinzelelastizitat. 
1 ef 6. 00,20 
Opi O. 0. 0 
Os 220 be 5 Oe (34 d) 
Cet Oot la eed 
my 02 Ua 0 
| (eet 0 OyeOy vet 
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Ubertragungsmatrix WU). Starrer masseloser Kreisbogen. 


cos @ 0 | sn® j|e(1—cos®)| 0 0 
0 cos @ 0 0 0 — sin @ 
SE ET 
o sin 0 0 
(34e) 
1 0 0 
0 1 o sin ® 
0 0 cos ® 
| O 0 ORS] 
0 0 |—sin@® 
0 0 0 
(34f) 
1 0 0 
0 1 0 
| 0 0 | cos® 


Fall B: Schwingungen senkrecht zur Ebene der Zentrallinie. 


Entsprechend wie in Fall A erhalt man hier mit der Koeffizientenmatrix % des 
Differentialgleichungssystems (20) die Ubertragungsmatrix fiir einen homogenen 


Kreisbogen (— = k = konst., « = konst., 6 = konst., ¢ = konst., Zentriwinkel ®) 
zunachst in der Form: 


U = exp %U-s = exp @ (35) 
mit 
s=90= 2 und A=porb ot. (36) 
Die charakteristische Gleichung der Matrix % @ ist 
det (E A — Ao) = Af 4 Atle A“* +Agiel ys 
da? C 2,2 C tye ty? c © ig? 
PALA = Ae at ae ga ee 


Aus den bereits im Fall A genannten Griinden verzichten wir auch hier auf die 
Herleitung der allgemeinen Ubertragungsmatrix fiir den homogenen Kreisbogen. 
Die Ubertragungsmatrizen fiir Sonderfille erhalt man auf dem gleichen Wege wie 
im Fall A durch Grenziiberginge in (35) und anschlieB8ende Umformung mit Hilfe 
des Ersatzpolynoms. Dabei haben die Exponentenmatrizen 8 in einander entspre- 
chenden Sonderfallen dieselben charakteristischen Gleichungen und charakteristischen 
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Zahlen, so daB man die Koeffizienten der Ersatzpolynome nicht neu zu berechnen 
braucht. 

Wir kénnen uns daher hier damit begniigen, die endgiiltigen Ubertragungsmatrizen 
anzugeben. Die in diesen Ubertragungsmatrizen enthaltenen Funktionen F,(®) sind 
bereits in (33) erklart. 

Die Hinzelelastizitaten e, und e, sind auf die gleiche Weise wie die Einzelelastizitiat e 
im Fall A eingefiihrt. Die Ubertragungsmatrix (38d) wird aus (35) durch den Grenz- 
tbergang 9 Db ->e,, oDc —-e, mit © +0 gewonnen. 


Ubertragungsmatrix 11,. Masseloser Kreisbogen. 


b+e | 
0} sn® | e—,—@sin@ | oe (b +c) F,(G) 
0 0 sin @ o (1 — cos ®) 
le @sin® | be F,(%) 
1| esin® | 008 F,(®) 
—¢ oe F.(®) vaaiiae (38a) 
bow dO) | tea gs 
ot" @sino 0| cos® ae, ye ae 
SUEY) | 6 FO) 
; cos @ o sin ® 
Oar 0 1 


Ubertragungsmatrix U,. Starrer Kreisbogen. 


cos ® 0 0 sin ® On| 0 a 
| po oi, F,(®) Uo oe) 
+ po ei2Fs(®) | cos® |uw2o2(P—sin®) meee (a +%y)") sin® | 9 (1 — cos ®) 
— pw 0° - F,(P) + Osin® 
— (1 — cos ®) cae ‘ gore <= | (38 b) 
—sin @ 0 0 cos @ : : 
—po? o F;(P) po 01,2 F,(®) 
+ po? 9 (ta + 4y)"| _ sind |yoro(1—cos®)| +uore(ig+  cos® osin ® 
+ Osin® + 0°) F,(B) 
— po o (® — sin ®) 0 uw o® uw @ (1 — cos ®) Oe 1 = 


Ubertragungsmatrix U,,. Hinzelmasse. 
[ 1 OF O 0 


_moiz1! 0 0 0 


(38c¢) 
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Ubertragungsmatrix ,. Hinzelelastizitat. 


0 0 
OOO 
Beda i (38d) 
Ye 
Qin die oO 
Cena 
Ubertragungsmatrix U,). Starrer masseloser Kreisbogen 
sin @ 0 0 
0 sin® 0 (1 — cos ®) 
osin ® 0 0 
: ee (38e) 
cos @ 0 0 
0 cos @ o sin ® 
0 0 1 
Ubertragungsmatrix U,. Eckmatrix. 
[ cos @ pg fe re sn® | 0 OLasetel 
0 0 
2 : (38f) 
0 0 
0 0 
0 | 


B. Reziproke Ubertragungsmatrizen 


In den Anwendungen werden auch die Reziproken der Ubertragungsmatrizen 
benotigt. Fiir die in dem vorhergehenden Abschnitt aufgestellten Matrizen kénnen 
die Kehrmatrizen wegen der besonderen Bauart leicht durch nachfolgende Uber- 
legungen hergeleitet werden. 


Die reziproke Matrix der Ubertragungsmatrix ll = e%s = ele? fiir ein Kreis- 
bogenstiick ist 


Ut ane tee (39) 


(s. z. B. [3], S. 120). Da in den im Abschnitt A betrachteten Fallen die Elemente 
der Koeffizientenmatrizen % unabhingig von © sind, erhalt man die reziproke Matrix 
in der Gestalt einer quadratischen Matrix U-1 aus der zugehérigen quadratischen 
Ubertragungsmatrix U dadurch, daB man in dieser ® durch — @ ersetzt. Nun sind 
die Elemente der Ubertragungsmatrix fiir ein homogenes Kreisbogenstiick, wie man 
zeigen kann, entweder gerade oder ungerade Funktionen von @, die sich in den 
beiden betrachteten Schwingungsfallen wie folgt verteilen. 
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gujugug 
ugigugw 
gujugug 
(41) 


g Elemente der Ubertragungsmatrix, die gerade Funktionen von @ sind. 
u Elemente, die ungerade Funktionen von © sind. 


Man erhalt also die zugehGrige reziproke Matrix aus der jeweiligen Ubertragungs- 
matrix bereits dadurch, daB man in dieser bei denjenigen Elementen das Vorzeichen 
umkehrt, deren Stelle in der Matrix mit ,,w‘‘ gekennzeichnet ist. Diese Vorschrift 
und das Schema (40) bzw. (41) gelten auch fiir die Bildung der Kehrmatrizen in den 
Sonderfallen (34a bis f) bzw. (38a bis f). 


CG. Ubertragungsmatrix fiir einen stiickweise aus Kreisbégen 
zusammengesetzten Stabzug 


Fir das folgende beachten wir, daB wir den Zustandsvektor 3 bei der Behandlung 
der eben gekriimmten Stabe auf das rechtwinklige x, y, z-Koordinatensystem bezogen 


Abb. 5 Abb. 6 


haben, bei dem die positiven Richtungen der x, y, z-Achsen der Reihe nach in Richtung 
des Tangentenvektors, Hauptnormalenvektors und des Binormalenvektors weisen 
(vgl. Abschnitt III und die dort getroffenen Vereinbarungen). Dieses System nennen 
wir kurz das Normalsystem. 

Wir beschaftigen uns zuerst mit der Aufstellung der Ubertragungsmatrix fiir 
ein Stabstiick, dessen Zentrallinie sich aus zwei Kreisbogenstiicken zusammensetzt, 
das heiBt wir fragen nach der Matrix, die den Zustandsvektor an einem Ende eines 
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solchen Stabes auf sein anderes Ende transformiert. Dabei ist bei Schwingungen in 
der Ebene der Zentrallinie der Zustandsvektor 


3 se {u, Qe Vv; P; pti pes = Ort (42) 
bei Schwingungen senkrecht zur Ebene der Zentrallinie 
3 = {6, M,, — W, Y) iM, Q.}. (43) 


Wir betrachten zwei Falle, die in den Abb. 5 und 6 dargestellt sind und sich dadurch 
unterscheiden, da8 beim Uberschreiten der Stelle 2 im 1. Fall (Abb. 5) der Krimmungs- 
sinn der Zentrallinie beibehalten wird, wahrend er sich im 2. Fall (Abb. 6) umkehrt. 


1. Fall (Abb. 5): Hier gelten die Beziehungen 
32,2= Ur a1, (44) 
3 = Unt: 32,r (45) 
mit den Ubertragungsmatrizen Uy, Uy; nach A. Die Indizes ,,2, I‘ und ,,2, r“‘ weisen 
darauf hin, daB sich der Zustandsvektor auf eine Stelle unmittelbar links bzw. un- 


mittelbar rechts der Stelle 2 bezieht. Fordert man fiir die Komponenten des Zustands- 
vektors aus mechanischen Griinden Stetigkeit an der Stelle 2, so gilt 


32,r—= €+ 32,  (€ = Hinheitsmatrix). (46) 
Aus den Gl. (44), (45) und (46) folgt nun durch Einsetzen 
gs = Un © Up jy = Un Ung. = U- 4. (47) 


U = Uz: Uz ist somit die gesuchte Matrix, die den Zustandsvektor von der Stelle 1 
auf die Stelle 3 tibertragt. 


2. Fall (Abb. 6): Hier gelten zunachst die Beziehungen (44) und (45). Fordert 
man wieder Stetigkeit des Zustandsvektors an der Stelle 2, so muB bei der Aufstellung 
der Beziehung zwischen 3.,, und 4.,, beriicksichtigt werden, daB sich 3,,, und 3,,, 
auf verschiedene Normalsysteme beziehen. Diese Koordinatensysteme haben an 
der Stelle 2 eine gemeinsame x-Achse und das eine folgt aus dem anderen durch 
Drehung um die x-Achse um 180°. Die Beziehung zwischen 3,,, und 3,,, lautet jetzt 


32,7 = C+ 30.1, (48) 
die sowohl fiir den Zustandsvektor (42) als (43) gilt, 


mit G* = = G1, (49) 
Aus (44), (45) und (48) folgt nun 
ds = Un © Uys jg =U (50) 
mit der gesuchten Ubertragungsmatrix 
U = Uy - G* - Uy. 


Oft ist es zweckmaBiger, den Zustandsvektor auf ein von der Anderung des 
Kriimmungssinnes unabhangiges Koordinatensystem zu beziehen (Abb. 7): die x-Achse 


Beitrag zur Kinetik und Statik ebener Stabwerke mit gekriimmten Staben 211 


weist in Richtung des Tangentenvektors, die y-Achse ist beim Fortschreiten in 
Richtung wachsender Bogenlinge immer nach derselben Seite gerichtet (in Abb. 7 
nach oben) und die z-Achse ergiinzt das System zu einem Rechtssystem. 


Zur Unterscheidung wollen wir auch die auf dieses System bezogenen Zustands- 
vektoren durch die Querstriche unten kennzeichnen. Im Abschnitt I (Abb. 7) ist 
dieses neue System mit dem Normalsystem identisch, so da8 gilt 


Sl Ot Aaa, toast 
und 

32,1 = Uy bi: (51) 
Dagegen miissen wir die Zustandsvektoren 
im Abschnitt II auf das neue Koordinaten- 
system umrechnen. Hs gilt 


32,7 = &* * }2.0 
$3. E* + 35 


(52) Abb. 7 


mit ©* nach (49). Die Beziehung (45) rechnet sich mit Hilfe von (52) auf das neue 
Koordinatensystem wie folgt um 


O45 — Un E*- 4a, 
3s — (G*-1 . Unt ° E* - 30.5 
und wegen &* = (*-1 [vgl. Gl. (49)] 
33 = &* - Ur - C* go, = Un* - 32,, (53) 
mit der im neuen See geltenden Geran cones 
Un? = Gein €. 
Die Produktmatrix U;,;* lautet ausfiihrlich: 


i 101 070.70" 0 Us, Ure) 3 Ure 
0 170 00 0 


Us, Use 


U23 Ura 


(54) 


Sie entsteht aus der Matrix My, einfach dadurch, daB man in dieser bei den Elementen, 
welche die Bedeutung von Kopplungsfaktoren haben, das Vorzeichen umkehrt. 
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Mit der Stetigkeitsforderung 


der =F gat (55) 
und den Beziehungen (51) und (53) folgt 
gs = Un* Ur- jg. = Ua (56) 
Die gesuchte Ubertragungsmatrix lautet hier 
| U = Unt - Uy. 


Nachdem wir gezeigt haben, wie man die Ubertragungsmatrix fiir ein Stabstiick 

aufstellt, dessen Zentrallinie aus zwei Kreisbégen besteht, kénnen wir nun leicht durch 

schrittweises Vorgehen die Uber- 

tragungsmatrix fiir einen beliebi- 

45 d, gen Stabzug berechnen, dessen 

Zentrallinie aus Kreisb6gen zu- 

sammengesetzt ist, wobei biege- 

7 steife Ecken zugelassen sind. Fiir 

6 den Stabzug in Abb. 8 z. B. be- 

rechnet sich die Ubertragungs- 
matrix wie folgt: 

Abb. 8 1. Weg: Man behalt das ur- 

spriingliche Koordinatensystem 

~ (das Normalsystem) bei und benutzt beim Ubergang in das anschlieBende i-te 

Stabstiick je nach Kriimmungsverhalten die Beziehungen 


(byte C- 3c... oder ji.7 GF: bi. te (57) 


2 x 


Man erhalt dann 
aq = U,: © Uy CU, CFU, So, <9, = (58a) 
2. Weg: Man benutzt das oben definierte, von dem Kriimmungsverhalten un- 
abhingige Koordinatensystem. Dann sind fiir die Stababschnitte, in denen die 
y-Achse in Richtung der Hauptnormalen weist, die Ubertragungsmatrizen U1, zu 
nehmen, im anderen Falle 11,*. Man erhilt: 


oG a" US Ue a a ols a H- ba: (58b) 


D. Durchfiihrung des Reduktionsverfahrens bei Stabziigen ohne 
Zwischenbedingungen 


In diesem Abschnitt beschaftigen wir uns mit der Einarbeitung der Randbedin- 
gungen und der Aufstellung der Frequenzgleichung bei den unter C niher gekenn- 
zeichneten einfachen Stabziigen. Fiir die Darstellung beschrinken wir uns auf den 
Fall der freien Schwingungen in der Ebene der Zentrallinie des in Abb. 8 skizzierten 
Stabzuges. Die Beziehung zwischen den Zustandsvektoren an den Randstellen 1 
und 6 dieses Stabzuges ist durch (58) gegeben und lautet ausfiihrlich geschrieben: 


U Ti1 Tr. | Tis Tia Tis 116 U 


Q. Tor Teo} Tes Toa Ve5 126 Q. 


(59) 
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An den Randstellen 1 und 6 seien homogene Randbedingungen vorgegeben. Diese 
sind im allgemeinen so beschaffen, da je drei Komponenten der Zustandsvektoren 25 
und 3, verschwinden. Setzt man die geltenden homogenen Randbedingungen in die 
Beziehung (59) ein, so erhilt man sechs homogene Gleichungen fiir die tibrigen 
Komponenten der Zustandsvektoren an den Randstellen. Drei von diesen Gleichungen 
(die Nullbedingungen fiir drei der Komponenten von 36) ergeben ein homogenes 
Gleichungssystem fiir drei der unbekannten Komponenten von 3, allein. Die Bedingung, 
daB die Koeffizientendeterminante dieses Gleichungssystems, deren Elemente vom 
Parameter w abhingen, im Falle einer nichttrivialen Lésung verschwinden muB 


»|=0 . (60) 


liefert die Bestimmungsgleichung fiir die Eigenfrequenzen. (60) wird als Frequenz- 
gleichung und A(w) als Frequenzdeterminante bezeichnet. 


Ist etwa der Stabzug in Abb. 8 an seinem linken Ende (Stelle 1) eingespannt 
und an seinem anderen Ende (Stelle 6) gelenkig gelagert, so gelten bei Betrachtung 
der Schwingungen in der Ebene der Zentrallinie die Randbedingungen 


a, = 0, 2, = 0; == 0); 
1 1 vA (61) 
esx), UV, = 0, Mo, == 0, 
und die Frequenzdeterminante ist 
Tre N45 116 
A(w) =!T32 35 16 (62) 


| Ts. Vs5 56 


Die praktische Durchfithrung des Reduktionsverfahrens zur Bestimmung von 
Kigenfrequenzen gerader Stabe ist in der Literatur ausfiihrlich beschrieben (vel. 
z. B. [5], [7], [9]). Mit den im A schnitt A bereitgestellten Ubertragungsmatrizen und 
den Ausfiihrungen im Abschnitt C verliuft die Rechnung bei gekriimmten Staben in 
genau entsprechender Weise. Nur in sehr einfachen Fallen ist es méglich, ohne erheb- 
lichen Rechenaufwand den Parameter w in der Rechnung bis zum Schlu8 mitzufiihren 
und somit A(w) geschlossen darzustellen. Im allgemeinen wird man daher die Rechnung 
fiir eine Reihe von fest angenommenen Parameterwerten w aus dem Frequenzbereich, 
fiir den das Schwingungsverhalten des Stabzuges interessiert bzw. in dem Kigen- 
frequenzen zu erwarten sind, jedesmal neu durchfihren (die Elemente der Ubertragungs- 
matrizen sind dann reine Zahlen) und die zugehérigen Funktionswerte A(w) berechnen. 
Die Nulldurchgiinge der Frequenzdeterminante findet man anschlieBend durch Inter- 
polation. Die Berechnung des Matrizenproduktes in (58) erfolgt zweckmaBig in einem 
fiir Matrizenmultiplikation tiblichen Schema (vgl. z. B. [9]). 


Fir wo =o, (w, = Higenfrequenz) besitzt das homogene lineare Gleichungs- 
system, welches zur Frequenzgleichung fiihrte, eine nichttriviale Losung. Durch 
diese sind die Zustandsvektoren an den Randstellen vollkommen bestimmt bis auf 
eine multiplikative willkiirliche Konstante, tiber die man etwa so verfiigen kann, 
daB man einer der Komponenten der Zustandsvektoren einen festen Wert erteilt, 
z. B. gleich ,,eins setzt. Ist nun der Randvektor 3, bekannt, so lassen sich die 
Zustandsvektoren an Zwischenstellen leicht berechnen, wodurch der Verlauf der 
Verschiebungs- und SchnittgroBen festgelegt ist. 
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Fir das betrachtete Beispiel (Abb. 8) folgt 
C2 Uy: bi 


gs = UU 
; (63) 


pe as ee il i, Uy" 3: 


E. Elastische Lager als Rand- und Zwischenbedingungen 


An einer Stelle i sei der Stabzug durch Federn und Drehfedern elastisch gelagert. 
Die Feder- bzw. Drehfederreaktionen lassen sich nach dem Kraft-Dehnungs-Gesetz 
fiir lineare Federn durch die zugehérigen VerschiebungsgréBen ausdriicken. Es gilt 


Q2,i = — Cy Ue, 
M2 = —Co* Vis | (64a) 
Qy..0 = — by M%, 
My,i = — 8° 4;, 
My, 5 = — Cy" Vis | (64b) 
Q2.4 = — Cw" Wi, 


wobei Cinaex die entsprechenden Federkonstanten bedeuten. Die Schnitt- und Ver- 
schiebungsgroBen, die sich auf die Stelle unmittelbar links (rechts) von der Stelle ¢ 
beziehen, kennzeichnen wir durch einen Index ,,/“ (,,7°‘). 


Fall A: Schwingungen in der Ebene der Zentrallinie (Abb. 9a). 


Zwischen den Schnitt- und VerschiebungsgréBen an den Stellen J und r gelten 
hier unter Beriicksichtigung von (64a) die Beziehungen 


Pr = Pr Up = U1; 
Me ilies + Cp Pr A as (ee + Cy* tis (65a) 
Quin = Qe.1 + Co MD Ur = Vy. 


Mit dem Zustandsvektor 3 nach (42) und der Matrix 


J 
SS s 
Cu, Cyr ly < a) ay 


Abb. 9a 


(66a) 


lassen sich die Beziehungen (65a) kurz als eine Matrizengleichung schreiben 


b= 8° at (67) 


Die Matrix § bezeichnen wir als Ubertragungsmatrix fiir ein elastisches Lager. 
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Fall B: Schwingungen senkrecht zur Ebene der Zentrallinie (Abb. 9b). 
Hier gelten unter Beriicksichtigung von (64b) die folgenden Beziehungen 


6, > 61, Yr = Yo 
M,,, = My, + 65° 6, My,» = My + ey, (65 b) 
W, = Wi, Qe.r = 92,1 + Cw? Wie 
Diese lassen sich wieder in der Gestalt (67) schreiben, wobei jetzt 3 nach (43) und 
roe 


Se Qo => 


zu setzen ist. 

Elastische Lager als Rand- und Zwischenbedingungen bringen somit fiir die Durch- 
fiihrung des Verfahrens keine Schwierigkeiten mit sich. Sie lassen sich nach den 
obigen Ausfiihrungen in genau entsprechender Weise durch Ubertragungsmatrizen § 
erfassen wie Stababschnitte mit bestimmten Massen- und Steifigkeitseigenschaften 
durch die Ubertragungsmatrizen MU. 


F. Vorgehensweise bei allgemeinen inneren Bedingungen 


Sind auch andere innere Bedingungen vorhanden, wie z. B. Gelenke oder feste 
Zwischenlager, so sind die Verhialtnisse nicht mebr so einfach. An den Stellen des 
Stabzuges, an denen solche innere Bedingungen vorgeschrieben sind, treten neue 
unbekannte GréBen hinzu, so z. B. bei Lagern die Lagerreaktionen. Die Anzahl der 
hinzukommenden Unbekannten stimmt jedoch mit der Anzahl der inneren Bedin- 
gungen iiberein, so daS man insgesamt ebenso viele homogene Gleichungen wie 
Unbekannte hat. Zur Erlauterung sind nachstehend einige Beispiele fiir innere 
Bedingungen bei Betrachtung der Stabschwingungen in der Ebene der Zentrallinie 


angegeben. 
= SS 
x VF . 
—- hee chee 
Abb. 10 


a) Tangential verschiebliches Lager (Abb. 10a). 
Bedingung: v = 0. 
Hinzukommende Unbekannte: Lagerkraft S,. 
b) Gerbergelenk (Abb. 10b). 
Bedingung: M, = 0. 
Hinzukommende Unbekannte: Winkel «. 
c) Unverschiebliches (festes) Lager (Abb. 10c). 
Bedingungen: u=— 0, v = 0. 
Hinzukommende Unbekannte: Lagerkrafte S, und S,. 
15* 
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An einer Stelle koénnen héchstens drei mechanisch sinnvolle Bedingungen vor- 
geschrieben sein. Sind an einer inneren Stelle 7 des Stabzuges genau drei Bedingungen 
vorgeschrieben (z. B. bei Einspannung an der Stelle 7 und Betrachtung der Schwin- 
gungen in der Ebene der Zentrallinie: u; = 0, v; = 0, y; = 0), so hat man schon bei 
der Rechnung bis zu dieser Stelle geniigend Gleichungen zur Bestimmung der bis 
dahin aufgetretenen Unbekannten. Dies bedeutet, daB das Schwingungsverhalten 
der Teile des Stabzuges links und rechts von der Stelle 1 voneinander unabhangig 
ist und jedes der beiden Teile gesondert fiir sich behandelt werden kann. Von Interesse 
ist somit nur der Fall, da an einer inneren Stelle eine oder zwei Bedingungen zu 
erfiillen sind. 

Fiir die praktische Durchfiihrung der Rechnung bei Stabziigen mit inneren Be- 
dingungen bestehen grundsitzlich zwei Méglichkeiten: Entweder eliminiert man die 
hinzukommenden Unbekannten im Laufe der Rechnung sofort bei ihrem Auftreten, 
so daB man zum SchluB ein homogenes Gleichungssystem nur fiir die Unbekannten 
an den Randstellen hat, oder man nimmt sie in der Rechnung bis zum SchluB mit 
und erhalt dann ein erweitertes homogenes Gleichungssystem fiir die Unbekannten 
an den Rand- und Zwischenstellen (vgl. [9], Abschnitt 4). In den beiden Fallen liefert 
die gleich Null gesetzte Koeffizientendeterminante des Gleichungssystems die 
Frequenzgleichung. 

Bei allgemeinen (verzweigten) Stabwerken zerlegt man das Stabwerk durch 
Schnitte in einzelne Stabziige und fiihrt die Schnitt- und Verschiebungsgréfen an 
den Schnittstellen als Unbekannte ein. Die Stetigkeitsforderungen fiir die Schnitt- 
und VerschiebungsgréBen an den Schnittstellen liefern die zusatzlichen Gleichungen 
fir die Bestimmung dieser Unbekannten. 

Fir spezielle Stabwerke, wie z. B. fiir einen durchlaufenden Stabzug auf festen 
Stiitzen, lassen sich zweckmaBige Rechenschemata in ahnlicher Weise aufstellen, wie 
es fiir gerade Durchlauftraiger in [10] ausgefiihrt ist. 


G. Behandlung beliebiger ebener Stabwerke. Wahl geeigneter Ersatz- 
systeme 


Nach dem in den vorhergehenden Abschnitten entwickelten Verfahren lassen sich 
freie Schwingungen von beliebigen ebenen Stabwerken mit beliebiger Masse- und 
Steifigkeitsbelegung naherungsweise berechnen, wenn man dem wirklichen System 
ein Ersatzsystem zuordnet. Dieses Ersatzsystem mu so beschaffen sein, daS man 
bei seiner Durchrechnung mit den im Abschnitt A bereitgestellten Ubertragungs- 
matrizen fiir Sonderfalle auskommt. Bei der Aufstellung eines Ersatzsystems miissen 
daher die Zentrallinien der einzelnen Stabziige des Stabwerkes stiickweise durch 
Kreisbégen angenadhert und die gegebene Masse- und Steifigkeitsverteilung durch 
eine stiickweise konstante oder diskrete ersetzt werden. 

Von den verschiedenen Moéglichkeiten, die Masse-Steifigkeits-Eigenschaften des 
gegebenen Systems zu erfassen, sind die folgenden drei die wichtigsten: Entweder 
zieht man die Massebelegung in diskrete Einzelmassen zusammen und nihert gleich- 
zeitig die Steifigkeitsbelegung durch eine stiickweise konstante an (M-Methode), 
oder umgekehrt, man ersetzt die Massebelegung durch eine stiickweise konstante 
und die Steifigkeitsbelegung durch Einzelelastizitaten (S-Methode), schlieBlich kann 
man sowohl die Masse- als auch die Steifigkeitsbelegung durch eine diskrete annahern 
(MS-Methode). Die nach der M- und S-Methode gebildeten Ersatzsysteme erfassen 
das wirkliche System etwa mit gleicher Genauigkeit. Im allgemeinen wird man 
jedoch bei Schwingungsuntersuchungen die nach der M-Methode aufgestellten Ersatz- 
systeme bevorzugen, weil sich dabei meistens rechnerische Vorteile gegeniiber den 
Ersatzsystemen nach der S-Methode ergeben. Das Ersatzsystem nach der MS-Methode 
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besteht aus Einzelmassen und Einzelelastizitaten, die durch starre masselose Kreis- 
bogenstiicke miteinander verbunden sind, und hat aus erklarlichen Griinden eine 
geringere Genauigkeit als die Ersatzsysteme nach den beiden anderen Methoden. 
Man hat aber bei der Rechnung mit einem solchen Ersatzsystem den Vorteil der 
besonders einfach gebauten Ubertragungsmatrizen. Selbstverstandlich kénnen bei 
der Festlegung des Ersatzsystems auch alle drei Methoden nebeneinander verwendet 
werden. 

Da das Ersatzsystem mit dem Reduktionsverfahren exakt erfaBt wird, liegt der 
Naherungscharakter des Ergebnisses begriindet nur in der Ersetzung des wirklichen 
Systems durch ein Ersatzsystem. Die erforderliche Genauigkeit der Ergebnisse laBt 
sich durch eine entsprechend feine Gliederung des Ersatzsystems stets erreichen. 


H. Beispiel 


Wir betrachten einen kreisbogenférmigen homogenen Stab mit dem Zentri- 
winkel « bei verschiedenen Lagerbedingungen an seinen Enden. Fiir die freien 
Schwingungen dieses Stabes in seiner Ebene bei den jeweils geltenden noch naher 


Abb. 11 


festzulegenden Randbedingungen sollen die beiden niedrigsten Eigenfrequenzen @, 
und w, berechnet werden. Die gesuchten Eigenfrequenzen haben die Gestalt (vgl. 
ZB. {h}) 
_ Oe, f) 5 
NE ae 
Kennt man somit eine Eigenfrequenz w; fiir einen speziellen Stab, so konnen nach 
der obigen Beziehung die entsprechenden Higenfrequenzen auch fiir Stabe mit anderen 
Werten von mw, a und @ bei festem « und f berechnet werden. 
Es sollen hier der Rechnungsgang und die Ergebnisse fiir das spezielle Beispiel 


a = = 60°; f=a03 pe ligkes | ft | ee 
mitgeteilt werden. 
Zur Durchfiihrung des Reduktionsverfahrens wurde fiir den homogenen Kreis- 
bogenstab ein Ersatzsystem nach der M-Methode (vgl. Abschnitt G) mit aquidistant 
angeordneten Hinzelmassen m gewihlt. Das Verfahren wurde fir zwei solche Ersatz- 
systeme durchgefiihrt, die aus den Abb. lla und 11b zu ersehen sind. 


Beim Ersatzsystem a ist (Abb. lla) ® = sth 15° und die Einzelmassen haben 


. It 
mit «= 1 die Werte m =e an den inneren Stellen des Stabes, m —-— an den 


24 
Randstellen. 


9 Beziiglich der Bezeichnungen vgl. Abschnitt III. 
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Beim Ersatzsystem b ist (Abb. 11b) ® = si = 10° i = an den inneren 


Stellen des Stabes, m = an an den Randstellen. 


Bei der Rechnung mit dem Ersatzsystem b wurde die Rotadinstraghatt der Stab- 
elemente beriicksichtigt, waihrend sie bei a vernachlassigt wurde (2, = 0). 


Mit den oben angegebenen Daten wurden die Ubertragungsmatrizen U, und U, 
nach (34a) und (34c) bzw. die fiir die numerische Rechnung zweckmafigeren 
Matrizen IN = U,,- Ul, zahlenmaig aufgestellt. Aus diesen wurden dann fiir sieben 
Parameterwerte w — 20, 25, 30, 35, 40, 45, 50 die Ubertragungsmatrizen fiir den 
ganzen Stab nach Abschnitt C berechnet und anaehiiefend die Frequenzdeterminanten 
fiir folgende Randbedingungen gebildet: 


1. gelenkig gelagert — gelenkig gelagert, 
2. eingespannt — eingespannt, 
3. eingespannt — gelenkig gelagert. 


Mit Hilfe der berechneten Frequenzdeterminanten wurden dann die Frequenz- 
kurven gezeichnet und die gesuchten EHigenfrequenzen als Nulldurchgange der Frequenz- 
kurven ermittelt. In der nachstehenden Tabelle sind die zahlenmaBigen Ergebnisse 
zusammengestellt. Sie enthalt in den Fallen 1 und 2 zum Vergleich auch die mit 
Hilfe des Ritzschen Verfahrens erhaltenen Resultate (vgl. [1], Zahlentafeln 18 bis 20). 


Ritzsches Verfahren .............. 19°80 | 33°63 25°99 53°76 
Reduktionsverf., Ersatzsystem a .. 19°8 32°5 25°5 48°9 21°5 40°7 
Reduktionsverf., Ersatzsystem b 19°7 | 31°6 25°4 48°6 21°2 40°2 


Fiir die Randbedingungen 1 und 2 stimmen die nach dem Reduktionsverfahren 
berechneten Werte fiir die erste Eigenfrequenz wm, mit denen nach dem Ritzschen 
Verfahren gut tiberein. Dagegen treten in den Ergebnissen fiir w, in beiden Fallen 
Abweichungen auf, und zwar liefert das Ritzsche Verfahren die gréBeren Werte. 
Diese Abweichungen sind darauf zuriickzufiihren, daB bei der Bestimmung von w, 
nach dem Ritzschen Verfahren die Dehnung der Stabachse vernachlassigt, wihrend 
sie beim Reduktionsverfahren beriicksichtigt ist. Eine Kontrollrechnung mit dem 
Ersatzsystem a fiir den Fall des beiderseits eingespannten Stabes, bei der in den 
Ubertragungsmatrizen der Einflu8B der Dehnung vernachlassigt wurde (f = 0), ergab 


W, = 52°4. 


Dieser Wert stimmt mit dem Ergebnis nach dem Ritzschen Verfahren erheblich 
besser tiberein. 


V. Erzwungene Schwingungen eben gekriimmter Stibe 


Wir betrachten einen Stabzug, der durch periodisch wirkende Krafte zu er- 
zwungenen Schwingungen angeregt wird. Wegen des Superpositionsprinzips geniigt 
es, nur den Fall zu betrachten, daB alle Erregerkrafte bzw. Momente harmonisch sind 
und die gleiche Frequenz und Phasenlage besitzen, also die Gestalt’ P sin Q t haben. 

Fiir den Sonderfall der eben gekriimmten Stibe (y = 0; ee Oe re at ee 

o b) 2 2 b) 


k, = 0) erhaélt man aus dem allgemeinen Gleichungssystem (13) mit den Tragheits- 
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widerstinden nach (14) fiir die beiden Schwingungsarten folgende Bewegungs- 
gleichungen: 


Fall A: Schwingungen in der Ebene der Zentrallinie. 


w= Q. +k», / 
Qo’ =kQ, + uti — pe sin Qt, 
Bee Mint (68a) 
yp =aM,, 
M,’ = —Q, + ni2 @—d, sin Qt, 
Q,’ = —kQ. +n i —p, sin Qt. 
Fall B: Schwingungen senkrecht zur Ebene der Zentrallinie. 
d=cM,+ky, 
M,/ =k M, +pi2 6— d, sin Qt, 
w' = — 9, 
y =—kd+bM,, bid 
M, =—kM,+9Q, +ui?~y —d,sin Qt, 
Q, = uw — p, sin Qt. 
Mit dem Lésungsansatz 
u(x, t) = u(x) sin QF, 
Q.(x, t) = Q,(x) sin Qt, (69) 


oe eee eee wee eee eee eee 


erhalten wir aus den Gl. (68) gewodhnliche Differentialgleichungssysteme, die bei 
Verwendung der Matrizenschreibweise lauten: 


uU U 0 
QO. Q. —De 
v l= -| v t+] 0 (70a) 
p P 0 
Ms M, ze z 
zQY —Q, Dy 
6 0 
M,; —d, 
— Ww" = : , (70 b) 
y" 9 
M, —dy 
Q, = Bs 


Sie haben die Gestalt 
Aires Wieraete B, (71) 
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Sind die Elemente der Matrix 1 konstant, so ist die allgemeine Lésung von (71) gegeben 
durch (vgl. z. B. [3], 8S. 121): 


3(s) = e%# 3 + e%# | e-Xe E (a) do 
0 


oder 
,= is go +Y (72) 
mit 
U = eds (73) 
und 
i ells | eX £ (a) do, (74) 


0 


wobei 11 = e%* die von der Behandlung der freien Schwingungen her bekannte, 

dort als Ubertragungsmatrix bezeichnete Matrix ist. Die unter dem Integralzeichen auf- 
tretende Matrix e-% * ist als Kehrmatrix zu e% *nach Abschnitt B 
leicht zu bilden. Die Beziehung (72) stellt den Zusammen- 
hang zwischen den Zustandsvektoren 3 an der Stelles = 0 

7 und einer beliebigen Stelle s her. Besteht der Stabzug aus 
zwei aneinanderstoBenden Abschnitten, in denen jeweils 
verschiedene Differentialgleichungssysteme gelten (vgl. die 
Skizze in Abb. 12), so ergibt sich zunachst 


2 = Us: a1 + F t,.) 
3s = Us 3s Bi J 
Fordert man Stetigkeit fiir den Zustandsvektor an der Stelle 2, so folgt fiir den 
Zusammenhang zwischen den Zustandsvektoren an den Stellen 1 und 3 aus den 
Gl. (75) die Beziehung 
3 = Us*(Ui- 31 +01) + te = Up Uy: 31 + Uy ty + Te. (76) 
Fiir die praktische Rechnung ist es vorteilhaft, unter Hinzunahme einer weiteren 
identischen Zeile der Gl. (72) die folgende Gestalt zu geben: 


2 


Abb. 12 (75) 


3 Wore 0 
| 1 od tela ee (77) 
6S u 30 


U soll ,,erweiterte Ubertragungsmatrix‘‘ und 3 ,,erweiterter Zustandsvektor“ heiGen. 
Durch Ubergang zur Komponentenschreibweise iiberzeugt man sich leicht, daB die 
Gl. (72) und (77) dieselben linearen Beziehungen zwischen den Komponenten der 
Zustandsvektoren darstellen. Wenn man nun mit Hilfe der Darstellung (77) die 
Ausdriicke (75) umschreibt und die Beziehung zwischen 3, und 3, einfach durch 


Multiplikation der siebenreihigen Matrizen 11, und ll, erklart, also durch 


poe oe van (78) 
oder ausfiihrlicher geschrieben 


33 U, [te U, i 31 Uy: U, | U,-r, + re val 
= Caer ese ‘| oo t= 7 (79) 


1 0 Fl O Fl 1 


so tiberzeugt man sich, da in der Tat die Beziehungen (78) und (76) iibereinstimmen. 
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Wir bemerken, da8 die Produktmatrix §$ dieselbe Bauart wie die Faktoren i, und 
U, hat. 
Besteht der Stabzug aus mehr als zwei Abschnitten, so berechnet sich die Be- 


ziehung zwischen den Zustandsvektoren an seinen Enden, z. B. 3, und 3, beim Stab- 
zug in Abb. 8 nach 


jo oll, - Ug Mos Uys Uy: 9. Re ay (80) 
wobei die Produktmatrix 9 von gleicher Bauart wie die von U, ist. Greifen in einem 


Abschnitt keine Erregerkrifte an, so sind in der zugehérigen Ubertragungsmatrix U, 
die Komponenten des Vektors r; gleich Null zu setzen. 


Bei Betrachtung der Schwingungen in der Ebene der Zentrallinie z. B. lautet 
die Beziehung (80) ausfithrlich 


Ug = Ty My + 1290.4 + 11301 + 1a Pr + 15 Ma, + 116 (—Qy.1) +112 
Qe. 6 = Toi Uy + 22 Qe,1 +123 %1 + 14 Pr + 5 Me. 1 + 126 (—Qy.1) + 1223 
Ue = Ts1 Uy +1 32Q9e,1 + 19501 +134 %1 +135 Me,1 +126 (—Q,1) +222 
Po = Tar Uy + 1429.1 +143 % + Tae Pr + 145 Mz. + 146 (—Qy.1) +142, (81) 
Mag =1%51 Ua + 152 90,.1 + 153% +154 91 +155 Mz +156 (—Vy.1) +157» 
—Qy.6 = er M1 + 162 Qe.1 +763 U1 + Tea Pr + 165 Me,1 +166 (—Vy.1) + Ter 
[=r 


Ist der Stabzug etwa an seinem linken Ende (Stelle 1) eingespannt und an seinem 
anderen Ende (Stelle 6) gelenkig gelagert, so gelten die Randbedingungen 


u,= 0; v, = 0; 7, = 9; | 
Uze=9; v=90; M,,,=09. J 


Setzt man die Randbedingungen (82) in die Beziehungen (81) ein, so erhalt man 
das inhomogene Gleichungssystem 


0 = 1y2Qe.1 +115 Ma. +116 (—@v.1) 

Qe,¢6 = 122Qe,1 +125 Mz.1 + 126 (—Qy,1) 
0 = 139 Qe,1 +135 Mz,1 +136 (—Qy.1) +137 

) 

) 


(82) 


Pe = 142 Qe,1 + Tap Mss 1 Tee ( Qy.1 
O == 55 Vas 55 tle. oe ee (— 1) 
—Qy,6 = 762 V2.1 +15 Mie fo — @y,1) tater 


aus dem man die noch unbekannten Komponenten der Zustandsvektoren an den Rand- 
stellen bei nichtverschwindender Determinante 


T12 "15 £4 
A= |132 T35 Ts6 (84) 


Ts2 "55 "s6 
bestimmen kann. 

Mit dem nun bekannten Zustandsvektor 3, lassen sich die Zustandsvektoren an 
beliebigen Zwischenstellen leicht berechnen und somit der Verlauf der Verschiebungs- 
und SchnittgréBen ermitteln. 

In der Praxis kommt oft der Fall direkt vor, daB die Erregung durch periodische 
Einzelkrifte eingeleitet wird. Ferner kommt man auf diesen Fall bei der Aufstellung 
eines Ersatzsystems, wenn man kontinuierlich verteilte Erregerkrafte durch punkt- 
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formig angreifende ersetzt, ihnlich wie man fir eine kontinuierliche Masseverteilung 
eine Einzelmassenbelegung waihlt. Die Ubertragungsmatrix fiir den Ubergang tiber 
eine Stelle, an der Einzelerregerkriifte angreifen, erhilt man durch den Grenziibergang 
n (72) bzw. (77) mit s +0 und gleichzeitig: 


peo PRP, de), 5. 
bei Betrachtung der Schwingungen in der Ebene der Zentrallinie, bzw. 
d,s =De djs Dy ps 7, 


bei Betrachtung der Schwingungen senkrecht zur Ebene der Zentrallinie. Dabei 
bedeuten P,, D,,... die Amplituden der Erregerkrafte. 

Fir die gesuchten erweiterten Ubertragungsmatrizen in den beiden betrachteten 
Schwingungsfallen folgt durch diesen Grenziibergang 


(85a) 


(85b) 


0: 0).0.20%,00/0 jpeyiliiast 


Zum Schlu& bemerken wir noch, daB die resultierende erweiterte Ubertragungs- 


matrix §t in (80) den folgenden Aufbau besitzt: 


Ros Us° Uy- Ug -U,- Uy 


(86) 
mit 
t= Us Uy: Uy: a> t, Us: Uy> Uy- tb Uys Uy: ty US ey ty. (87) 

Will man das Verhalten eines Stabwerkes fiir verschiedene Gruppen von Erreger- 
kraften gleicher Erregerfrequenz 2 untersuchen, so braucht man die Produkten- 
matrix U;-U,-U;-U,- Ul, nur einmal zu berechnen. Bei verschiedenen Erregerkraften 
andern sich im Gleichungssystem (83) lediglich die rechten Seiten, so da nur diese 
nach (87) jeweils neu zu berechnen sind. Ferner erkennt man dabei, da® die Deter- 
minante (84) genau die von den freien Schwingungen her bekannte Frequenzdeter- 
minante ist. Ihr Verschwinden bedeutet, daB die Erregerfrequenz Q mit einer Higen- 


frequenz w, zusammenfallt und das Gleichungssystem (83) dann entweder tiberhaupt 
keine Lésung besitzt oder nicht eindeutig lésbar ist. 
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VI. Statik eben gekriimmter Stabe 


Da die Statik als ein Sonderfall der Dynamik aufgefaBt werden kann, kénnen wir 
alle Beziehungen, die wir bei der Behandlung der erzwungenen Schwingungen auf- 
gestellt haben, sofort auf den statischen Fall iibertragen, indem wir in diesen alle Aus- 
driicke, die von den Tragheitswiderstanden herrithren, gleich Null (rein formal also 
22 = 0) setzen. 


Die Statik der eben gekriimmten Stabe wird durch folgende Differentialgleichungen 
beherrscht, die aus (13) fiir den Sonderfall 4 = 0; + 0 br a OO hae O) 
folgen: 

Fall A: Beanspruchung in der Ebene der Zentrallinie. 


w’ | 0 u 0 
Q. Q» — Pz 
vy |= Vat a1 0 (88a) 
Pp P 0 
M, M, —d, 
| Qy/ —Q, Py 
Fall B: Beanspruchung senkrecht zur Ebene der Zentrallinie. 
[8 OAre a RO ban), 0 é 0 
M,/ 0 0 0 M, a, 
sty 0 0 OF BL MeOy 2.0 —w a 0 (88b) 
y —-k 0 0-0 b..0 y 0 
M,/ 0—-k Oper 0a OF OL M, —d, 
0. | 0 0 00 00 Onl le 


Diese Gleichungen haben die Form (71), und ihre allgemeine Lésung ftir den Fall 


ors é ee 1 : 
eines Kreisbogens mit konstanten SteifigkeitsgréBen = 5, und a bzw. 6 und c ist 


durch (72), (73) und (74) gegeben. Dabei ist die Matrix 11 = e%s fiir den Belastungs- 
fall A in (34a) und fiir den Belastungsfall B in (38a) angegeben. Der Zusammenhang 
zwischen zwei Zustandsvektoren, die sich auf verschiedene Stellen beziehen, kann 
wieder in der Form (77) mit Hilfe einer erweiterten Ubertragungsmatrix 11 geschrieben 
werden. Fiir den Ubergang iiber Stellen, an denen Einzellasten angreifen, erhalt 
man erweiterte Ubertragungsmatrizen, die rein formal mit denen in (85) tiberein- 
stimmen, nur daB jetzt P,, D,,... Einzelkrafte und Einzelmomente bedeuten. 


1 
Ist ein homogenes Kreisbogenstiick (ar = konst., a = konst. bzw. 6b = konst., 


C— konst. mit dem Zentriwinkel ® und dem Radius o durch eine beliebige konti- 
nuierliche Belastung beansprucht, so lautet der fiir die Aufstellung der erweiterten 
Ubertragungsmatrix in diesem Fall bendtigte Vektor 

8 @ 


q = J eM ¥(a)da = @| ee? f(g £) dé, (89) 
0 0 


(s=0P; «=. 6) 
oder ausfiihrlich geschrieben: 
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Fall A: Beanspruchung in der Ebene der Zentrallinie. 


(90a) 


®D 
o?- \ |p. (aff, +aeF,) +d,ao0(f —sin¢) —p, (Sf fsin€—a oF, )|aé, 
0 
@D 
e[[—pscost +p,sin f] dé, 
G c 
e 
o?-\|p.(Safsin¢ —ag?F,) —d.ao(1 — cost) —p, (af tao) Fylde, 
0 ’ 
® 
o* | [peae (6 —sin£) +d,a¢ + pyae (1 — cos ¢)] dé, 
0 
@ 
o°|[—pee (1 — cos 6) —d, — pygsin £] dl, 
0 
@D 
Q -| [p. sin ¢ + p, cos ¢| dé. 
ee 
Fall B: Beanspruchung senkrecht zur Ebene der Zentrallinie. 
© 
ot \|de(eP, —bF,) —d, yb +0) Csin £ +-p,0(b +0) Faldé, 
0 
® 
0 -\[—d, cos ¢ +d,sin ¢ — p, 0 (1 — cos ¢)] dé, 
0 
i 1 
g-|[-d,(b +e) F,—d, (5b Csin€—cF,) +-p,0(bF, —oFe)|ae, 
ane (90 b) 
1 : 1 ihe 
e-\\a.5 b+e)sing +d, (bF, —cF,) —p.o(ybésin s—eF,\|de, 
0 
; 
Q° \ [—d,sin € —d,cos¢ + p, osin 6] dé, 
0 
@D 
-9'| p, de. 
6 


Die Integrale in (90a) und (90b) werden am besten nach einem numerischen oder 
graphischen Verfahren ausgewertet. 


Die Durchfiihrung des Reduktionsverfahrens, insbesondere das Einarbeiten der 
Randbedingungen, erfolgt bei den statischen Berechnungen mit den bekannten Uber- 


tragungsmatrizen U, nun genau so wie bei der Behandlung der erzwungenen 
Schwingungen. 
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Bemerkung zu: , 
Warmeleitprobleme mit zeitlich variabler Ubergangszahl 


Von F. Selig und H. Fieber, Wien 


Herr Paterson von den Imperial Chemical Industries Ltd., Stevenston, machte 
uns freundlicherweise darauf aufmerksam, daB in der im Titel genannten Arbeit? 
die Uberfiihrung eines Problems mit zeitlich veranderlichem Ubergangskoeffizienten h(t) 
in ein Problem mit bewegtem Rand und konstanter Ubergangszahl « durch die 
Transformation 

lees h ue By Sie h ou 1, t= | 42() do 


0 
die Homogenitaét 2. Grades der Laméschen Koeffizienten g; notwendigerweise voraus- 
setzt. Diese Forderung kommt in der Gleichung 


ee Aas 1 h(t f 
96,0") = gi'(E, 0,0) 


zum. Ausdruck, wurde aber nicht explizit ausgesprochen. Um eventuellen MiBverstand- 
nissen zu begegnen, erscheint es uns zweckmafig, hierauf naher einzugehen. 


Sei x = 9(&, 7) eine homogene Funktion 2. Grades, so ist es stets méglich, eine 
homogene Funktion 2. Grades y = y(é, 7) so zu bestimmen, daB die Schar der Kurven 
& = konst. die Schar 7 = konst. orthogonal schneidet. Setzen wir zunachst y als 
bekannt voraus und fassen wir 

a = 9&7), | (1) 
y= En) J 


als Transformation der 2x,y-Ebene auf die é,7-Ebene auf, so sind die Laméschen 


Koeffizienten 
91 = 92 + ve, 


Jo = Pn + Vn 
homogen 2. Grades und 
Giz = PEPn + YE Pn = 9, 


1 Osterr. Ing.-Arch. 11, 37 (1957), wird im folgenden mit I bezeichnet. 
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was unmittelbar einzusehen ist. Zur Bestimmung von »(é, 7) steht uns also das System 
Eye + Yn = 2Y (2) 
PE Pn + Pe Yn = O (3) 


zur Verfiigung. Ist 2 = A(é, 1) homogen 0. Grades, so wird (2) und (3) durch den 
Ansatz 


1 
Ye=AGn Yo=— ZV 
befriedigt, wenn fiir 4 die aus der Integrabilitétsbedingung folgende Gleichung 
; 1 1 
dn Pn + A Qnn = Ga Ae PE — | VEE (4) 
erfiillt ist. Unter Beriicksichtigung der Homogenitatseigenschaften 
p=ro(t) wma 2-a(2 


und der Bezeichnung (...)’ = (5 oc) etna ee = usw., ergibt sich aus (4) fiir 


22 = eine Abelsche Differentialgleichung 
, 22D’ —_2z2@ o” 26— 22D’ + 2G” 
Hae @’ | ei 2|u ae o’ , (5) 


die mittels der Substitution? 


( 22 O’ —2z2@ 
ee 


"am ad (6) 


auf die Form 
vv = [6 (2z—20')@” + 2(22D’—_D)@']v + 42D’ (P’ —z2 6”) (O@—zG@’) (2H—2zG’) (7) 
gebracht werden kann. Hieraus folgt 


+f 


Als elementar integrierbare Falle ergeben sich: 


PEP 5 
N(,n) dg — yi" dy] mit. N# =F v— Loegy (8) 


a) O'=20" ‘fuhrt auk io Ce +9, | 


(= . &% — oe) 7, 
b) O=2z@’  fiihrt auf p= C, én, 
C,2 
p= ze ca, | 


c)} 20—2@’ fihrt auf gy = C177, 
aa =e. 
hierbei bedeuten «, Cj, C, beliebige Konstanten. 


d) Der Fall, da®B der Koeffizient von v in (7) verschwindet, fiihrt auf eine 
Differentialgleichung 2. Ordnung fiir ©, die auf das System 


z—f’ (9) 


File (10) 


* Vel. E. Kamke: Differentialgleichungen, Lésungsmethoden und Lésungen. 
3 Auf eine analoge Transformation fihrt der Fall &’ = 0. 
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zurickgefiihrt werden kann. Die Lésung von (10) lautet: 


ie lim a 


daher ergibt sich aus (9) 


dz 
z—f* 
Die auftretenden Integrationen sind elementar durchfiihrbar. 

Da fiir die Liésung einer partiellen Differentialgleichung die Einfiihrung separabler 
Koordinaten von Vorteil ist, wollen wir noch den Fall diskutieren, da8 die Trans- 
formation (1) eine konforme Abbildung der 2,y-Ebene auf die é,y-Ebene darstellt. 
Neben (2) und (3) miissen dann noch die Cauchy-Riemannschen Differentialgleichungen, 
das hei8t 4 = 1, gelten. Dieser Fall ist in b enthalten und fihrt auf die parabolischen 


Koordinaten «= én, y = 5 (7? — &). 


Aus dem Gang der Rechnung in I ist ersichtlich, daB die Homogenitat 2. Grades 
von g und »y und damit von g,; nicht wesentlich ist, sondern auch homogene g; von 
n-tem Grade durch eine entsprechende Anderung der Zeitskala auf ein Ersatzproblem 
mit konstantem Ubergangskoeffizienten fiihren. Als separable Koordinaten ergeben 
sich Real- und Imaginarteil der logarithmischen Potentiale ¢” = (£ + 77)"; die zu- 
gehorigen Laméschen Koeffizienten sind homogen vom Grad 2 (n — 1). 

SchlieBlich sei noch bemerkt, daB die Homogenitat aller g; unter Umstanden nicht 
erforderlich ist, nimlich dann, wenn das Problem Symmetrieeigenschaften aufweist, 
wie dies am einfachsten Beispiel der Zylindersymmetrie in I gezeigt wurde. 


& = C,exp | 


(Hingegangen am 6. Juni 1957) 


Zur theoretischen Bestimmung des Dellendruckes 
bei abgeléster Strémung’ 


Helmut H. Korst, Urbana, Ill., U. S. A.2, dzt. Wien? 
Mit 9 Textabbildungen 


Zusammenfassung. Hin einfaches Strémungsmodell, das zum Studium abgeléster Strémungen 
entwickelt wurde, erweist sich als offenbar geeignet, bei Auftreten von hinlanglich ebenen 
Mischungszonen die wesentlichsten EinfluBgr6Ben, hier Machsche Zahl und — in einer Grenz- 
betrachtung — die Reynoldssche Zahl, zu erfassen. Asymptotisches Verhalten bei dinnen Grenz- 
schichtdicken vor der Ablésungsstelle wird aufgezeigt und als den niedrigsten Totwasserdriicken 
zugeordnet gefunden. Diese Grenzdriicke kénnen in Fallen, wo ,,geschlossene Totwassergebiete“ 
auftreten, ohne empirische Information gerechnet werden, und zwar sowohl fur laminare als auch 
fur turbulente Strahlmischung. Bei ,,offenem Totwasser“‘ ist dazu eine einzige VersuchsgréBe, 
iiber die einige Kenntnis vorliegt, erforderlich. 


Bezeichnungen: 


G, pro Breiteneinheit und Zeiteinheit in das Totwasser langsam einge- 
blasener Massenstrom. 

G, pro Breiten- und Zeiteinheit zwischen der Ablésungsstromlinie und 
der Stromlinie s gefiihrter Massenstrom. 


1 Die Forschungsarbeiten des Verfassers wurden unter teilweiser Unterstiitzung durch die 
U. 8. Airforce OSR, Fluid Dynamics Division, unter Contract AF 18 (600)-392, ,,Transonic 
and Supersonic Flow of a Real Fluid at Abrupt Increase in Cross Section“, durchgefuhrt. 

2 Professor of Mechanical Engineering University of Illinois. 

3 Als Senior Postdoctoral Fellow der National Science Foundation, Washington D.C. 
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absoluter Druck. 

absolute Temperatur. 
Geschwindigkeitskomponente in der x-Richtung. 
natiirliche Koordinaten. 

Bezugskoordinaten. 

Integrationsvariable. 

Grenzschichtdicke. 

(scheinbare turbulente) kinematische Zihigkeit. 
Stromlinienneigung. 

spezifische Masse. 


lane 
DM DBemenNree RN 


Indizes kennzeichnen: 


+ = 1, 2, 3,4, Strdémungsquerschnitte (s. Abb. 1). 
a die dem Scherstr6mungsbereich unmittelbar benachbarte freie AuBen- 
stromung. 
die diskriminierende Stromlinie. 
die Ablosungsstromlinie. 
Staubedingung. 
groBen Wert von y, wo Scherstroémung und freie Stromung sich nur 
mehr unwesentlich unterscheiden. 
co ungestérte Anstrémungsverhiltnisse. 
red reduzierte Anstr6mungsverhiltnisse. 
s Stromlinie, durch Gl. (16) definiert. 


Peete) aes foie 


Dimensionslose GroBen: 
1 
O'My Cina = (1 on ae) * Croccosche Zahl. 
M Machsche Zahl. 
P,/P, Totwasserdruckverhaltnis (base pressure ratio). 
Re Reynoldssche Zahl. 
¢, = y/6, dimensionslose y-Koordinate. 
€ = y/6, dimensionslose y-Koordinate. 
n =(C-np dimensionslose Koordinate. 
np = 1/(2 Vé) Positionsparameter. 
€ transformierte x-Koordinate. 
gy dimensionslose Geschwindigkeit, bezogen auf die freie AuBenstromung 
im gleichen Stro6mungsquerschnitt durch y; = u,/U;q. 
y = «x/6, dimensionslose x-Koordinate. 
o Ahnlichkeitsparameter fiir die homogene Koordinate y/s. 
k =c,/c, Quotient der spezifischen Warmen. 


I 


I. Einleitung 


Die klassische Behandlung des Ablésungsproblems beruht auf der Betrachtung 
des Grenzschicht-Strémungsprofiles, das man unter Beniitzung des von der AuBen- 
strémung aufgeprigten Druckverlaufes mit Hilfe der Grenzschichtgleichungen strom- 
abwarts verfolgt. Obwohl der Charakter der Grenzschichtgleichungen parabolisch 
ist, wird durch den Druckverlauf, der wahlweise entweder von der Potentialtheorie 
oder dem Experiment entnommen wird, dem Einflu8 der Gesamtstrémung Rechnung 


getragen. Die Ergebnisse solcher Rechnungen erweisen sich jedoch vielfach als un- 
befriedigend. 
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Demgegeniiber ist in neuerer Zeit angestrebt worden, den Ablosungsvorgang aus 
der Wechselwirkung zwischen der nahezu schubspannungsfreien AuBenstromung und 
der Strahlmischungszone entlang des Ablésungsgebietes unter besonderer Beriicksichti- 
gung der erforderlichen Drucksteigerung am Ende der Delle zu erklaren. 

Crocco und Lees* haben bereits ein derartiges Konzept beniitzt und eine 
Methode zur umfassenden Behandlung solcher Probleme entwickelt, die qualitativ 
befriedigende Aufschliisse fiir eine Reihe bedeutsamer Falle im Uberschallbereich 
liefert. Von gréBerer physikalischer Anschaulichkeit erwies sich jedoch ein Strémungs- 
modell, das D. R. Chapman zunichst zur halbempirischen Auswertung von Dellen- 
druck- (base pressure) Problemen vorschlug®. Durch Hinzufiigung von als wesentlich 
erkannten Verfeinerungen, die den Strahlmischungsvorgang betreffen® und durch 
eine energiemafig formulierte Bedingung fiir die Drucksteigerung am Ende der 
Mischungszone, wodurch die Masse im Totwasser konstant gehalten wird, wurde 
dieses Modell vom Autor zur quantitativen Behandlung von abgelésten Stromungen 
brauchbar gemacht’. * und speziell auf Fille turbulenter Strahlmischung ange- 
_wendet®: 1°. Unabhangig davon ist D. R. Chapman die quantitative Auswertung 
seines Modells fiir laminare Strahlmischung gelungen!!: 22. 

Im folgenden sollen die der Ablésungstheorie zugrunde liegenden Gedanken er- 
értert werden, und zwar zunachst im Sinne des vom Autor vorgeschlagenen Modelles 
fiir Uberschallstrémung®. 


Il. Modell des Ablésungsmechanismus 


1. Das Ablosungsmodell hat vier Str6mungsgebiete (s. Abb. 1): 

a) Die Stromung stromaufwarts der Ablosung, oberhalb des Querschnittes 1. 

b) Die Str6émungsablenkung an der Ablosungsstelle, zwischen 1 und 2. 

c) Die Stromung entlang des Totwassers, mit besonderer Beriicksichtigung des 
Strahlmischungsvorganges, zwischen 2 und 3. 

d) Das Gebiet der Drucksteigerung am Ende des Totwassers, zwischen 3 und 4. 

2. Der Stromungsmechanismus wird durch fiinf Bedingungen formuliert: 

a) Statische Driicke im Bereiche der Scherstr6mung werden durch die benachbarte, 
nahezu gleichformige, schubspannungsfreie AuBenstro6mung eingepragt. 


P,/P;, = 1. (1) 

4 L. Crocco und L. Lees: A Mixing Theory for the Interaction Between Dissipative Flows 
and Nearly Isentropic Streams. Journal of the Aeronautical Sciences 19, 649—676 (1952). 

5 D. R. Chapman: An Analysis of Base Pressures at Supersonic Velocities and Comparison 
With Experiment. NACA TN 2137 (1950), auch NACA Report 1051 (1951). 

6 H. Korst: Auflésung eines ebenen Freistrahlrandes bei Berticksichtigung der urspriing- 
lichen Grenzschichtstr6mung. Osterr. Ingenieur-Archiv 8, 152—157 (1954). 

7 H.H. Korst: Transonic and Supersonic Flow of a Real Fluid at Abrupt Increase in Cross 
Section. Contract No. U. 8. Trust AF 18 (600)-392, Quarterly Progress Report A-PR-4, July 31, 
1953 (ASTIA). 

8 H. H. Korst: Comments on the Effect of Boundary Layer on Sonic Flow Through an 
Abrupt Cross-Sectional Area Change. Journal of the Aeronautical Sciences 21, 568—569 (1954). 

9 H. H. Korst, R. H. Page and M. E. Childs: A Theory for Base Pressures in 'Transonic 
and Supersonic Flow. University of Illinois, ME-TN-392-2, OSR-TN-55-89, Contract No. AF 
18 (600)-392, March, 1955. 

10 H. H. Korst: A Theory for Base Pressures in Transonic and Supersonic Flow. Journal 
of Applied Mechanics 23, 593—600 (1956). Ag : 

1D, R. Chapman, D. M. Kuehn and H. K. Larson: Preliminary Report on a Study of 
Separated Flows in Supersonic and Subsonic Streams. NACA RM A 55 L 14 (1956). 

12 PD. R. Chapman, D. M. Kuehn and H. K. Larson: Analysis and Experiments on Separated 
Flows in Supersonic and Subsonic Streams. Paper presented at the IX International Congress 
for Applied Mechanics, Sept. 5—13, 1956, Briissel, Belgien. 
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b) Die Ablenkung der freien Uberschallstrémung an der Ablésungsstelle zwischen 1 
und 2 ist durch eine eindeutige Zuordnung von Druckanderung und Richtungséinderung 
gekennzeichnet. Dies sei in funktioneller Schreibweise dargestellt, wobei das Ver- 
haltnis der spezifischen Warmen als Parameter nicht ausdriicklich mitgefiithrt wird: 


P,/P, = P,/P, (M, 0. = 0"): (2) 
Es ist dann auch.die Anderung der Machschen Zahl der freien Stromung durch die 


Richtungsiinderung bestimmt. (Zum Beispiel ware in Abb. 1 ein solcher Zusammenhang 
durch die Prandtl-Meyer-Expansion gegeben.) 


M,, = M, (MM, 02, ca 0, ,)- (3) 


a 


2a 
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Abb. 1. Strémungsmodell fiir ebene abgeléste Uberschallstrémung 


Fiir den Bereich der Scherstromung kann man formal schreiben: 


P2 (CS) = Po [Pr (C1); MM, P,/P,] (4) 
und 


02/0, = 6/0; [1 (C1), M, P,/P,). (5) 


c) Die Strahlmischung entlang des Totwassers wird als ebenes Problem betrachtet 
und erfolgt bei konstantem Druck, eine Forderung, die annaihernd nur dann erfiillt 
sein wird, wenn die Dicke der Mischungszone geniigend klein bleibt gegeniiber allen 
wesentlichen Lingen, die mit dem Totwassergebiet assoziiert werden kénnen. Unter 
diesen Voraussetzungen kann fiir turbulente zweidimensionale kompressible Strahl- 
mischung die vom. Autor entwickelte Theorie’ verwendet werden. Es handelt sich 
dabei um eine Impulsintegralmethode, die auf solche Geschwindigkeitsprofile ange- 
wendet wird, die tatsichliche Anfangs- und Randbedingungen erfiillen und aus der 
Integration einer weitgehend vereinfachten Bewegungsgleichung 

ou ee (20) 
ibe Up, Oy? (6) 


gewonnen werden, in welcher die scheinbare turbulente Zaihigkeit durch eine sinn- 
gemaBe Erweiterung des Prandtlschen Ansatzes erfaBt wird: 
‘Le 
aire) U2, Y Ty): (7) 


Die zur Erfiillung des Impulsintegrals erforderliche Beweglichkeit der Losung 
wird durch ihre Deutung in einem natiirlichen Koordinatensystem (,,intrinsic system 
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of coordinates‘) (x, y) erzielt, welches gegen das Bezugskoordinatensystem (,,reference 


system of coordinates) (X, Y) quer zur Stromungsrichtung verschoben, aber nur 
unwesentlich verdreht ist: 


Xe; | 
Y= Y — Ym (x), wobei | (8) 
Wg. (Oje= 0: 


Das Bezugskoordinatensystem folgt der durch die Potentiallésung bestimmten, den 
tatsichlichen Druckverhaltnissen entsprechenden, freien Stromlinie (,,boundary of 
the corresponding inviscid jet“). Eine Beschrankung auf ahnliche Mischungsprofile 
ist nicht vorausgesetzt, wodurch der Einflu8 der Reynoldsschen Kennzahl nicht 
ausgeschlossen wird. Zufolge des turbulenten Mischungsmechanismus tritt eine 
derzeit nur empirisch erfaBbare Zuordnung von Strémungsprofilen zu Lingen- 
koordinaten auf, da eine Integraltransformation, welche auf die physikalische Liingen- 
koordinate angewendet wird, 

wy 

f= 5 | ily) dy, (9) 


0 


in ihrem Kern den Schubspannungsansatz enthalt. In dimensionsloser Darstellung 
formuliert man dann das Mischungsproblem folgenderma8en: 


O~ << op ) 
GE GES 
p(0, Cy) == 0 fiir 7 32 aa <0; \ (10) 
90,0) =9,(6) fir 0 <o<1, | 
p(0,¢)=1 fir 1<f€<+00, | 
und erhalt als allgemeine Losung 
n 
ee I u Hp _p2 
? = 7 [2 (5), NP; n] = [1 + ert (n — np)] 4 vies lel )e dp, (11) 
n—np 
worin der Positionsparameter 
np = 1/(2 Vé) (12) 
und 
7 = C Np. (13) 


Von groBer Bedeutung ist die Betrachtung der in groBeren Entfernungen strom- 
abwiarts asymptotisch erreichten Losung, 7p —> 0, welche das voll ausgebildete, turbu- 
lente Strémungsprofil in der Mischungszone darstellt. Die einer quantitativen Aus- 
wertung schwer zugiinglichen, durch Gl. (4) und Gl. (5) ausgedriickten Beziehungen 
werden in diesem Falle gegenstandslos. AuBerdem wird aber auch die in Gl. (9) ent- 
haltene Beziehung nicht mehr bendtigt. | 

Physikalisch entspricht dies dem Fall einer diinnen Grenzschichtdicke vor der 
Ablésung und kann als Grenzfall fiir hohe Werte der Reynoldsschen Kennzahl, be- 
zogen auf den umstrémten Korper, gedeutet werden. Fiir diese asymptotische Losung 
erhalt man 


p=5(l tert 7), (14) 
wobei 1 = o(y/x). Sodann ist fiir isoenergetische Mischung, wo 


T/T, =1 (15) 
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nach Anwendung des Impulssatzes fiir kompressible Stroémung zufolge der Konti- 
nuitatsgleichung, mit 
Ys 


Ns 
LO Uy [1—C, ?] y dn 
G,=\o-w-dy = = aren a \ 2g? (16) 
Yj ) “ 
nR c 1R 
. Pp ay ey oG, ae tices <e dyn 
1—O,7¢ 1 20, um, I—C,?) =| 1—O, 2g?” 0 
Ng —co 


worin fiir y = y(n) nach Gl. (14) einzusetzen ist. Fiir die Ablésungsstromlinie (jet 
boundary streamline), fiir welche G, = 0 ist, kann man schreiben: 

1R QR 

(  Pdy gp? dn 

| ea |i (18) 


"5 CS 


wonach sich 4; und daher nach Gl. (14) auch g; als von C,, abhingig bestimmen 


lassen. Die Ergebnisse der am ILLIAC™ durchgefiihrten Berechnungen"™ sind in 
Abb. 2 graphisch dargestellt, insbesondere das Hilfsintegral: 


L, (p(n), C, 7] faa 


Sa 
Ss 
& 


Die Verwendung der asymptotischen Losung fiir die Mischungskomponente des 
Ablésungsmodelles wurde vom Autor und seinen Mitarbeitern als _,,eingeschrankte 
Theorie“ (,,restricted theory‘) bezeichnet (yp = 0). Zunachst ist es bemerkenswert, 
daB die mit Hilfe der eingeschrankten Theorie errechneten Totwasserdriicke eine 
untere Grenze fiir den allgemeineren’ Fall darstellen. Von gro8tem Interesse ist aber 
die Tatsache, daf im Rahmen der ,,eingeschrankten Theorie“ analytische Lésungen 
fiir das Ablésungsproblem erhalten werden, die im Falle von ,,geschlossenen Totwasser- 
gebieten* keinerlei empirische Konstanten enthalten, so dafi durch Vergleich zwischen 
Rechnung und Versuch eine scharfe Bewertung der Theorie méglich wird. Treten 
,offene Totwassergebiete“ auf, das heiBt liegen soleche Falle vor, in welchem Masse 
in das Ablésungsgebiet eingeblasen oder vom Ablésungsgebiet abgesaugt wird 
(,,bleeding“‘), so wird die theoretische Losung von dem empirischen Faktor o abhangig, 
auf deren Bedeutung als kennzeichnend fiir die Zunahme der Mischungszonenbreite 
schon vorher in der Literatur hingewiesen wurde. Es steht fest, daB o von der Mach- 
schen Zahl der freien Stromung abhiangt, doch ist das derzeit vorhandene Versuchs- 
material noch zu diirftig, um den Funktionsverlauf genau festzulegen. Es scheint 
jedoch, da8 eine von Tripp neuerdings vorgeschlagene Form,-o = 12 + 2578 M, 
befriedigende Resultate fiir isoenergetische Strahlrander liefert. 


d) Die Drucksteigerung am Ende des Totwassers ist — in Ubereinstimmung mit 2, a 
— von der benachbarten freien Stromung aufgedriickt, wobei wieder eine eindeutige 


13 Blectronic Digital Computer, Engineering Research Laboratory, University of Illinois, 
Wi Sale ci 

44H. H. Korst, R. H. Page and M. E. Childs: Compressible Two-Dimensional Jet Mixing 
at Constant Pressure-Tables of Auxiliary Functions for Fully Developed Mixing Profiles. Uni- 
versity of Illinois ME-TN-392-3, OSR-TN-55-99, Contract No. AF 18 (600)-392, April, 1955. 

16 Wilson Tripp: Analytical and Experimental Investigation of the Base Pressure Behind 
a Blunt Trailing Edge for Supersonic Two-Dimensional Flow — Approaching Streams Having 
the Same Stagation Temperature but Dissimilar Mach Numbers and Stagnation Pressures. 
PhD Thesis, University of Illinois, 1956. 
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Zuordnung zum Ablenkungswinkel vorausgesetzt wird (z. B. durch einen schiefen 


StoB) 
P,[P; = P,/P3 (M2, 6; eakeeie (19) 


Hat man Falle mit geneigten Wanden vor der Ablésung (boattailing) und (oder) nach 
Wiederanlegen der Strémung zu behandeln, so erlaubt die Einfiithrung des Konzeptes von 
,reduzierten Mach-Zahlen“ und ,,reduzierten Dellendriicken“ (,,reduced Mach number 
concept‘) eine vereinheitlichte Darstellung von 
theoretischen und _ versuchsmifigen Ergeb- 
nissen!®, gs. auch Abb. 3. 
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Abb. 3. Ablésung ebener Uberschallstr6mung von einer Wandstufe. a) Geometrische Beziehungen 
fiir die AuSenstro6mung. b) Theoretische Lésung fiir diinne Grenzschichtdicke vor der Ablésung, 
verglichen mit Egginks experimentellen Ergebnissen (ref. 11 in 1°) in ,,reduzierter“ Darstellung 


e) Fir das Totwasser muB die Kontinuitatsbedingung erfillt sein: Zunachst 
definiert man die ,,diskriminierende Stromlinie“’ mit der Koordinate y, (bzw. 7.) 
und der Geschwindigkeit yz derart, daf sie am Ende des Totwassers (Querschnitt 3) 
gerade geniigend mechanische Energie (ausgedriickt durch den Ruhedruck) besitzt, 
um durch isentropische Umwandlung der kinetischen Energie den statischen Druck 
nach der Drucksteigerung am Ende des Totwassergebietes (Querschnitt 4) zu er- 
reichen (,,escape concept‘). 

Po,4/Ps = P,/Ps, (20) 


Poml lam Us> VaiGs aan ae. (21) 


Fiir ein ,,offenes Totwasser“ ist es erforderlich, daB ein von auBen zugefiihrter Massen- 
strom G, am Ende des Totwassergebietes zwischen der Ablésungsstromlinie und der 
diskriminierenden Stromlinie abgefiihrt wird, 


G+ Gz. = 0, (22) 
wobei 
Ti 
Ga=\e-u-dy 
Yj 
ist. 


Bei _ ,,geschlossenem Totwasser“ fallen diese beiden Stromlinien zusammen, 
Ya = Yi- 

3. Die gleichzeitige Erfiillung der durch die oben angefiihrten Gleichungen aus- 
gedriickten Beziehungen ergibt eine eindeutige Lésung fiir das Dellendruckproblem, 
wobei die Stabilitaét der Losung an Hand des monotonen Charakters der in Abb. 7 
gezeigten Kurven erlaiutert werden kann®. Im speziellen Fall, wo die K6rpergeometrie 
und die Anstrémverhaltnisse gegeben sein sollen, mu8B diese Lésung durch Probieren 
mit verschiedenen Werten von P,/P, gefunden werden. Fiir allgemeinere Betrach- 


~ 
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tungen, in welchen der EinfluB der Machschen Zahl untersucht werden soll, wird mit 
Vorteil M 2, (bzw. C,,) als unabhingig Veranderliche fiir die Berechnung gewahlt, 


wodurch eine direkte Auflésung des Gleichungssystems erméglicht wird. 


Il. Anwendungsbeispiele 


1. Eingeschrankte Theorie fiir diinne Grenzschichtdicke vor der Ablésung. 
a) Geschlossenes Totwasser. Ubereinstimmung zwischen theoretischen und ge- 
messenen Dellendriicken ist bemerkenswert gut im Hinblick auf die Tatsache, daB 
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Abb. 4. Vergleich von gemessenen Totwasserdrucken mit der vorliegenden Theorie. —-— Hin- 

geschrankte Theorie (7p = 0); schraffiert: Bereiche von gemessenen Mittelwerten nach K. L. Goin 

(ref. 12 in ?°) und D. R. Chapman et al. (ref. 10 in 1°); — im freien Flug gemessene Werte (ref. 3 
in Goins Arbeit) 


die analytische Lésung keinerlei empirische Werte enthalt. Dies ist durch die Abb. 3 
illustriert. Der asymptotische Charakter der ,,eingeschrankten Losung“ tritt in Abb. 4 
zunichst nur fiir soleche Mach-Zahlen klar hervor, die nicht schallnahe sind. Allerdings 
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Abb. 5. Totwasserdruck bei einer ebenen, plotzlichen Querschnittserweiterung nach einer konver- 
genten Diise. Flachenverhaltnis 2°275. — Eindimensionale Theorie; - - theoretischer Totwasserdruck 
fiir das ebene Problem; —:— theoretischer Totwasserdruck fiir axiale Symmetrie; © Versuchswerte 
fiir axiale Symmetrie®; /\ interpoliert zwischen Wicks Versuchswerten fiir vier Flachen- 


verhdltnisse (4 = 1°76, 2°36, 2°95 und 3°54) 
1 
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liegt die Vermutung nahe, daB die Art der Versuchsdurchfiihrung zu unerfaSten 
Interferenzerscheinungen mit der Dellenform fiihrte. Die Theorie gibt naémlich aus- 
gezeichnete Ergebnisse in Fallen, wo die Beeinflussung der Dellenform durch auBere 
Stérungen theoretisch beriicksichtigt werden kann, s. Abb. 5 und 6. 
b) Offenes Totwasser. Die Wirkung von Masseneinblasung in das Totwassergebiet 
zeigt Abb. 7 in dimensionsloser Auftragung. Wiahrend die dort gezeigten theoretischen 
Kurven auf dem fiir niedrige Mach-Zahlen 
gut bestatigten Wert von o = 12 beruhen, 
hat Wu? spater gezeigt, daB die Uber- 
einstimmung von Versuch und Rechnung 
durch die von Tripp vorgeschlagene Be- 
ziehung fiir o noch verbessert wird. Ab- 
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Abb. 6. Totwasserdruck bei rotationssymmetri- Abb. 7. EinfluB von langsamer Massenein- 


scher plétzlicher Querschnittserweiterung nach 
konvergenter Diise. — Theorie fiir M 1, — | und 


ip = 9; O Versuchswerte von Cole und Mills 
(ref. 36 in °); A eigene Versuchswerte® 


blasung auf den Totwasserdruck bei Ablésung 
von Uberschallstr6mung von einer Wandstufe. 
a) Strémungsbild der AuBenstrémung. b) Theo- 
retische Ergebnisse unter Zugrundelegung von 


o = 12 und Vergleich mit Versuchsergebnissen 
bei einer Mach-Zahl von 2 (Definition der 
lésungsdellen, welche von freien Stro6mun- dimensionslosen ,,bleed number‘, s. 1°) 
gen mit verschiedenen Mach-Zahlen und 
Ruhedriicken, jedoch gleichen Ruhetemperaturen begrenzt werden, fallen ebenfalls in 
diese Gruppe und lassen sich befriedigend durch die Theorie erklaren: 1’, s. Abb. 8 und 9. 
Das Einblasen von hei®er Luft in das Totwassergebiet erfordert die Erweiterung des 


Stro6mungsmodelles im Sinne einer ,,nicht isoenergetischen, kompressiblen“ Mischungs- 


16 Chen-Yuan Wu: The Influence of Finite Bleed Velocities on the Effectiveness of Base 
Bleed in the Two-Dimensional Supersonic Base Pressure Problem. PhD Thesis, University of 
Ilhnois, February, 1957. 

1 H. H. Korst and Wilson Tripp: The Pressure on a Blunt Trailing Edge Separating Two 
Supersonic Two-Dimensional Air Streams of Different Mach Numbers and Stagnation Pressures, 
but Identical Stagnation Temperatures. Paper presented at the Midwest Conference on Solid 
and Fluid Mechanics, University of Michigan, Ann Arbor, Michigan, April 1—2, 1957. 
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theorie fiir turbulente Strémung, die gleichfalls befriedigende Ergebnisse brachte!8. 19, 
Gestiitzt auf die so gewonnenen Erkenntnisse ist jetzt auch die Vorausberechnung 
der durch die AuBenstrémung unter Uber- 
schall-Flugbedingungen bedingten  Kin- 
fliisse auf die in Héhenpriifstanden er- 
mittelten Eigenschaften von Strahltriebs- 
einheiten”® theoretisch méglich geworden. 


| 
GS bie 20 40 60 CO 700 
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Abb. 8. Einflu8 des Ruhedruckverhaltnisses auf Abb. 9. Theoretische Bestimmung der Stré- 
den Totwasserdruck bei einer von zwei Uber- mungsablésung von einem stark angestellten 
schallstromungen gebildeten Ablésungsdelle Uberschallprofil, verglichen mit Versuchswerten 
(hier M,, = M,/ = 1:98). a) Stromungsbild von A. Ferrinach NACA TM 946 (1940) (Modell 
der Au&enstrémung. b) — Theorie fiir yp = 0; GU 3, Moo = 213, «x = 22°). a) Stromungsbild 
© Versuchswerte von W. Tripp” der AuBenstromung. b) --- Theorie, Charak- 
teristikenmethode; — theoretischer Totwasser- 


‘ : ; : POs Vv h t 
2. VerhaltnismaBig wenig theoretische ee Ot: <y ereuchew enue 


Anhaltspunkte sind gegenwiartig fiir solche : 

Falle verfiigbar, wo die Voraussetzungen der ,,eingeschrankten Theorie“ nicht zu- 
treffen. Allerdings ist eine qualitative Erorterung des stabilisierenden EHinflusses der 
Grenzschicht vor der Ablésungsstelle méglich, wodurch diese sich auch im Falle von 
glatten, konvex gekriimmten Wanden als festgelegt erweist®. Die Untersuchungen 


18 R. H. Page: The Non-Isoenergetic Turbulent Jet Mxing of a Two-Dimensional Supersonic 
Jet With Consideration of Its Influence on the Base Pressure Problem. PhD Thesis, University 
of Illinois, February, 1955. 


19 R. H. Page and H. H. Korst: Non-Isoenergetic. Turbulent Compressible Jet Mixing 
With Consideration of Its Influence on the Base Pressure Problem. Transactions of the Forth 
Midwestern Conference on Fluid Mechanics, Purdue University, September, 1955. 


20 EK. M. Cortright Jr.: Some Aerodynamic Considerations of Nozzle-Afterbody Combinations. 
Paper presented at the 24th Annual Meeting of the Institute of the Aeronautical Sciences, Jan. 
23—26, 1956. Preprint No. 614. 
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von Wu! 2! iiber den Einflu8 von endlichen Einblasegeschwindigkeiten in Tot- 
wassergebiete ergaben zwar einen tberraschend engen Bereich fiir eine empirische 
Abklingfunktion von urspriinglich gestorten Mischungszonen, doch miissen erst weitere 
systematische Versuche durchgefiihrt werden, um feststellen zu konnen, ob brauchbare 
Niherungslésungen durch eine beschrankte Verallgemeinerung dieser Funktion 
erwartet werden k6nnen. 

3. Ersetzt man: die Bedingungen II, 2, b und II, 2, d durch die Annahme, dai 
am Ende der Delle der statische Druck der (gegebenenfalls ,,reduzierten“) Haupt- 
stromung erreicht werden muB, so kann fiir Fille mit geschlossenem Totwasser eine 
besonders einfache Lésung fiir den Totwasserdruck gegeben werden. Diese An- 
nahme gestattet aber auch sofort, die Berechnung von Totwasserdriicken auf den 
Unterschallbereich auszudehnen. D. R. Chapman und seine Mitarbeiter haben mit 
diesem Ansatz befriedigende Resultate in Fallen mit laminaren Mischungszonen 


erzielt1!, 12, 


21 C.-Y. Wu and H. H. Korst: The Influence of Finite Bleed Velocities on the Effectiveness 
of Base Bleed in the Two-Dimensional Supersonic Base Pressure Problem. Paper presented at 
the Midwest Conference’ on Solid and Fluid Mechanics, University of Michigan, Ann Arbor, 


Michigan, April 1—2, 1957. 
(Hingegangen am 18. Juni 1957) 


Zur Biegungsbeanspruchung eines Seiles 
auf einer elastisch gefiitterten Rolle 


Von E. Engel, Wien 


Mitteilung aus dem Institut fiir EKisenbahnbau an der Techn. Hochschule Wien 


Mit 3 Textabbildungen 


Zusammenfassung. Es wurde eine Naherungsberechnung ftir den kleinsten Kriimmungs- 
halbmesser eines Seiles auf einer elastisch gefiitterten Rolle gewonnen, die vom theoretischen 
Standpunkt als grob betrachtet werden mu8. Auf die genaue Ermittlung der Druckverteilung 
zwischen Seil und-Rolle, die eine co-fach statisch unbestimmte Aufgabe darstellt, wurde dabei 
verzichtet, und statt dessen wurde ein parabolischer Ansatz fur jene verwendet. Fur die Be- 
rechnung der Formanderungen der Seilkurve wurde von zusatzlichen Naherungen Gebrauch 
gemacht. 

Wie die Zahlenbeispiele zeigen, werden die Ergebnisse um so genauer, je kleiner der Seil- 
ablenkwinkel, je gréBer die Biegesteifigkeit des Seiles und je harter das Rollenfutter ist. Doch 
in Anbetracht der grofen Unsicherheit bei der Bestimmung der Biegesteifigkeit eines Seiles mit 
Beriicksichtigung der inneren Reibung und auch im Hinblick auf die ubrigen, stark idealisierten 
Voraussetzungen, die der vorliegenden Berechnung und auch der strengen Theorie von Czitary 
zugrunde liegen, kénnen die Ergebnisse fiir praktische Zwecke als hinreichend brauchbar betrachtet 
werden. Ferner erkennt man, dai eine VergréBerung des Rollendurchmessers die gréBten Biege- 
beanspruchungen des Seiles auf der Rolle verhaltnismaéBig wenig vermindert; der gleiche Gewinn 
kénnte auch durch ein entsprechend weiches und dickes Rollenfutter erzielt werden. 


I. Aufgabenstellung 


Die aus der Erfahrung bekannte Tatsache, da8 sich eine Fiitterung der Seiltrag- 
oder Ablenkrollen mit einem weichen, elastischen Werkstoff auf die Verwendungs- 
dauer der Drahtseile sehr giinstig auswirkt, wurde theoretisch von E. Czitary! unter- 
sucht. Es wurde dabei eine genaue Berechnung der Kriimmungsverhiltnisse und 
der spezifischen Pressungen im Beriihrungsbereich zwischen Seil und Rolle durch- 
gefiihrt. Es zeigte sich, da der kleinste Kriimmungsradius auf einer elastisch ge- 


1 E. Czitary: Verhalten eines Drahtseiles auf einer Rolle mit elastischem Futter. Osterr. 
Ing.-Arch. 10, H. 4 (1956). 
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fiitterten Rolle bei geringer Seilablenkung wesentlich gréBer ist als auf einer nicht 
verformbaren Rolle. 

Die praktische Ermittlung der Biegespannungen bei Zugseilen von Seilschwebe- 
bahnen, Standseilbahnen, Sesselliften usw. wird bisher. stets so durchgefiihrt, daB als 
mafgebender kleinster Kriimmungshalbmesser des Seiles der Halbmesser der Antrieb- 
scheibe, Seiltrommel usw. in Rechnung gestellt wird. Es besteht aber kein Zweifel, 
daB in manchen Fallen auf den Ablenkrollen des Zugseiles trotz der begrenzten Seil- 
ablenkung von héchstens 8 bis 10% viel kleinere Kriimmungsradien und damit 
groBere Biegungsspannungen auftreten kénnen. Ebenso ist der Einflu8 des Rollen- 
durchmessers auf die Biegebeanspruchung des Seiles noch nicht endgiiltig geklart. 
Die strenge Theorie von Czitary lat eine geschlossene Darstellung der Ergebnisse, 
insbesondere des kleinsten Kriimmungshalbmessers des Seiles nicht zu und erfordert 
liberdies bei der Auswertung eine erhebliche Rechenarbeit. Es erscheint daher 
winschenswert, eine Niherungsberechnung zu versuchen, mit dem Ziel, einen 
einfachen Ausdruck fiir den kleinsten Kriimmungshalbmesser eines Seiles auf einer 
elastisch gefiitterten Rolle zu gewinnen, der noch hinreichend brauchbare Ergebnisse 
liefert?. 

Als Vergleichsgrundlage fiir die Beurteilung des Fehlers, den die folgende 
Naherungsberechnung ergibt, soll die strenge Theorie von Czitary dienen; es werden 
daher dieselben Voraussetzungen: gemacht, die jener zugrunde liegen: 

1. Der Seilablenkwinkel 2 « an der Rolle sei klein, so daB tg « <1 und cos «= 1 
wird. 

2. Die Seiloberfliche sei vollkommen glatt. 

3. Die Biegesteifigkeit B des Seiles (analog zur Biegesteifigkeit eines geraden, 
homogenen Stabes, die das Produkt aus Elastizitaétsmodul und Traigheitsmoment 
des Stabquerschnittes bedeutet) sei konstant, das heiSt unabhangig von der 
Kriimmung. 

4. Die Rolle mit dem Halbmesser R wird im Beriihrungsbereich mit dem Seil 
durch eine quadratische Parabel ersetzt; ihre Gleichung (s. Abb. 2) lautet: 

2 
Ye = (1) 

5. Die Zusammendriickung des Rollenfutters 7 sei propor- 

tional dem spezifischen AnpreBdruck p je Langeneinheit: 


=e p. (2) 
Fiir die Konstante c erhalt man entsprechend Abb. 1 
Abb. 1. Querschnitt 
ie wh BL x9 (3) durch ein Rollenfutter 
— 2d Hg 


wenn 6 die Dicke des Rollenfutters, Hg den Zusammendriickungsmodul desselben 
und d den Seildurchmesser bezeichnet. (Das Rollenfutter erhalt meist eine Rille, 
die zwei Drittel des Seildurchmessers breit ist.) 


Il. Die Seilkurve bei nichtverformbarer Rolle 


Die geometrischen Verhiiltnisse eines biegesteifen Seiles auf einer nicht verform- 
baren Rolle wurden zuerst von Isaachsen® untersucht und kénnen als hinreichend 
geklirt betrachtet werden, so daB die ausfiihrliche Herleitung der Ergebnisse unter- 
bleiben darf. Der Ablenkwinkel 2 « des Seiles an der Rolle, der gleich ist jenem Winkel, 


2 Die Anregung hiezu verdankt der Verfasser ebenfalls Prof. E. Czitary. 
8 Igaachsen: Die Beanspruchung der Tragseile. Z. Ver. dtsch. Ing. 1907, 655ff. Siehe auch 


FuBnote 1. 
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den die Asymptoten an die Seilkurve miteinander einschlieBen, ist durch den Rollen- 
druck V und durch die Seilspannkraft S bestimmt: 


tga =e. (4) 


Abb. 2. Form des Seiles auf einer starren (strichliert) und elastisch gefiitterten Rolle 


Mit den Bezeichnungen der Abb. 2 erhilt man die Gleichung der Biegelinie des 
Seiles in der Form 


[s. 
ee 5” + 2 tg o, (5) 
wobei der Abstand des Scheitels der Seilkurve vom Asymptotenschnittpunkt 
A* = Oo" te x | (6) 


betrigt. Fiir den Kriimmungsradius des Seiles 9,* an der Berithrungsstelle mit der 
Rolle erhalt man das bekannte Ergebnis 


1 gp £ 
ae =|) Fee. (7) 


Dieser Wert ist vom Rollenhalbmesser R unabhangig. Dies gilt allerdings nur fir 
0o* > FR, solange sich also das Seil an die Rolle nicht anschmiegt. 


ILI. Die gréBte Kriimmung des Seiles auf elastisch gefiitterter Rolle 


Besitzt die Rolle ein elastisch verformbares Futter, so wird der Rollendruck nicht 
punktformig, sondern lings einer Strecke 2] auf das Seil ttbertragen, wobei sich 
das Seil in das Rollenfutter eindriickt. Man kann sich statt dessen ebensogut das 
Seil in jener Lage, die es auf einer starren Rolle einnimmt, festgehalten (wie in Abb. 2 
strichliert dargestellt) und dafir die geftitterte Rolle in das Seil eingedriickt denken. 
AuBerhalb des Beritithrungsbereiches bleibt die Seilkurve auch nach der Verformung 
des Rollenfutters unverandert, wahrend sie sich im Beriihrungsbereich selbst wegen 
der veraénderten Druckverteilung verflachen wird. 

VoraussetzungsgemaB wird die nichtverformte Rolle im Beriihrungsbereich durch 
eine quadratische Parabel ersetzt. Ebenso kann man sich die wirklich auftretende 
Seilkurve in der Umgebung des Scheitels (x = 0) in eine konvergente Potenzreihe 
entwickelt denken, von der auch wieder nur die Glieder bis zur zweiten Potenz be- 
ricksichtigt werden. Dies bedeutet noch keine Verschlechterung der Genauigkeit. 
Die Zusammendriickungen 7 sind die Differenz aus den Ordinaten der Seilkurve und 
denen der nicht verformten Rolle und haben infolgedessen gleichfalls einen parabolischen 
Verlauf langs der Berithrungsstrecke mit dem GréBtwert 7, so dab 


= 1p (1 — +) (8) 
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Daraus folgt dann fiir den Druck p zwischen Seil und Rollenfutter, der voraussetzungs- 
gemaB der Kindriickung 7 proportional ist, ebenso ein parabolischer Verlauf 


P =po(1— Fr] (9) 


mit dem GroBtwert p, im Scheitel. Aus der Gleichgewichtsbedingung fiir die Rolle 
in der y-Richtung 


4 
V=2\pde= Fpl (10) 
0 
ergibt sich 
3 V 
Po = ay 
und mit (4) 
3S tg « 
HL oes oe (10a) 


Die groé&te Eindriickung in « = 0 betrigt daher 


iio — Capa eee TS ae (11) 


Die Berechnung des kleinsten Kriimmungshalbmessers 0, des Seiles im Scheitel 
erfolgt aus der Gleichung der elastischen Linie 


alah Shi 
alice (12) 
Dabei ist M, das Biegemoment im Scheitel. Aus Abb. 2 entnimmt man 
1 
M,=S-h—|pxdz. (13) 
0 
Mit (9) und (10a) liefert die ausgefiihrte Integration 
M, =8(h— S1tg«). (13a) 
Damit erhalt man aus (12) 
B 3 
Sg = hel tga. (14) 


Bei der Eindriickung der Rolle in das Seil hat sich der Abstand h* des Scheitels der 
Seilkurve vom Asymptotenschnittpunkt auf A vergréBert, wihrend der Kriimmungs- 
halbmesser von 0,)* auf 0, angewachsen ist. Fiir die Berechnung von h und / werden 
folgende zwei Naherungen verwendet: 


1. Wegen des linearen Zusammenhanges von h* und 0,* gemaf (6) wird auch fiir 
die Anderung der beiden GréBen der Ansatz versucht 


h — h® = (05 — Qo*) tg? «. (15) 


Diese Beziehung ware nur dann exakt giiltig, wenn nach der Eindriickung die Seil- 
kurve innerhalb und auch auSerhalb des Beriihrungsbereiches durch eine Gleichung 
von der Form (5) und (6), nur mit 0, anstatt 0)*, dargestellt werden kénnte. Dies 
bedeutet aber im Hinblick auf das vorher Gesagte eine Naherung. Mit (6) und (15) 


erhalt man dann 
h = 0, tg? x. (15a) 


Diese Niherung ergibt fiir h stets etwas zu groBe Werte. Bei einer genaueren Erfassung 
miiBte naher auf die Form der Biegelinie des Seiles eingegangen werden. 
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2. Die Ordinaten der Seilkurve y, sind gemiB Abb. 2 die Summe der Ordinaten 
des Kreises und der Eindriickungen: 


Ys=h— Mo + Ye +0: (16) 
aii 6 
Pir ule 
ees 20 
und wegen (1) 
" 1 
Yi: mayer 
und wegen (8) und (11) 
2 3c Stg x 
Tae [2 No aa 13 ? 


1 th 3cS tg « 
Poi R B 
oder 
fe V/ ee (17) 
R Qo 


Die GréBe der Berihrungslinge ist aber fiir den Wert des kleinsten Kriimmungs- 
halbmessers 9, von geringerer Bedeutung*, und man kann in guter Naherung / so 
berechnen, als ware das Seil bei der Eindriickung in die Rolle unverformt geblieben, 
das hei8t man kann in (17) @, durch @,* ersetzen: 


_°/ 3cStg oa 
a os (17a) 
Pe Rae 


Diese Naherung liefert etwas zu grofe Werte fir die Beriihrungslange; sie wird um 
so genauer, je groBer die Biegesteifigkeit des Seiles ist. 


Die beiden Naherungen (15) und (17a), die beide etwas zu groBe Werte liefern, 
gleichen sich aber bei der Differenzbildung in (14) gegenseitig wieder teilweise aus. 
Man erhalt mit ihnen aus (14) nach geringer Umformung 


8 tg x B 
Oo” cw ane Ste? a 0, (18) 


also eine quadratische Bestimmungsgleichung fiir den kleinsten Kriimmungshalb- 
messer 0). Die Auflosung ergibt das angestrebte Ergebnis 


16 tg « 1 Ry eae 53 
Cory dail sesh bos: (18a) 


Dabei ist / aus (17) bzw. (17 a) zu berechnen. 


Fiir eine nicht verformbare Rolle (c = 0) erhalt man aus (18) durch Spezialisierung 
B 
20 reve 0? = S tg? «’ 
und dieses Ergebnis ist identisch mit der Formel von Isaachsen. Es mu8 daher er- 
wartet werden, da die Naherungsformel (18a) um so genauere Werte fiir den 
Krimmungshalbmesser Qo liefern wird, je hirter das Rollenfutter ist. 


* Nach der strengen Rechnung diirfte auch die Voraussetzung 4 nicht mehr verwendet werden, 
denn die Ermittlung der Berithrungslange 2 / beruht auf der genauen Bestimmung eines schleifenden 
Schnittes zweier Kurven! 


R. Bruniak: Uber die Ablésung der Grenzschicht beim VerdichtungsstoB 243 


Zahlenbeispiele 


Ue Es werden zwei Beispiele mit verschiedener Biegesteifigkeit des Seiles gerechnet und die 
Ergebnisse der Naherungsrechnung denen der genauen Theorie von Czitary gegenubergestellt. 


Annahmen: 


Litzenseil der Bauart h + 6 (1+ 6+ 12), 

Drahtdurchmesser dy) = 1:9 mm, 

Seildurchmesser d = 30 mm, 

Anzahl der Drahte 7 = 114, 

Elastizitatsmodul des Drahtmaterials EH = 2200000 
kg/em?, 

Seilspannkraft S = 6000 kg, 

Rollenhalbmesser R = 250 mm, 

Dicke des Rollenfutters 6 = 60 mm, 

Elastizitaétsmodul des Futters Hg = 850 kg/cm?,, 


Ablenkwinkel an der Rolle tg 2 « = 0°08. 


750 


700 


50 


Abb. 3. Kleinster} Kriimmungshalb- 
messer des Seiles in Abhangigkeit 


Damit ergeben sich folgende kleinste Kriimmungs- 
vom Rollendurchmesser 


halbmesser des Seiles an der Rolle: 


Kriimmungshalbmesser in cm 


Biegesteifigkeit in cm* Rolle gefiittert. Fehler % 
Rolle starr |——— = 
| N&herung genaue Werte 
Ss Oise | 
Bit E aa eae 0:0073 41°0 85°0 94°4 10 
B—=10 Be = 0:073 130°0 165°6 169-0 2 


2. Der EimfiuB 


des Rollenhalbmessers auf die Biegebeanspruchung des Seiles laBt sich auf 


Grund der Ergebnisse (17) und (18a) beurteilen. Je groBer die Biegesteifigkeit B des Seiles und 
je harter das Rollenfutter ist, desto unbedeutender wird die GréBe des Rollenhalbmessers R. 
In Abb. 3 ist der kleinste Kriimmungshalbmesser des Seiles fiir die beiden Werte der Biege- 
steifigkeit aus Beispiel 1 graphisch dargestellt. 

(Hingegangen am 26. Juni 1957) 


Uber die Ablésung der Grenzschicht beim VerdichtungsstoB! 
Von R. Bruniak, Wien 


In dem genannten Aufsatz wurde versucht, mit Hilfe des Impulssatzes unter 
gewissen Vernachlissigungen und mit der Annahme von Pohlhausen-Profilen eine 
Beziehung fiir die Ablésung der Grenzschicht beim VerdichtungsstoB8 zu finden. Es 


hat sich folgendes ergeben: 


6 
! (9 uv) udy — J (ou) udy = (pz, — py) 0. 

Das erste Integral ist an der Stelle zu nehmen, wo der VerdichtungsstoB entsteht, 
das zweite knapp nach der Stelle, wo die Geschwindigkeit uw, nach dem Sto8 in der 
freien Stroémung bereits erreicht ist. 6 ist die Grenzschichtdicke und wird in dem 
kurzen Bereich als konstant angenommen. Wenn man nun die Stromdichte 9 u = 0 u 
auch iiber die Grenzschichtdicke als konstant annimmt, wie es geschehen ist, so ist 
dies eine gewisse Naherung. 


1 Yusatz zur gleichnamigen Originalarbeit, Osterr. Ing.-Arch. 10, H 2/3, 129ff. 
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y 


Setzt man 7 = so hat man 


6? 
1 1 
| ow ddn — \ 9 u? ddy = (p2 — Px) 6 
0 0 : 
oder 
r oe 
| uedn — | @wdn = pz — Pr. 
6 6 
Fir die Geschwindigkeiten w und u gelten die Gleichungen 


uU = |(29 — 297 + 7) 4 a (7 Bp +319 — 79)| u, = [F(y) + 4G(q)] 4 


wobei 


Pq) = 2H 27 7, 


) 
G(n) == (n 7 el bcm 


ist. Ebenso 


As 


UG = |(2y — 29> + 8) +E (yn — 89? + 37% — m)| We = [F(q) + An G(m)] We 


Wenn 7, die Stautemperatur ist, so gelte 7’) = konst. in und auBerhalb der 
Grenzschicht. Fiir die turbulente Grenzschicht ist dies richtig und fiir die laminare 
sind die Abweichungen von 7’, = konst. gering?. 

Wird keine Wirme zu- oder abgefiihrt, so gilt fiir einen Stromfaden 


Toes 
CO a ag nl Cy 


wenn c, die Ruheschallgeschwindigkeit und c¢ die 6rtliche Schallgeschwindigkeit be- 
deuten. Somit erhalt man 

On x—1 wv 
Ps ta > 


Nun ist aber ? =xg RT und — = M, daher 


ee: aurea te 
Siciee MS Resa 


wobei M die ortliche Mach-Zahl ist. Da zwischen der Grtlichen und der kritischen 
Mach-Zahl die Beziehung besteht 


rem ia 
iit hie 2° 
5 3 M 
so hat man 
Dot ee as SE: 
jie: ea 9 M? eae 5) ea ephee| a 
2 2 me 
x+] 
as aay 1 x+1 
ppd gol UO re ae eee 
3 3 Vise 
Nun 1&8t sich die Dichte g nach der Gasgleichung bestimmen und man erhiilt 
Pp Pp %*+1 


i 


eC GkT. GRE, (etl) oe 


* A. Busemann: Gasdynamik. Handbuch der Experimentalphysik von Wien-Harms, 
Bd. IV, 1. — E. Pohlhausen: Der Warmeaustausch zwischen festen Kérpern und Fliissig- 
keiten mit kleiner Reibung und kleiner Warmeleitung. Z. angew. Math. Mechan. 1 (1921) 
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Nimmt man an, da der Druck in und auBerhalb der Grenzschicht konstant ist, so 
ist vor dem Sto p =p, und nach dem StoB p = p,. Es ist daher 


0 Pi = x+1 9 Pi x+l1 
Rae a's yee [| 2 

erae eel) ez 1) an Cokes VaR eee 1)( | M® 

d “s 
un 

0 Pe “n+ ] 
== a : 
Gi dis Ah (x + 1) (% 1) = M,*? 
2 


Geht man nun mit diesen Werten von 0 und 0 in die Integrale hinein, so erhalt man 
iL 


Oe ee - (LF +2, GP u2dy — 


Grivel. aS Pe) 
\ a M, D1 + GP 


Ray ey : *[F + 1,GP uP dyn = pe. — py. 
7” [F + AP 


Ferner 


dn 


t—— J MM, (F + 1,47 


1 
[F + AGP 
ay} 3% x + 1 


Aen, G}? 
M,(F + A, GP 


dy = Pes Ps 


i 


eee 


ey 
byl 
Sh, 
o 
> = 
| 


und 


ae M,* [F + 1,GP 


1 
? 2 
p, U2 ae fy + A, G4] lah z 
3 Lp is 
0 


Os Shame 
eo : M,*?[F + As GP 


Pees ’ cot ee — dry = [1 — Bi] 2 Po. 
he 
(0) 


Es ist somit ie 
Pr Uy (x sol es Tees EY ees Ler Cat OS 7 AE | 
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Setzt man 


1 1 
re er 2 ee ee 
so hat man 
(6— M,*) ) M,** ( cs ; a F 6 ( ME 1) MM, 
6M? —1- oe Hf | 1 6.0 81s 
Ge 5 x 6 M,* 
ft) 0 


Zur weiteren Abkiirzung sei noch 


(6 — M,*?) M,*4 6 (M,*2 — 1) M,*? 


1, = 6 M,*? —1 und Z, = 6 M,*? —1 
Es bleibt nun : 
te 
Z, <a ee dn \e ots iN = 22 
; 6 6 M,*? 
0 0 
Der Ausdruck F + AG wird mit A gréBer und somit Lie mit wachsendem 4A 
kleiner. Man setzt nun fiir 2, = + 12 und fiir A, = — 12 und findet ftir die A- bzw. 
B-Werte: 
A=+12 ee A A=—12 n | B 
OOS Tes 0 Co 
Ol 8°964 O'l | 472°6 
O72 | 3°035 0-2 36°3 
0°3 1°807 0:34} 8°86 
0-4 | 1/358 0-4 | 3°67 
05 | 1-181 0:5 | 216 
0-6 | 1:103 06 | 1:53 
0:7 1:050 OF 1:24 
0-8 | 1°016 0°38 | a Oy 
0-9 | 1:016 0°9 1:02 
1 I Lact I! 
Es sei nun M,** — 1°6. Im ersten Integral, das nach der Simpsonschen Regel 
ausgewertet werden soll, sei A = +12. Im zweiten Integral sei 4 = — 12. Wenn 
1 1 
is 1 te 1 
= M*? dy und [y= 1 dn, 
4-4 ) 2— sire 
0 0 


so erhélt man J, = 0°895 und J, = 0°516. Ferner ergibt sich Z, 7, = 1:172 und 
I, +Z, = 1:186, somit 
ie ae no ae 


Nun sei M,** = 15. Dann erhalt man auf die gleiche Weise Z, J, > J, + Z,. Will 
man im ersten Fall die Gleichheit herstellen, so mu8 man J, EN eh ‘Aw atk A. 
Das heiBt auch bei einem Parameterwert, der in der seabilen Grenzschicht nicht 
auftritt, wirde es zu einer Ablésung der Grenzschicht nach dem Sto8 kommen. Im 
zweiten Fall miiBte man J, und somit 4 verkleinern, das hei®t ein Parameterwert 
kleiner als + 12 wiirde zu einer Ablésung fiihren. Aber 4 = + 12 mu noch keine 
Ablosung ergeben. Fir M,*? = 1:5 ist M, = 1-291. 


(Eingegangen am 9. Juli 1957) 
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Buchbesprechungen 


Physikalische Chemie als Einftihrung in die chemische Technik. Von O. Fuchs. (Grundlagen 
der chemischen Technik: Band 2.) Mit 98 Textabb., 4968S. Aarau und Frankfurt/Main: 
Verlag H. R. Sauerlander & Co. 1957. Geb. sfr. 37-—. 


Im Hinblick auf die grundlegende Bedeutung, die der physikalischen Chemie sowohl in 
stofflicher wie auch in methodischer Hinsicht in der chemischen Technik und Verfahrenstechnik 
zukommt, ist ein Werk mit der Zielsetzung des vorliegenden Buches sehr zu begriiBen. Bei 
der erforderlichen Abgrenzung des Stoffes auf die wichtigsten Teilgebiete wurden vor allem 
der Darstellung der chemischen Gleichgewichtslehre und der chemischen Kinetik Rechnung 
getragen, wobei Phasengrenzflachen, Mischphasen, Reaktionen mit festen Phasen sowie elektro- 
chemische Vorgange eingehender besprochen werden. 

Der titberwiegende Teil des Inhaltes ist altbekannt, wenn auch auf gewissen Gebieten neuere 
Arbeiten beriicksichtigt sind. In den ersten Kapiteln (Die chemischen Elemente, Praktische 
Stéchiometrie, Energetische Grundlagen), die zum Teil recht elementar und breit — wenigstens 
fiir die Chemiker — gehalten sind, kommt besonders zum Ausdruck, daB das Buch sich auch an 
den Ingenieur wendet. Ansonsten ist das Werk ahnlich wie ein Repetitorium einer Einfiihrungs- 
vorlesung tiber Physikalische Chemie, besitzt jedoch den erheblichen Vorteil, daB dem Leser 
vor allem die groBen Zusammenhange und grundlegenden GesetzmaGigkeiten nahegebracht 
werden. Einschrankend mu8 man jedoch feststellen, da&8 Vorkenntnisse an Begriffen und 
Definitionen nétig sind, um das Buch mit Vorteil gebrauchen zu kénnen, was sich schon daraus 
ergibt, daB auf die Begrindung von GesetzmaBigkeiten nicht eingegangen wird. 

Fur Studierende macht sich das Fehlen des didaktischen Momentes unangenehm bemerkbar. 
Es ist em Mangel, wenn wichtige Begriffe, mit denen in dem Buch operiert wird, nicht mit den 
Namen bezeichnet werden, unter denen sie allgemein bekannt sind (z. B. Pauli-Prinzip, Halb- 
wertszeit, Kirchhoffscher Satz); oder wenn Begriffe nicht genitigend prazise (etwa Warmekapazitat 
— Molwarme — spezifische Warme; oder in der Kinetik: Unterschied zwischen Molekularitat 
und Ordnung) oder erst an spaterer Stelle definiert werden, nachdem sie zuvor schon mehrfach 
verwendet werden. Die Verwendung des Diminutivzeichens ,,d‘* sowohl fiir exakte Differentiale 
wie auch fur kleine Betrage von Arbeit oder Warme ist nicht zweckmaBig. Auch sind manche 
Formulierungen nicht ganz gegliickt bzw. unrichtig: z. B. 8.59 fiir die Energiedquivalente; 
S. 60 fiir die mit zunehmender Temperatur gréBer werdende Geschwindigkeit thermisch bewegter 
Teilchen; S. 94. ftir die Boyle-Temperatur; S. 351: statt Wasserstoff-Normalelektrode sollte es 
nur Wasserstoffelektrode heiBen; auch die angeschriebene Reaktion macht CO, als Saéure im 
Sinne Broénsteds nicht deutlich. 

Als vorteilhaft darf hingegen hervorgehoben werden, da’ am Ende der einzelnen Kapitel 
Literaturzusammenstellungen angegeben sind, die ein vertiefteres Studieren erleichtern. Bei 
den Zitaten im Text sollten immer die Namen der Autoren angegeben werden. 

Die beispielsweise angefiihrten Mangel, die teils sicher auf das Eimengen eines groBen Stoffes 
auf kleinem Raum und die dadurch nétige Kiirze zurtickzufiihren sind, beeintrachtigen den 
Wert des fliissig geschriebenen Buches kaum. Den gré8ten Nutzen aus dem Werk werden vor 
allem Praktiker, sowohl Chemiker wie Ingenieure, ziehen. 

Druck und Ausstattung sind ausgezeichnet. H. Nowotny-A. Maschka, Wien 


Particulate Clouds: Dusts, Smokes and Mists. Their physics and physical chemistry and industrial 
and environmental aspects. Von H. L. Green und W. R. Lane. Mit 123 Textabb., XIX, 436 8. 
London: E & F. N. Spon Ltd. 1957. £3.10s.0d. : 


Das ziemlich umfangreiche Werk ist eine ausfiihrliche Zusammenstellung der auf verschiedenen 
Gebieten der Technik, Physik und Chemie in den letzten 30 Jahren entstandenen wissenschaft- 
lichen Arbeiten iiber die Staubtechnik, welche die drei Hauptgebiete: Staub (dust), Rauch (smoke) 
und Neben (mist) umfassen. Das Buch ist in zwei Teile gegliedert: 1. Grundlagen der Physik 
und Physikalischen Chemie, 2. Behandlung des Problems vom industriellen Standpunkt. Die 
Grundlagen umfassen das Einzelverhalten von Staubteilchen bzw. einer Wolke in Gasen, 
physikalische und optische Eigenschaften von Staub, Rauch und Nebel, wobei der atmospharische 
und der industriell entstandene getrennt behandelt wird. Ferner Staubabsetzung und Filterung 
(mechanische und elektrische). Im zweiten, praktischen Teil werden. industrielle Erfahrungs- 
werte, MeBverfahren zur Bestimmung der TeilchengréBe besprochen und auch gesundheitliche 
Schidigungen. Nach jedem Abschnitt finden sich ausfithrliche Literaturangaben und am SchluB 
des Buches ein Autoren- und Sachverzeichnis. Dadurch wird das Buch besonders wertvoll fiir 
alle an den in neuerer Zeit immer mehr an Bedeutung gewinnenden Problemen der Staubtechnik 
interessierten Technikern. F. Magyar, Wien 
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Flichentragwerke. Einfiihrung in die Elastostatik der Scheiben, Platten, Schalen und Falt- 
werke. Von K.Girkmann. Vierte, verbesserte und vermehrte Auflage. Mit 316 Textabb., 
XIX, 5968. Wien: Springer-Verlag. 1956. Geb. 8S 396-—, DM 66:—, $ 15°70, sfr. 67°60. 


Das im Inland und Ausland bestens bekannte Werk liegt nunmehr bereits in der vierten Auf- 
lage vor. Gegentiber der vorangegangenen Auflage sind verschiedene Erganzungen und Er- 
weiterungen durchgefithrt worden, die den Umfang des Werkes um etwa 30 Seiten vergroBerten. 
Die Zusatze finden sich in allen Abschnitten des Werkes und gewahrleisten den Anschlu® an 
den gegenwiartigen Stand der Forschung. Sie sollen im einzelnen nicht durchbesprochen werden. 

Es eriibrigt sich wohl, auf die Vorziige dieses bekannten Standardwerkes genauer einzugehen, 
sie sind dem grofen Leserkreis, der sich fiir die Statik der Flachentragwerke interessiert, wohl- 
bekannt und dieses Werk wird sicher viele Jahre hindurch seine fithrende Position im Schrifttum 
behaupten. : G. Heinrich, Wien 


Die Messung von Flammen- und Detonationsgeschwindigkeiten bei der explosiven Zersetzung 
von Azetylen in Rohren. Von P. Hélemann, R. Hasselmann und Grete Dix. (Forschungs- 
berichte des Wirtschafts- und Verkehrsministeriums Nordrhein-Westfalen: Nr. 382.) Mit 
7 Textabb., 26S. Kéln und Opladen: Westdeutscher Verlag. 1957. DM 8-10. 


Im vorliegenden Forschungsbericht werden interessante Messungen mitgeteilt, wobei die 
Methode und die Versuchsanordnung sehr gut erdacht und zusammengestellt wurden. Da die 
Flammenfront leitend ist, kann man diese Eigenschaft dazu beniitzen, die Geschwindigkeit der 
Flamme zu bestimmen. Eilt némlich die Flamme durch ein elektrisches Feld, so wird sich die 
Leitfahigkeit bemerkbar machen. 

Die Verfasser haben es verstanden, nicht leicht zu machende Versuche in sehr geschickt 
erdachten Experimenten auszufiihren und damit einen interessanten Beitrag auf diesem Gebiet 
zu leisten. R. Bruniak, Wien 


Analytische und konstruktive Differentialgeometrie. Von H. Kruppa. Mit 75 Textabb., VII, 
191 S. Wien: Springer-Verlag. 1957. Geb. S 234-—, DM 39:—, $ 9°30, sfr. 39°90. 


Der Verfasser, gegenwartig der hervorragendste Vertreter der international angesehenen 
,,Wiener Schule*‘ der Geometrie, die im Rahmen einer gediegenen Ausbildung von Lehramts- 
anwartern fiir Darstellende Geometrie schon zahlreiche namhafte Mathematiker hervorgebracht 
hat, legt als Ernte einer dreiBigjaéhrigen Lehrtatigkeit an der Technischen Hochschule in Wien 
ein ausgezeichnetes, den Stempel seiner Persdénlichkeit tragendes Lehrbuch der Differential- 
geometrie im dreidimensionalen euklidischen Raum vor, das mathematische Strenge mit geo- 
metrischer Anschauung aufs gliicklichste paart. 

Die erste Halfte des Buches legt unter Verwendung der Vektorrechnung das analytische 
Fundament und bringt den klassischen Bestand der Kurven- und Flachentheorie sowie die wesent- 
lichen Tatsachen aus der Geometrie der ein-, zwei- und dreiparametrigen Strahlmannigfaltig- 
keiten, wobei besonders auf die vom Autor entwickelte Theorie der Regelflachen hingewiesen 
sel, die mit dem Sonderfall der Torsen auch die tibliche Lehre von den Raumkurven umfaBt. 
Die zweite Halfte des Werkes verhilft dann der anschaulichen Seite zu ihrem Recht, und zwar 
durch Betonung der im erweiterten Sinne des Verfassers zu verstehenden ,,konstruktiven‘’ Be- 
trachtungsweise, die nicht bloB auf dem Zeichenbrett, sondern auch unmittelbar im Raum operiert ; 
hier finden natiirlich auch alle mit der Darstellung von geometrischen Gebilden verkniipften 
differentialgeometrischen Fragen ihre Behandlung, tiberdies wird auf manche wichtigen speziellen 
Kurven und Flachen eingegangen, Hervorzuheben ware noch das der Kinematik gewidmete 
Schlu8kapitel, in dem der Verfasser — wie an vielen anderen Stellen — durchaus eigene Wege 
geht. Uberall liegt ihm die Vermeidung des unzulassigen Jonglierens mit dem ,,unendlich Kleinen‘ 
a Herzen, an dessen Stelle saubere Begriffsbildungen und strenge Grenziiberginge zu treten 

aben. 

Das gehaltvolle, duBerst klar geschriebene und mit zahlreichen Figuren geschmiickte Werk 
ist als Lern- wie als Handbuch gleichermaBen geeignet und wird sich einen bleibenden Platz in 
der einschlagigen Literatur erobern. W. Wunderlich, Wien 


Higentiimer, Herausgeber und Verleger: Springer-Verlag in Wien J, Mélkerbastci 5. — Fiir den Inhalt verantwortlich: Prof. 
Dr. Franz Magyar, Wien IV, Technische Hochschule, Karlsplatz 13. — Druck: Manzsche Buchdruckerei, Wien IX, Lust- 
kandlgasse 52 ; 


SPRINGER-VERLAG IN WIEN 


Die Mikrophotographie 


Dr. Kurt Michel 


Aalen/Wiirttemberg 


Mit 550 teils farbigen Textabbildungen. XXXII, 740 Seiten. Gr.-8°. 1957 
(Die wissenschaftliche und angewandte Photographie, Band X) 


Ganzleinen S 840.—, DM 140.—, sfr. 143.40, $33.30 


Das Mikroskop ist schon lange nicht mehr nur das Handwerkszeug einer kleinen Gruppe von 
Forschern. Es ist heute in praktisch jedem Laboratorium ein unentbehrliches Hilfsmittel und 
in sehr vielen Fallen spielt es auch bereits in der industriellen Fertigung eine wichtige Rolle. 
Zum Festhalten der Beobachtungsergebnisse bedient man sich in der Regel der photographischen 
Verfahren. Die Kombination von Mikroskopie und Photographie, die Mikrophotographie, hat 
sich zu einer wichtigen Technik entwickelt, und es wurde schon lange als schwerwiegender 
Mangel empfunden, daB seit beinahe 20 Jahren eine zusammenfassende systematische Darstel- 
lung dieses Gebietes in deutscher Sprache fehlt. Im vorliegenden Band unternimmt es ein seit 
langem bekannter Fachmann, der als Leiter der Abteilung ftir Mikroskopie der Firma Cart ZEISS 
mit den Geraten besonders vertraut ist, eine solche, den neuesten Stand beriicksichtigende 
Darstellung zu geben. Der Band gliedert sich in drei Teile. Im ersten wird ein konzentrierter 
Uberblick iiber die theoretischen Grundlagen der Mikroskopie und Mikrophotographie gegeben, 
ohne deren Kenntnis eine Beherrschung der Praxis immer Stiickwerk bleiben mu8. Der zweite 
Teil behandelt die Gerate und der dritte schlieBlich die mikrophotographischen Methoden. Die 
Ausfiihrungen werden durch eine groBe Zahl zum tberwiegenden Teil fiir das Buch eigens neu 
angefertigter Abbildungen erlautert. 


Induktionsmaschinen. Von Dr.-Ing. W. Sechuisky, Konsultierender Ingenieur bei der All- 
manna Svenska Elektriska Aktiebolaget (ASEA) in Vasteras, Schweden. Mit 370 Textab- 
bildungen. XII, 501 Seiten. Gr.-8°. 1957. Ganzleinen S 534.—, DM 89.—, sfr. 91.10, $ 21.20 


Galvanotechnik. Von Dr. phil. Jean Billiter, emer. a. 0. Professor an der Universitat Wien. 
Zweite, neubearbeitete und erweiterte Auflage. Mit 154 Textabbildungen. VIII, 441 Seiten. 
Gr.-8°. 1957. Ganzleinen S 396.—. DM 66.—, sfr. 67.60, $ 15.70 


_Wirtschaftslehre des Bergbaues. Von Dr. rer. pol. Dipl. Ing. Paul Ippen, frither Leiter 
der Obersten Bergbehérde in Wien. Mit 33 Textabbildungen. VII, 229 Seiten. Gr.-8°. 1957. 
Ganzleinen 8 144.—, DM 24.—, sfr. 24.60, $ 5.70 


Vom Gesetzesstaat zum Richterstaat. Recht als MaB der Macht. — Gedanken iiber den 
demokratischen Rechts- und Sozialstaat. Von René Marcic. XV, 548 Seiten. Gr.-8°. 1957. 
Ganzleinen S 288.—, DM 48.—, sfr. 49.10, $ 11.45 


Das Wahrheitsproblem und die Idee der Semantik. Hine Hinfiihrung in die Theorien 
von A. Tarski und R. Carnap. Von Wolfgang Stegmiiller, Tit. a. o. Professor an der Univer- 
sitét Innsbruck. X, 328 Seiten. Gr.-8°. 1957. Ganzlemen 8 198.—, DM 33.—, sfr. 33.80, $ 7.85 


Zu beziehen durch Ihre Buchhandlung 


SPRINGER-VERLAG IN WIEN 


Soeben erschien das erste Lehr- und Nachschlagebuch iiber Atomreaktoren 
in deutscher Sprache 


Physik und Technik 


der Atomreaktoren 


Von 


Ferdinand Cap 


Tit. ao. Professor an der Universitit Innsbruck 


Mit 100 Textabbildungen. XXIX, 487 Seiten. Gr.-8°. 1957 = 
Ganzleinen 8 384.—, DM 64.—, sfr. 65.50, $15.25 


Seit der Genfer Atomkonferenz ist das Interesse fiir Atomreaktoren in Kreisen der Wissenschaft und Industrie sehr stark 
angewachsen. Gleichzeitig wird es jedoch immer schwieriger, aus der Flut von Einzelarbeiten und aus den sparlichen Buch- 
veréffentlichungen, die bisher in Amerika erschienen sind, ein abgerundetes Bild tiber Theorie, Physik und Technik der 
Atomreaktoren zu bekommen. Professor Cap, der wohl als einer der ersten in Europa schon im Studienjahr 1950/51 an der 
Universitit Innsbruck Vorlesungen iiber Reaktorphysik gehalten hat und der in verschiedenen Reaktorzentren Europas Er- 
fahrungen sammeln konnte, gibt mit diesem Lehrbuch, dem ersten seiner Art in deutscher Sprache, einen Uberblick iiber 
den derzeitigen Stand des theoretischen, experimentellen und technischen Wissens tiber Atomreaktoren. Ausgehend von den 
kernphysikalischen Grundlagen, behandelt der Verfasser Theorie, Bau, Betrieb und Sicherheit solcher Anlagen. Ein besonderer 
Abschnitt bringt in iibersichtlicher Darstellung eine eingehende Beschreibung der bisher entwickelten Reaktortypen. Daran 
schlieBen sich Kapitel iiber die Anwendung der Reaktoren, iiber Strahlungsschutzrecht, Organisationsfragen und Atom-Vélker- 
recht. Die wissenschaftliche Literatur konnte bis zum Ende des Jahres 1956 beriicksichtigt werden. 

Das Buch, das nicht nur als Lehrbuch, sondern auch als Handbuch und Nachschlagewerk gedacht ist, verfiigt iiber ein sehr 
ausfiihrliches Literatur- und Sachwortverzeichnis (deutsch—englisch) und gibt tber alle wichtigen, bisher Lopate No Kern- 


physikalischen Daten Auskunft. Mi 


Soeben erschien: 


Die lsolierstoffe der Elektrotechnik 


Von 


Dr. techn. Wilhelm Oburger 


Dipl.-Ing., Zivilingenieur fiir Elektrotechnik, Wien 


Mit 356 Textabbildungen und 3 Ausschlagtafeln. VII, 542 Seiten. Gr.-8°, 1957. 
Ganzleinen S 660.—, DM 110.—, sfr. 112.60, $26.20 


Durch ihre von den metallischen Baustoffen verschiedenen Higenschaften begrenzen die elektrischen Isolierstoffe die Aus- 
niitzbarkeit elektrischer Maschinen und Gerdte. Die Kenntnis der elektrischen, thermischen, mechanischen und chemischen 
Eigenschaften, ihrer Abhangigkeit von verschiedenen Parametern und ihrer Auswirkung sind fiir den richtigen Hinsatz der 
Tsolierstoffe von groBer Bedeutung. 

Der Verfasser hat sich bemiiht, nach einer kurzen Darstellung des Aufbaues eine Ubersicht iiber die erwahnten Eigenschaften 
zu geben und anschlieSend die gasformigen, fliissigen und festen Isolierstoffe sowie die Lacke mit ihren charakteristischen 
Merkmalen und Higenschaften, Verarbeitungsweisen und Verwendungsmdéglichkeiten zu besprechen. Bei den festen Isolier- 
stoffen ist ein Kapitel jenen auf der Basis von Naturprodukten, ein weiteres den Kunststoffen gewidmet. Die Schilderung 
der Herstellung aller Arten von Spulen schlieBt das Buch ab. 

Da Isolierstoffe praktisch in allen elektrischen Hinrichtungen vorkommen, wendet sich das Buch an den Blektromaschinen- 
und Transformatorenbauer ebenso wie an den Konstrukteur von Hoch- und Niederspannungsschaltgeriten und Stark- ae 
Schwachstromapparaten, zum Teil auch an den Hochfrequenztechniker. 

Die Zusammenfassung aller Daten in zahlreichen tibersichtlichen Tabellen, eine groBe Anzahl von Abbildungen, die auch in 
die Isolierstoffherstellung und -anwendung Hinblick gewiihren, ein umfangreiches Literaturverzeichnis und ein nach allgemeinen 
und Werkstoffbegriffen geteiltes Stichwortverzeichnis dienen der weiteren Unterrichtung des Lesers und ermoglichen ihm, 
rasch eine ausreichende Information tiber bestimmte oder seinen Zwecken entsprechende Isolierstoffe zu erhalten. 
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